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5 Ultrasonic Relaxation Processes in Pure Liquids 
P- 
4 By J. H. ANDREAE anp JOHN LAMB 
- Electrical Engineering Department, Imperial College, London 
Communicated by Willis Jackson; MS. received 27th Fune 1951 
:? 


» SUMMARY. The assumptions made by Hall and Herzfeld in previous theories of 
a absorption are interpreted on the basis of a new and more general analysis. 
formulae given by them are discussed, and it is suggested that the conditions implicit 
jin Herzfeld’s formula should be modified. By this modification a new relation is derived, 
j which provides a satisfactory explanation of the observed relaxation effects in those liquids 
which have been studied. 
‘2 A method of ‘curve fitting’ is employed to compare experimental results with 
_ theoretical requirements. Finally, it is concluded that the ultrasonic behaviour of acetic 
acid and propionic acid is governed by the excitation of bending vibrations of carbon— 
hydrogen bonds. Measurements in carbon disulphide indicate that in this case also a 
_ bending vibration is responsible. Some incomplete data for toluene suggest that here 
a stretching vibration of the carbon—hydrogen bonds in the methyl] group is involved. 


§1. INTRODUCTION 

CONSIDERABLE amount of experimental data has been published on the 
absorption and velocity of ultrasonic waves in pure liquids. In most 
cases, however, measurements have been taken at single frequencies or 

temperatures giving little systematic information. A satisfactory correlation of 
| €xperimental results with existing theory requires a knowledge of the frequency 
and temperature dependence of the absorption and velocity over a complete 
Fae of dispersion: only in the case of acetic and propionic acids (Lamb and 


i A 


... 


i 
' 


Pinkerton 1949, Lamb and Huddart 1950) has this been provided as yet. 


It was shown that the measured dispersion in acetic and propionic acids 
ean be interpreted on the basis of a relaxation of the transfer of energy between 
two postulated states; but in calculating the difference of free energy between 
| these two states a large entropy factor had to be invoked, which is difficult to 
explain. It is well known that the carboxylic acids exhibit association as double 
| molecules and so it was natural to suspect (Spakowski 1938) that the equilibrium 
| between the monomeric and dimeric states might be responsible for the 
| telaxation. An examination of the results on this assumption does not, however, 
| appear to account for the entropy and, further, there is evidence of dispersions 
an ethyl and methyl acetates (Biquard 1936, Beyer and Smith 1946) and in 
€thyl formate (Liebermann 1949), where double molecules do not occur. 
Hall (1948) has calculated the effect of the relaxation of an isothermal change 
in the structure of water, but it is questionable whether the stringent isothermal 
) Zestriction is applicable in the general case. It seems more likely that a process 
__—Ss« PROC. PHYS. SOC. LXIv, I2—B 3¥ 
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allowing for a change of specific heat together with a change of volume would 
be responsible for the observed relaxations. The following general theoretical 
treatment was prompted by recent analyses given by Kneser (1950) and 
Hall (1948). 
§2. THEORETICAL ANALYSIS 
2.1. Derivation of the General Equation 

It will be assumed that the molecules of the liquid can exist in either of two 
states, designated 1 and 2, and that the transition of molecules between these 
two states can be described by the first order reaction equation: ) 


| 

| 
= == R37, [x-=]. af) Thoth & © ieee (1) : 
n, and m, are the numbers of molecules per mole in states 1 and 2, respectively; 
ky. and ky; are the appropriate rate constants and K( =Aj:/k.;) is the equilibrium | 
constant. The total number of molecules per mole, N(=,+ 7 ), is unaffected | 
by the transitions. In the presence of a sound* wave of frequency w/2z7, 
equation (1) becomes jw dn, =k2,n,{dK — d(nz/n,)}, from which is derived the - 


relation NK d(In K) 
a= (14K) l+jer’ 


where 7 =1/(Ry5 + Rp). 
Now the effective adiabatic compressibility 8 at frequency w/2z is governed | , 
by the above reaction and can be separated into two parts as follows: 


p= = (20) == (h(i) _ a 
Taos V\0n,/,\OP]3 V\ 0x] n,\OP/5 
where P, V, S have their sae thermodynamic significance and the quantity a - 


is defined in terms of that contribution to the compressibility 8,, which is not / 
affected by transitions between the two states: thus 


ae Res 
oie V Ox Ns oP s. 


B.=B-Ba=~ 7(5-) (33) a oe 


is called the ‘relaxing component of the compressibility’, and by virtue of fi 
equation (2) will be shown later to have the form 


B,=Bi/(1+jor) =(Bo—Bao)/((L+jor), sean (4) 
where fy is the static adiabatic compressibility. 
A complex velocity of the sound wave V, is conveniently defined as 


Pi onl ber el .aC a 
mea 2 (1%), (3) | 


c WwW 


= 


where c is the phase velocity and « the (pressure) attenuation coefficient. The} 
absorption per wavelength y is given by 
i} 1 


B=2raclw=—7Im(e/V2),  ~— ~ wav (6) Ii 


* The word ‘sound’ is frequently used in the wider sense embracing both audible and|) 
ultrasonic waves. 


+ Kneser (1950) appears to make the error of putting 


ree ohn | 1) (72) : 1 /o Ong 
Po— Ban (+ s \aP i Of BiaiPs Sis (=), 2). peat 
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provided («c/w)?<1. It follows from the above, if p is the density, that 
1/V*=pB=p{B..+Bi/(1+jor)} whence p= {w7/(1+w*7)}7B,/(Bo8..)"”, and 


7=7(8../Bo)"*. When w7=1, has a maximum value 


=47BM(BoBo)M@=SB ly = we eee (7) 


if 8;<fy. The associated fe Cate in phase velocity is given by 
(Cao — €o)/€o = (Bo/ Bao)? —1= py, /77 
where c?=1/pBo, Ca? =1/pBo 
It is now required to derive an expression for f,/8, to account for the 


measured values of ,,. The factor (0”,/0P), appearing in f, (equation (3)) is 
obtained from equation (2) as 


vse), =| Cae), * Car) Ge) lereepreiee 


Using the Van’t Hoff equation [0(InK)/éT]»=AH/RT? and its isothermal 
analogue [0(In K)/0P],;= —AV/RT the previous relation becomes 


ey Fa od : 
N\OP/]s— RT RT?\0P/s_|{1+exp(—AH/RT?)}1+jwr ~ "°° (3) 
where K=exp (—AH/RT) 

and AH =AER+P*AV=ClAT AAP PAV) ieee (9) 


is the heat of reaction per mole for the transition from state 1 to state 2; 
AE is the change of internal energy, at constant volume, and AV is the isothermal 
volume change accompanying the reaction; the internal pressure P, is defined by 
P,+ P= P*=T(S/oV)r. 

The heat of reaction for the transition can be expressed in terms of an entropy 
change dS where 

as 0s oP OV dny 
TdS= (7), dV + (55) 4P= T (Fr), dV —T (57), dP=AH — N° 


((QV\ __ aH, (av) y(aP 
Hence N (5), = T(OP/OT)s 4c (55),4 (5). state sire (10) 


Putting N(dP/dn,),=AP, and combining equations (3), (8) and (10), 
Pa —poAP | AV eA. =) | exp (—AH/RT) 
Pra | prerere— RT RT? {1+exp (—AH/RT)}*~ 


Finally, rearrangement of this equation gives: 


Bi = aay ae nies P*AV Gi 2m rl 
B= [o-D-wePt sr lpi ar ese a 

AH exp (—AH/RT) ‘ 
where C= R & 7) {1 — exp (—AH/RT)}* aueiaueta ce (12) 


is the relaxing component of the specific heat and y (=C,/C,) is the ratio of the 
static specific heats. 

The formula (12) for C; was obtained on the assumption that there is an 
equilibrium between two states. In certain cases, where the heat of reaction 
may result from the excitation of vibrational energy levels of the molecules, there 
is the possibility that the higher energy states may become important. The 

See 
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maximum contribution which the energy levels of a particular vibrational mode 
can make to the relaxing specific heat C,, is given by the Einstein relation 

AE\? exp (—AE/RT) (13) 
RT} {1—exp(—AE/RT)P 


where AE=hy and » is the fundamental frequency for the mode in question. 

In the case of a degenerate vibration the Einstein specific heat has to be multiplied 
by the degeneracy ; when only the lowest two states are considered the reaction 

equation (1) applies with K =n exp (—AE/RT), giving 

AE exp (— AE/RT) i 

= eet tts 1 | 

“ mR (er 7 (T+ nexp (—AE/RT)} is 


C=R( ar 


for the relaxing specific heat. 

In addition to the above there is the possibility that transitions between two | 
postulated states will be accompanied by a change in the number of molecules | 
present. For example, the dissociation of the double molecules of acetic and 
propionic acids into single molecules would be described by the reaction equation 
dnz/dt =ky,(n, + g)n{K —ng?/(ny+ng)n,} where the equilibrium constant 
K=h,»/Ro, as before; m, and n, are the numbers of double and single molecules, 
respectively, and N(=2n,+n,) the total number of single molecules including 
those forming double molecules is constant. It can be shown that the relation 
‘corresponding to equation (2) for this case is 


1, NK?» dink) 
4(1+ K/4)?? 1+jear 
and the contribution to the relaxing specific heat is 
AH exp (—AH/RT) 
7) Tee aH RI ye 
where AH is half the heat of dissociation per mole. 


| 


adn, = 
Gis =-R(Fr 


2.2. Special Cases 


(a) Hall’s isothermal expression for AE=0. The formula for an isothermal | 
process follows from the general equation (11) by putting AE=0 and 
—V-l(0V /0P)p= — V-(AV/AP); this gives 


Ba ae eR (Re) eee | 
isothermal EG RT | : 


Boi Pe yl Gy aT {1+ exp (— AA/RT)P 


where e¢ is the thermal expansion coefficient, and finally 


pe an 1 | 

BENG Wi: 2{1+ cosh (AH/RT)}? | 

which is identical with that given ve Hall (1948). | 
(6) Herzfeld’s expression for AV=0. By assuming that the isothermal | 


compressibility is unaffected by the relaxing process in the liquid, Herzfeld (1941) | 
has derived a formula equivalent to 


Biba CIC, re (17) 


This relation is obtained from the general equation (11) by putting AV =0 and | 
ZT(GP/0T),=TAP/AT. | 


| 
| 
| 
| 
| 
| 
| 
| 
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(c) Alternative expression for AV =0. The assumption 27(0P/0T) = TAP/AT 
upon which Herzfeld’s formula (17) appears to be based seems somewhat 
arbitrary. By analogy with the condition —V—\(0V/0P),)=—V—AV/AP used 
above to give the isothermal relation, it is suggested that in the case where AV =0, 
the conditions TAP/AT=T(@P/dT), or TAP/AT=T(@P/0T), would be more 
appropriate.* ‘The former does not permit a contribution to the compressibility 


{i.e. 8; =0), but the latter gives 
Bi/By = C,/C). ec aL 8) 


In the following sections it will be shown that this equation leads to the best 
agreement with experiment. 

In order to facilitate a comparison of the equations derived in this section 
with experimental results, the chief points of the preceding analysis are 
summarized in Tables 1 and 2. 


SUMMARY OF THEORETICAL ANALYSIS 
General Equation 


iat Bren Bi 4 VAR PFAaV ie 
= ea ~ Lo DvP ar | LG ar — "Ve, 
AH=AE+ P*AV 
Table 1. Special Cases 


Assumptions Resulting equation 
(a) AE=0 
5 aT Ci 
a S = = — ; (5 ji Isothermal relaxation o = ae (16) 
(b) AV=0 
AP eP Isothermal compressibility unaffected Bi Cj 
—_— = = — =(y—1) — 17 
2 AF ma (; ), by relaxing process Bo (y eo 7) 
(c) AV=0 
AP aP 
Fe =() 
@ Tar (sr), Bi 
AP oP Bi C; 
ii = —— (18) 
a AT (=), ee eC 
Table 2. Relaxing Component of Specific Heat 
Assumptions Associated specific heat 
(1) Two-state reaction: general case. C=R AH\?  exp(—AAH/RT) (12) 
AH=heat of reaction ed {1+exp (—AH/RT)} 
(2) Two-state reaction: dissociation of AH 5 pC ATIRT 
AH=half heat of reaction 2 °XP 
(3) Two-state reaction: excitation of 
degenerate vibrational level. C2aR exp (— AE/RT) (14) 
AH=AE=hy TES A E/RT)}" 
n= degeneracy; v=frequency 
(4) Excitation of vibration : all energy Caer ae exp (— AE/RT) (13) 
levels contributing to Cj. Sa {1— exp (— AE/RT)}? 


* See Appendix, however, for alternative conditions leading to a similar result. 
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§3. INTERPRETATION OF EXPERIMENTAL RESULTS 
3.1. Curve Fitting 


The theoretical formulae (16), (17) and (18) which have been derived from — 
the general equation (11) under different assumptions are all of the same form — 


Bi Be=OCOW Sve a eee (19) 
where Q is a factor depending on temperature, and has the values y/(y—1)C,, 
(y—1)/C, and 1/C, for the three cases, respectively. The left-hand side of 


equation (19) is obtained from the experimental measurements by means of | 
equation (7) and it will also be dependent upon temperature in a known way. | 
The temperature dependence of C,, however, is controlled by the value of AH | 
(equation (12)), which is unknown. The following graphical procedure has been | 


developed for the purpose of testing the validity of the theoretical relations (16), — 


(17) and (18) against experimental results. 


The problem to be solved is essentially this: the value of £;/8)Q is known 


from experimental data for a range of temperatures, and it is suspected that C; 


must be of the form 
aT) exp (—AH/RT) 


C(m, nN, AH) =mR & {1 +nexp(—AH/RT)}2 Ss siete 


where m can only have the values —1, +1, +2, +3,.... but AH is unknown. 
It can be seen that m=n=1 gives equation (12), m=1, = —1 gives equation (13) 


and m=n gives equation (14). (It is found that the dissociation equation (15) is — 
practically identical with equation (12) in the range of chief interest.) Now, — 
if p=log,)C; is plotted as a function of g=log,) T for an artibrary value AH, — 
of AH (Figure 1), then moving the graph in the p-direction by an amount — 


Figure 1. Example of curve fitting for acetic acid. 
@ Experimental values. 
AH}? exp (AH/RT) : 
az | {1-+exp (—AH/RT))? for AH=4 kcal. per mole. 
Values derived from fitted curve (— — —) are m=2°6 (=antilog,,) 0:41), 
H=2°8 kcal/mole (=4/antilogy) 0-157). 


Full curve represents C. i=R| 


logio@ is equivalent to increasing m to am; similarly shifting the graph in the | 


q-direction by an amount log, is equivalent to multiplying AH, by b. Hence, 
if the value of log;,6;/8yQ is plotted on graph paper as a function of logig f3 
and curves of logy) C,(m,, AH,) are plotted, to the same scale, on a transparent 
celluloid sheet for the values AH,=4kcal/mole; m=n=1 m= —1 =e 
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m=1, n=2; m=1, n=3; etc., then, by sliding the transparent curves over the 
graph paper, positions for coincidence of the experimental (logy) f;/B)Q) and 
theoretical (logy) C,;) curves can be found. These positions will differ from the 
reference locations, where the axes coincide, by amounts log,,) a and logy) 6, which 
will give the required numerical quantities : m =a, AW =bAH, = 46 kcal/mole. The 
value of m, however, must be unity, except when the curves for n=2, 3,4,.... 
are used, and then m and » must be equal. The value of mis, therefore, a criterion 
for deciding whether the experimental results do, in fact, confirm the theory. 
Figure 1 indicates the principle of the method outlined here and shows the 
experimental values of C,(f;/B) for acetic acid coinciding with a portion of 
the curve peli AH Ney (oALIRD) 

.. (er) {1 + exp (— AH/RT)} 
plotted as described above. It can be seen that log,,a=0-41 giving a=m=2:6; 
log,, b= — 0-16 giving AW =bAH, =4/1:44 =2:8 kcal/mole. 


3.2. Acetic and Propionic Acids 


In order to apply the method of curve fitting described in the previous 
section, it is necessary to calculate from the measurements the values of the 
relaxing specific heat over the available temperature range corresponding to the 
three special cases given in Table 1 and expressed by equations (16), (17) and (18) 
respectively. Equation (17) given by Herzfeld is, however, an approximation 
for the condition C,<C,, and in what follows the exact equation 

B/Po=(y—- YE KC,— C;) or Cj= Cl viel) Og) ij ee es = « (17 a) 
will be used. 

The values calculated from the experimental data for acetic and propionic 
acids are given in Table 3, those of y=C,,/C,, being obtained from the formula 
y—1=«*c?T/C, where c is the measured velocity of sound and the values of C, 
and the expansion coefficient < are taken from the International Critical 'Tables 
(1928). 

The results given in Table 3 are now plotted in the manner described in §3.1 
and the method of curve fitting is employed to determine consistent values of 
m, n and AH, appearing in equation (20). Taking, first, »=1, the values of 
m and AH derived for each special case (‘Table 1) are given in Table 4. 


Table 3 
Cp 2/m (y= 1) Cy 24m 
OG. (Ct = ESE SE! C= ——————— j= Gy 
BOC) Hm (cal/mole/°c.) z ey ee “14 da(y—1/um) ; P 
Acetic acid (eqn. (16)) (eqn. (17 a)) (eqn. (18)) 
20 0-0510 29°25 1-214 0-168 3-855 0-951 
30 0-0545 29-80 1-208 0-178 4-255 1-034 
40 0-0583 30°38 1-202 0-190 4-725 el 32 
50 = 0:-0610 30-92 1:196 0-197 Spi 118} 1-202 
60 0-0641 31°50 1-190 0-205 5:570 1-287 
Propionic acid 
Mee 8 0-0107 © 32:97 1-272 0-0481 0-806 0-225 
21 0-0134 35°11 1:265 0:0593 1-033 0:283 
31 0-0154 33-20 1-260 0:0673 1-207 0:326 
41 0-0178 SocoL 1°253 0-0764 1-428 0:378 


51 0:0201 33°42 1-246 0:0846 1-652 0-428 
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Table 4 
AH\2 exp (—AH/RT) AH=—AE+P* 
ee =) ea = + PFAV 
nee ie 7 (i+ exp (—AH/RT)}* 
Case (a) \E=0 Case (6) AV=0 Case (c) AV=0 
__ 24m ¥~1 yp id Gi re eee Cc 
Preto + Cr 1+432(y—1)/Um a 
: 15 2°6 
Acetic acid Ad Ba 3-2 2:8 
Res > iar m 0:33 13 2°3 
Propionic acid AH 345 4-2 3:8 


According to the previous analysis the calculated value of m in Table 4 should 
be unity in case (a), and where it is other than unity in cases (4) and (c) agreement 
is only possible if a degeneracy be assumed. It has already been mentioned 
that the value of m in equation (20) is a criterion for deciding which theoretical 
approach is applicable. In the process of curve fitting the value of m which is 
obtained must equal unity or the assumed value of m, but in case (4) it is found 
that taking a larger value of m only increases m. However, the values of m for 
both liquids under case (c) suggest that the excitation of a vibrational energy 
level may account for the heat of reaction, AH =AE; in such a case a degeneracy 
of the order of three would be consistent with a vibration of carbon—hydrogen 
bonds. The closest agreement between the derived value of m and the assumed 
value of n is obtained by taking »=3 for both acetic and propionic acids: the 
values then obtained are given in Table 5. The correspondence between the 
frequencies seen in Table 5 and the two frequencies for the bending of 
carbon—hy drogen bonds given by Herzberg (1947), namely : 


(1) »,=1,000 cm} (2) m= 1,450 emt 
H H 
ia WA 
—C—H CH 
\ oN 
H H 


shows that these are probably connected with the process responsible for the 
observed behaviour, and further, the degeneracy of three is reasonable. It is 
concluded, therefore, that the measured dispersions in acetic and propionic 
acids can be interpreted on the basis of a relaxation of the transfer of energy 


Table 5 
AE\? _ exp (—AE/RT) 
C= Rie ie eee 2 
ih (el {1+3 exp (—AE/RT)}*’ AE=hy 
n m AE(kcal/mole) Frequency of vibration 
(v cm) 
Acetic acid 3 3-0 2.9 aa 
Propionic acid 3 2-8 4-0 ais 


between the ground state and a vibrational energy level corresponding to bending 
of the C-H bonds. Different bending modes are responsible for the behaviour 
of the two liquids, that for acetic acid being given in diagram (1) above - 
(¥,=1,000cm~'), and that for propionic acid by diagram (2) (v= 1,450 cm). 
Moreover, the analysis given by the authors with assumption (c) in Table 1. 
leads to self-consistent results, which is not so for assumptions (a) and (6d). 
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4 : It has been assumed that the value of C, to be substituted in equation (18): 


is the static specific heat, which may not be strictly correct if other relaxations. 
occur at lower frequencies. However, as can be seen by reference to columns 3 
and 7 of Table 3, it is unlikely that these would contribute a sufficient amount 
to the total specific heat to cause much difference in the values derived above. 


3.3. Carbon Disulphide 
Measurements of the absorption in carbon disulphide are shown i in Figure 2,, 
these being a combination of results obtained in this laboratory (Lamb and. 
) Andreae 1951) at frequencies up to 70 Mc/s. with some by Rapuano (1950) at. 
%) frequencies greater than 75 Mc/s. The curve drawn is based on the authors’ 
) measurements and is of the form 


af? = ETTIA, cad Se Fie aa, eee (21) 


math A=5,385 x10-!’sec?cm , B=428 x 107!" sec?cm-!, f,=72-1 Mc/s. The 
first term of (21) is characteristic of a single seein process and the 
corresponding expression for the absorption per wavelength, due to the process 
Wis w=(Ac/2n7)w7/(1+ wt) where 2n7=1/f,, 7 being the relaxation time and 

} jf. the ‘characteristic frequency’. Hence p,,=4+Acf,=478,/B, from (7). 
} Allowing fora velocity dispersion of 7°%, since p.,,/7 0-07, and taking Huddart’s. 
| (1950) value for the velocity at low frequencies (cy = 1-14 x 10° cm/sec. at 25° c.); 


@ Measured Values 


x1077 
+ Rapuano, interpolated to 25°C. 


6000 


4000 


a/f? (sec? cm=') 


2000 


2 5 lo 20 50 100 200 
Frequency (Mc/s) 


Figure 2. Absorption in carbon disulphide. 


gives the experimental value for pw, as 0-23. In a non-polar liquid such as: 
_ carbon disulphide it is unlikely that there will be any appreciable order among 
» the molecules at temperatures well above the melting point (—108-6°c.), so that 
) the existence of an equilibrium between, say, an ordered and a random 
arrangement of the molecules is discounted; further, the molecules are stable 
so that an equilibrium of dissociation is improbable. ‘The most likely cause of 
» the observed behaviour, therefore, is to be found in an equilibrium between 
_ molecules which have been excited to a higher energy level and those which are 
in the ground state. The relaxation times to be expected from the excitation 
| of rotational energy levels are too small to cause ultrasonic dispersion in the 
neighbourhood of 70Mc/s. (Figure 2). It remains to examine equilibria 
» between molecules excited to higher vibrational levels and those in the ground 
state. 
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There are four normal vibrations of the carbon disulphide molecule, two being 
degenerate: these and the corresponding fundamental frequencies given: by 
Herzberg (1947) are shown in Table 6. Assuming a two-state reaction it 1s possible 
from equation (14) to calculate the minimum relaxing component Cicnin) which 
-each mode can contribute to the specific heat (x =2 for the degenerate vibration 
and n =1 for the other two modes). Likewise the maximum possible contribution 
Cianax) i8 given by equation (13), where it is assumed that all the higher energy | 
-states of each vibration relax at the same rate. Substitution of these values in 
equation (18) gives the minimum and maximum values of yp, associated with! 
-each mode, assuming that the authors’ equation (18) represents the condition 
AV =0 (Table 1, case (c)). The results obtained in this manner are collected 
in Table 6. Comparison of the experimental value of ,,(=0-23) with the results| 


Table 6. Carbon Disulphide 
T=25° c.; Cp=18-37 cal/mole/°c. (International Critical Tables) 


SS 


Mode of vibration 1 2 $i 
‘Degeneracy, n 1 2 1! 
‘Frequency, v(cm~) 657 397 1523( 
hv/RT=AE/RT Spill 1-918 7:361 | 
_K= exp (—AE/RT) 0-042 0-147 0-00064 
NEN exp (—AE/RT) q 
hei —— ) ———_____~__ ¢al/mole/°c. 0:°772 1-282 0-0689' 
Citminy =" (ar) fitnern( ARRIVE Tee 
AE\? exp(—AE/RT) - ' 
Ciman=2R (az) Ti—exp (AER) cal/mole/°c. 0-914 2°951 0-069! 
Am (min) = 40 ianin)| Cp 0-066 0-11 0-0059" 
Hm(max) = SmCiqnaxy Cp 0 ‘078 0 25 0 -O0SS" 
Max. possible velocity dispersion = “™m) 2-507, 8% <F%i 
T 


‘given in the table leads to the conclusion that the degenerate bending vibrations 
(v,=397 cm") is responsible for the observed behaviour, and it would appear’ 
‘that the higher energy states of this mode do contribute to the relaxing specific: 
heat; this conclusion must be regarded as tentative pending the completion of) 
variable temperature measurements by the authors. 


3.4. Toluene 


Recent measurements by Moen (1951) have shown that the absorption int 
toluene is characteristic of a relaxation process. Moen’s results cover the: 
frequency region 150kc/s. to 1-2 Mc/s. and, when taken together  withi 
measurements at higher frequencies obtained in this laboratory by Huddart 
(1950), make it possible to extrapolate the experimental value of [, Corresponding; 
to 27°c.: this is found to be u,,=3-8 x 10-5. By assuming, on the basis of) 
previous success, the participation of a C-H bond stretching vibration (for which! 
Herzberg (1947) gives v=2,960 cm-) and calculating the expected value of ree 
for a degeneracy of three, from equations (14) and (18) there is obtained the 
value p,,~3-2x 10-5. Within the limits set by the accuracy of Moen’s 
measurements and the process of extrapolating the value of y,,, the two values, 
for », show sufficient agreement to strengthen the hypothesis that the 
‘experimental behaviour is due to the relaxation of the transfer of energy to the 
‘stretching vibration of the C-H bonds in the methyl group of toluene. | 


1 


Se NO EE SR 


— 


——_— a 


Ultrasonic Relaxation Processes in Pure Liquids 1031 


It is hoped to extend these measurements in this laboratory and to study 
‘the behaviour at different temperatures. As in the case of carbon disulphide, it is 
-essential to know the temperature dependence of p,, for toluene before any 
confidence can be placed in the suggestions put forward for this liquid. 
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APPENDIX 
It can be shown that the assumptions underlying Hall’s isothermal 
-expression (equation (16)) may be expressed in the form 


AK tee OV. OV 
AT=0; AP = (5), = (SB). ce sees (A 1) 
. : : ah IN GI oP 
while the authors’ assumptions AV =0; AT = (sr), =e (Fr). ane (A2) 


-do not exhibit a corresponding symmetry. However, if the compressibility is 
-expressed in an Tee form 


view BERL) CH). “BBB ot 


OV\ (OT oy 
po-- 7(5r), (5). (3). Motcicrs (A4) 
-then the final cas reduces to 
S AP Vie NaC, 
** Ge eal et ed 
e-[1+ yBoP “ar | [1 STENT Go ae (A5) 
In this equation i assumptions 
AP oP oP 
AV=0; AT (sr), > (57), oie oe! erate’ (A 6) 
‘lead to the result Bilis a =a Cy) Comma eee. 8) ela an (A7) 


‘Comparison with equation (18) shows that this result does not differ much from 
‘the result used in the paper and more extensive experimental data will be 
“necessary to decide between equations (18) and (A 7). 
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ABSTRACT. An attempt is made to define theoretically the conditions for the- 
appearance of cavitation in liquids subjected to alternating pressure changes. It is found 
that cavitation is restricted to a definite range of variations of the following parameters : 
(i) alternating pressure amplitude, (ii) frequency of the pressure wave, (iii) radius of the: 


bubble nucleus, (iv) hydrostatic pressure. Under certain conditions, the change from 
non-cavitating to cavitating conditions is found to be exceedingly sharp, and in these | 


cases it is shown that the threshold for cavitation can be accurately expressed in very 
simple terms. 


§1. INTRODUCTION 

N a previous paper by the authors (1950), a theory is presented which” 

describes the motion of a gas-filled cavitation bubble in a liquid subjected 

to alternating pressure changes. In this earlier work, the differential 
equation giving the motion of the bubble surface was solved for various arbitrary 
values of the parameters representing the size of the cavitation nucleus and the 
frequency and amplitude of the alternating pressure changes in the medium. 
The data obtained from these solutions were used to form an estimate of the 
intensity of cavitation action and its dependence on these parameters. It was 
found that this cavitation intensity, however it is measured, is very dependent 
on the nuclear size and that for any assumed nuclear radius an upper threshold 
in frequency must exist above which the conditions for true cavitation can never 
be realized. An increase in frequency for any given bubble size will result in the 
changes in bubble radius becoming less violent and finally degenerating into 
comparatively gentle oscillations as the frequency is increased above the 
resonant value of that particular size. 

A further series of deductions has now been made from the original theory, 
and we find that cavitation is restricted, not only to a finite range of frequencies w 
and nuclear sizes, radius Ro, but also to a fixed range of variations of hydrostatic’ 
pressure P, and alternating pressure amplitude P. In the present theoretical 
work, an attempt is made to define the conditions (in terms of the four 
parameters w, Ry, Px, Py) under which true cavitation can be said to.occur. Under 
any other conditions it is found that the bubbles undergo forced oscillations, of 
comparatively small amplitude, at the frequency of the impressed pressure field. 
The change from non-cavitating to cavitating conditions may occur abruptly 
or quite slowly. 

Some of the theoretical results obtained are of special importance. For 
example, in the previous work we were forced to accept the possibility that 
indefinitely small bubbles might serve as cavitation nuclei, but we now find a 
definite lower threshold in Ry below which cavitation does not occur, and this. 
threshold is, in general, extremely sharply defined and easily calculated in terms 
of Py. This limitation in R, sets, in turn, an upper limit to all theoretical estimates. 
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of the intensity of cavitation, however this is expressed. It is also found possible 
‘to deduce an absolute threshold in P) for any assumed w and P, below which 
cavitation can never occur whatever the distribution of nuclear sizes; the 
‘importance of this result lies in the fact that R, is the only one of the four 
‘parameters mentioned above which is, in most practical cases, not fixed by the 
-experimental conditions. The theoretical result therefore applies to any 
ordinary “gassy’ liquid such as would normally be used in practice as a medium 
‘for transmitting the pressure wave. A further, rather unexpected, result is the 
‘appearance of a theoretical upper threshold in P), independent of frequency, 
-above which cavitation is suppressed. 


§2. CAVITATION THRESHOLDS 
(1) Thresholds for Ry and w 


The excess pressure due to surface tension demands that the pressure of gas 

in very small bubble nuclei is extremely high, and it becomes difficult to see 
' how such bubbles can exist for long even in a liquid which is saturated with gas 
' under the existing conditions. For this reason, in the theoretical work already 
quoted, a somewhat arbitrary restriction was placed on the smallness of the 
cavitation nuclei considered. However, recent experimental work (for example, 
Blake 1949) suggests that extremely small bubbles with radii of approximately 
‘0-1 » are in fact stabilized in ordinary liquids, although the method of stabilization 
is not known. We have therefore extended our theoretical survey down to much 
‘smaller nuclear radii, and we find that for Ry sufficiently small, cavitation does not 
occur because surface tension forces prevent the nucleus from ever growing 
large enough. 

For example, radius—time curves have been obtained similar to that shown in 
Figure 1 of the previous work, corresponding to w=9x 104, P)»=4x 108; 
P,=10c.c.s. units and a wide range of nuclear sizes extending from 164 to 
0-15. We find that the curve for Ry =0-15 hardly rises above the time axis, 
while for R, =0-4 the maximum radius R,, is only a few per cent less than that 
_ for R, very large. The relation between R,, and Ry is given in Figure 1; 


Rm (MICRONS) 


0 ’ 2 3 0 
Ro (MICRONS) 


Figure 1. R,, (x) and t,, (©) as functions of Ro for Pyp=4X 10°, w=9X 104, P,=10°. 
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R,, occurs at time ¢,, which is also shown plotted as a function of Ro. The curve 
for R,, shows a well-defined threshold in Ry below which the bubbles do not 
expand significantly and therefore cavitation will not occur. In this region the: 


motion is controlled by surface tension, and the changes in radius and radial 


velocity over the cycle are quite small; we can infer from this that the spatial 
variations in fluid pressure round the bubble during its motion will also be- 
insignificant, so that the liquid tension at a great distance from the bubble at any 


instant, P,+P,sinwt, is almost the same as that in the liquid just outside the _ 


bubble surface. Under these simplified conditions we can readily deduce an 
approximate expression for the threshold R, as a function of the pressure 
amplitude. To this end we examine the conditions for the stable equilibrium — 
of the bubble surface. | 
The pressure at any time just outside the bubble is given by 
p=(P, 42S/Ro) (Rol RF —2ZS/Rx..) ys ieee (1) 
where S is the surface tension of the liquid, assuming the gas changes to be 
isothermal. From this equation, it follows that p is a minimum with respect 


to R when dp/dR=0, namely when R=Ry=[3Ro?(P, +2S/Ry)/2S}?. That | 


this value of R makes f, defined by (1), a minimum may be easily verified. The 
alternative expression for p, obtained by neglecting inertia forces, is that for the 
pressure in the liquid well away from the bubble, namely P, + Py sin wt, and 


unless the minimum value of p, determined by (1), is greater numerically than the: - 
minimum of P, +P, sin wt, or Py)—P,, we have an unstable situation where: - 
the gas pressure and surface tension forces cannot be balanced by the pressure: - 
in the liquid just outside the bubble. Inertia terms must now be included and a. - 
rapid expansion of the bubble occurs. On the other hand, if the minimum p * 
allowed by (1) is greater than P)— P, the bubble radius does not vary much from _ 
its initial value, and throughout its oscillations the gas pressure and surface: © 


tension forces are substantially in continuous equilibrium with the applied 
pressure in the liquid. The bubble radius is not stable in the sense that it’ 


fluctuates wildly if (Py— P,)>(P, +2S/Rp)(Ro/Ry)?—2S/Ry>0 or, substituting © 
for Ry, if (Po—Ps)>4S/3Ryp or >(8/9)[3S3/2R,9(P, +2S/R,)]}". If Py, is: 4 
regarded as constant, then for each value of Py greater than P, there isa minimum. © 


size of bubble given by Rp in this inequality for which the bubble is unstable.. 
Smaller bubbles do not satisfy it and are therefore stable. It may be verified 
that for the conditions depicted in Figure 1 R)=0-19, is the threshold value 
for Py)—P, equal to 3 atmospheres. This value of Ry is very close to the 
critical Ry which can be read off from the curves, and from this we infer that, in. 


this case, the simple approximate treatment given above is sufficiently accurate 


for the purpose. However, we cannot assume that this will be so in all cases: 
for lower values of P)— P, the change will be less abrupt. In the general case, a 


reliable estimate of the threshold can only be obtained from the relevant set 


of (R, t) curves. : 
As we have already shown (Noltingk and Neppiras 1950), an upper threshold 


for Ry exists, corresponding to any frequency w, above which true cavitating 


conditions are never realized. ‘This is given approximately by the equation 
pa® Ro? =3y(P,+2S/Ry)—2S/Ro 5 awe (2) 

which corresponds to mechanical resonance, where p is the liquid density 

and y the ratio of specific heats for the gas content of the bubble. As Ry is 
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increased above this value, the bubble executes a complex oscillation while the 
changes in Ry remain comparable with R,; but as R, is still further increased, 
the oscillations become more nearly sinusoidal. 

We therefore see that both for bubbles such that Ry is less than the lower 
threshold and for those for which R, is much greater than that given by 
equation (2) the changes in bubble radius are relatively small. Under these- 
simplified conditions, the equation for the bubble motion reduces to a particularly 
simple form. This can be found by making the substitution R= R,+AR in (1) 
and solving for AR on the assumption that AR/R, is so small that powers higher 
than the first can be ignored. We then obtain 

RoPo sin wt 

R are Ro a B(PAE4 SSR) ee en elie ereleve (3) 

R therefore varies sinusoidally with ¢ with the same period as the impressed 
pressure but 180° out of phase with it. If AR is less than Ry but higher order 
terms than the first power of AR/R, cannot be ignored, we may readily obtain 
an expression for R of the form R=R,)—a, sin wt+a, (sin wt)?. . ., where 
my, a3, ... . are functions of Ry, Py and P,. The (R, 2) curve is no longer 
sinusoidal but is still periodic with period 27/w. 

Summarizing, therefore, we see that, as Ry is increased from indefinitely 
small values, the changes in radius are at first small and the bubble motion is. 
truly sinusoidal but 180° out of phase with the impressed pressure. The (R, #) 
curve becomes distorted as the lower threshold for Ry is approached, and at the 
critical R, (corresponding to the particular pressure amplitude considered), a 
sudden expansion occurs, representing the start of true cavitation. Cavitation. 
conditions continue until Ry reaches and passes the resonant value (corresponding 
to the particular frequency considered) given by (2), after which the bubble 
motion again becomes sinusoidal and is again given by the simple equation (3).. 


(ii) Absolute Thresholds for Py and w 

We see, then, that there must exist a finite range of nuclear radii suitable for 
cavitation. As the frequency increases at any pressure amplitude, the range- 
is reduced by a lowering of the upper limit. After this has been reduced to: 
coincide with the lower threshold, cavitation can never occur, theoretically, 
whatever the distribution of nuclei. This is true for frequencies greater than. 
that given by the equations (2) and 

48/3 Ry =(8/9)[3.S3/2Ro3(P, +25/R,) 2 =K 
or VaR oRe hae OSha mee BM en ens (4) 
We obtain the relation between w and K by eliminating Ry between (2) and (4): 
and making the substitution Ry=4S/3K. 

The threshold frequencies are found to be high even for moderate values. 
of K. For example, corresponding to K =2 atmospheres, we obtain w/27 = about 
10 Me/s., taking P, =1 atmosphere in a water medium. 

(iii) Limiting Cavitation Pressures 

We have previously shown (Noltingk and Neppiras 1950) that the maximum 
liquid pressures and velocities generated during the collapse of cavitation bubbles 
are very dependent on the size of the bubble nucleus, increasing rapidly as Ro 
decreases. According to the present theory, however, there is a limit to the 
smallness of possible cavitation nuclei and this, in turn, sets a theoretical upper 
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limit to the magnitude of the collapse pressures, etc. For example, we have 
previously shown that the maximum liquid pressure during collapse is given by, 
{P4/81(P, +2S/Ry)}(Rm/Ro)® in terms of R,, and Ro, where P is the effective: 
liquid pressure at r= 00 during the collapse. Since R,, is a function of Ro, for 
a given P, and w there will be an optimum R, at which the cavitation pressure isa_ 
maximum. Obviously, this Ry is almost exactly that for which R,,/Rp is a 
maximum, the variation of P with R, being quite slow. If therefore in Figure 1} 
we draw a tangent from the origin to the curve, the R, coordinate of the point 
of contact is the required optimum value. It is very slightly greater than the> 
threshold Ro, as we would expect: | 


(iv) Thresholds for P, and Py 

We turn now to the effect on the bubble motion of changes in the parameters 
P, and Py. It is not immediately obvious how the rapidity of the collapse of ' 
a cavitation bubble will be affected by changes in the standing pressure Py. | 
The bubbles will certainly grow larger if P, is decreased, but the collapse will | 
occur at a lower effective hydrostatic pressure, and these changes affect the : 
‘cavitation intensity in opposite directions. But if the bubbles become very | 
large the time 7 required for their collapse may be so great that the pressure» 
begins to go negative before the completion of the collapse, which will then become 
less violent. This will be true for t,,+7>@1(27+sin 1 P,/P,)). A lower 
threshold must therefore exist for P, below which true cavitating conditions | 
do not obtain owing to the bubble not having time to collapse completely. 
Reduction of Py, below this limit will result in slower changes in the bubble 
radius. ‘This prediction has been verified by obtaining solutions of the equation 
of motion corresponding to P, small compared with P,. Strictly speaking, © 
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Figure 2, R,, and maximum fluid pressure as functions of P, for 
Py=4 x 108, w=9 x 104, Ry=1-6 X 10-4, 

for bubbles which grow so large, it will not be permissible to ignore the effects — 
of liquid vapour evaporated into the bubble during the motion as was done in | 
the derivation of the equation of motion. The graphs of Figure 2 show the 
variation of R,, and p,, (the maximum collapse pressure) with P, for the case 

=9x 104, Pyo=4x 10°, Ry=1-6 x 10-4 c.G.s. units. p,, was calculated on the 
seen he, that the bunbies at maximum radius contain vapour at pressure 
always taken to be 10* dyne/cm? (about 0-75 cm.Hg), the vapour content 
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remaining constant during the collapse. The effective liquid pressure P at which 
the collapse occurs has been estimated by a process of successive approximation, 
P finally being the mean liquid pressure at r= oo over the period ¢,, to f,,+rT. 
‘The broken line indicates the limiting value of P, below which the collapse 
will not be completed during the time that the external liquid pressure remains 
positive. We see that the peak of p,, occurs at a value of P, greater than one 
atmosphere. ‘lhe curve shows both upper and lower thresholds in P, outside 
which cavitation pressures will not be developed during the motion. It must be 
remembered that these curves refer to particular values of the parameters w, Py 
and Ry, so that only very limited conclusions can be drawn from them. 

Solutions of the equation of motion have also been obtained corresponding 
to P,>P,). They resemble, in all respects, the (R,¢) curves for Py less than 
the threshold pressure, as we would expect, for, unless P, is greater than P,, 
P, is always less than the threshold whatever R, may be. When P,/P, is small, 
the equation of motion is, of course, given by (3). As an example, the curves of 
Figure 3 show the variation of R,, and ¢,, with Py for Py =108, w=9 x 104, 
Ry =3-2 x10 c.c.s. units. The sudden expansion of the bubble occurs at 
P,=P,. We see, then, that cavitation can be suppressed either by reducing P) 
to be less than P,, approximately, or by increasing P, to be greater than P), even 
in ordinary ‘gassy’ liquids containing comparatively large bubbles. ‘These 
theoretical conditions correspond to well-known experimental facts. 
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Figure 3. R,, and tf, as functions of Py for Ry=3:2 x 10-4, w=9x 104, P,=10°. 


It was noted above that decreasing P, relative to Py) may result in the collapse 
time 7 becoming so large that the external pressure has become negative at 
t=t,,+7, and this implies a lower threshold in Px. It also implies an upper 
threshold in P, if P, is thought of as constant. For example, the curve of 
Figure 4 gives t,,+7 as a function of P, for the case w=9 x 104, P,=10%, 
R,=3-2 x 10-4 c.c.s. units, the dotted line showing the points where ¢,,+7 
becomes approximately equal to w~(27+sin“! P,/Po). At this point, P)=6°5 
atmospheres. It can readily be seen from the data previously given (Noltingk and 
Neppiras 1950) that this critical value of Py should be independent of frequency 
if R, is well below the resonant value (making R,, very much greater than Ro); 
also, we may note that the critical P) will not be very dependent on KR, if this is 
well inside the cavitation range of nuclear sizes, for R,, changes only slowly 
with R, in this range. 
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Figure 4. ¢,,+-7 as a function of P, for Ry=3-2 x 10-4, w=9 x 104, P, = 10°. 


S37 CO NCU UStON 
We have shown that cavitation can be suppressed by sufficiently large 
alterations, in either direction, in any one of the three parameters P,, P) and Ry 


when the others remain fixed; w also has an upper threshold, definable in terms. 


of P, and Ry. The theory indicates no lower threshold for w and we should 
therefore expect that the cavitation intensity would continue to increase as w was. 
decreased. However, at sufficiently low frequencies the theory will cease to be 
applicable to practical cases. The reason for this is that if the radian frequency w. 
of the face of the generator becomes small enough, liquid will have time to flow 
towards the low pressure regions, so that the liquid tensions necessary for the 
bubble growth will not be developed. In this case, no alternating pressure 
field would be set up corresponding to the motion of the vibrator. 

Sometimes the change from cavitating to non-cavitating conditions will be 
a gradual one and no precise point can be fixed for the threshold. On the other 
hand, the change can be extremely sharp, as with variations in P, for small Ry 
(see Figure 1). In some cases, simple approximate expressions can be deduced! 
for the cavitation limits. Summarizing, we find that cavitation will not occur: 
(i) for Rp<4S/3(Py)— Px) (or Po <P,+4S/3Ry) whatever the frequency, where 
2SRp=3P,Re+6SR"; (it) for w? >[3y(P,+2S/Ro)—2S/Ro]/pRy? (or Ro. 
greater than a root of the cubic pw*R,?=3yP,Ry+2S(3y—1)) approximately, 
although variations in Py will also play some part, the formulae being more 
accurate the smaller P); (iii) for P, > Py approximately, whatever w and R, may 
be; (iv) for P, so small (or Py so large) that the cavitation bubble has no time to 
collapse before the end of the pressure cycle. 

Combining (1) and (11) we obtain an absolute threshold for Po, for any w and 
P,, below which cavitation can never occur whatever the distribution of nuclear 


sizes. 
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ABSTRACT. The theory of the self-magnetic field in high current discharges, originally 
developed by Tonks, has been investigated omitting a special assumption used in the work 
of Tonks. The influence of the parameters describing the discharge conditions has been 
considered in some detail. A comparison is also made with a simple, though more 
approximate, version of the theory due to Thonemann and Cowhig. The two theories are 
shown to be in good agreement. 


Soe UN ROD U Carton 


HE theory of the self-magnetic fields in high current discharges (the 
‘Pinch effect’) has been given by Tonks (1939). ‘This is based on the 
equations of ‘Tonks and Allis (1937) governing the motion of electrons 

and positively charged ions under the action of electric and magnetic fields. 

The equations were used in a cylindrical case in which the main current is 
in the x direction and a radial current flows to the walls. The theory is, however, 
limited by a special assumption: the walls are assumed to be conducting and so 
highly charged that there is a radial positive current but no electronic current. 

We shall take, in this calculation, the more usual condition of a non-conducting 
wall. In this case, we should expect, in a stationary state, that the total current 
to the walls would be zero, i.e. the positive and negative currents cancel. We 
shall have to make some assumptions, in particular that the difference in density 
of the positive and negative charges is negligible compared with the magnitude 
of this density. Except in a certain limited region, this assumption will be 
shown to be consistent with the results of the calculations. Further, the charge 
density is taken to be small compared with that of the neutral molecules. 

The basic theory is that due to Tonks (1939). 


§2. BASIC EQUATIONS 
As indicated above, we consider a cylindrical case, the main current being 
in the x direction, with a fixed electrostatic field —@V/dx in this direction, 
i.e. OV dx is independent of 7, where r is the radial cylindrical coordinate in the 
plane perpendicular to the x direction. 
The equations for the components of the velocity of the electrons u, and of 
the positives u, (Tonks 1939) are 


(uj (ey 41 On. mu 

nor D, lef ar in, Or Bé ee eeee (1) 
Co peepee 16 0 ONG 

Sona atgreETi yo (2) 
(Ux)e dn , 1 On, 
ae £+ 8 aeoee 1) eee (3) 
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| Ye. gh. 4 3D, oH ay) 
with ee re laces peiy =) n= Ga me’ | 
s e OV _ Om =i) ole opaie tel (4) 
SET te bal ee 
D,=hC,; Dy=4lyC) ; 


where 2,=number of electrons per cm’, m,=number of positives per cm°, 
e=the electronic charge, m=mass of the slcouon c=velocity of light, C, = mean ; 
velocity of the electrons, C,, =mean velocity of the positives, /,=mean free path : 
of the electrons, /, =mean ftée path of the positives, k= Boltzmann’ s constant, 
T,=absolute temperature of the electrons, 7,,=absolute temperature of the | | 
positives, H = magnetic field. 

The magnetic field H =H, where ¢ is the cylindrical angle coordinate. From _ 
the symmetry of the problem, H, is the only component of the magnetic field — 
different from zero; also the densities n, and m, and the potential V, as well — 
as H, are functions of r only. The functions y, and y, are functions of h, given by 
Tonks (1939). 

The magnetic field is related to (u,), by the expression 


1 EH) tau.) = — PEED, | agp (G+ 7 FE) |. oe (5) 


A further relation which we shall need is the continuity equation 
d 
og {Ng 4p )o® } = TAs = age eae eee (6) 


where A is a parameter which allows for the multiplication by collision. 

In the solution of the equations we make the approximation y,=1=y). — 
‘This is a rough approximation, but it will be shown below that the results are : 
not significantly different when y, and y, are introduced as functions of h. 
Further, we use the fact that 7,<T,. It will be shown later that the results | 
-are largely independent of the ratio T ral ees : 

Taking n,=n,=n and (u,),)= Gare, we obtain u, by eliminating dy/dr © 
from (1) and (2) a DD, eee 1 z ; | 

1 cD eh Ds Ee 35 rs — Bé a ert (7) 


Using the continuity equation (6) to eliminate u, we obtain with T,<T, 


Ni f Te Es 1 dn |. 1 dn dh 
a@l(>t+>(1+h))=--—+- 
T, ee PDs ) rae ae au sia esies a | 
th  €& dh : | 
Fitton ndel Siomoos (8) | 
| 
where A=T.D,/T,D:... ial ee (9) 
The equations for the magnetic field (5) in terms of h (see equation (4)) lead to ) 
the form | 
led. 1 dn e 
ain (Ar) = yn Eas (q+ a a) | SMOG bo (10) 


with =-12neD2/mX{C Yaa (11) 


. 2. Re 
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Eliminating » with the help of (7) and (10), we obtain 


ee én énh? Hee oe i 
y dr T+ * +1 +40 +e] + TPAC SAY ae te) 


aad (8) and (12) are then the basic equations to be solved for and h as a 
function of r. 

These equations are solved by numerical methods. The boundary conditions 
at r=( can be found from the limiting forms of the equations. These can be 
seen to be consistent (for small 7) if 

n =n(1 —6r?), Se ete ie aul Heo (13) 
where f=tyn€é 


4b+2f€= a (7 Be a+ i) = 


and 7m, is the value of the charge density at r=0, of either the positive or the 
negative charge. The constant b is determined from the values of the parameters 
:, és i r, yNo- 


$3) DTHE CASE OF A=0 
For the purposes of numerical calculation, the gas was taken to be hydrogen, 
and its density 1016 molecules per cm’. The electron temperature was taken to 
correspond to 20ev. The value of y is then —3-9x10- c.c.s. units. All 
other parameters were subject to some examination, the magnitudes of these 
parameters being given below (Table 1). 


Table 1. Values of Parameters used 


yng —3-90 —39-0 —390 
é 4-0; 8-0; 48-0 0-5; 1-125*; 2-0; 4:5 0-2; 0-5; 1-0 
ine 160 160, 12* 160 


The values of ym, are obtained by choosing m, to have values of 101%, 1014 
and 10% molecules per cm*. The set of values of ym) does, however, allow one 
to see the effect of various electron temperatures and of the density of the neutral 
gas (since y depends on D.?). "The numerical calculation depends only on n/n. 
The variation of m and of A with 7T,/T, was examined in one case €=1-125,. 


_ T./T,,=12, indicated in Table 1 by an asterisk. ‘These calculations show a 


negligible change in n(r) and h(r); these are in fact too small to be shown in 
diagrams on the scale of Figures 1 and 2. It was, therefore, not thought necessary 
to extend the investigation of this parameter to other cases. 

The n(r) and A(r) curves based on the parameters given in Table 1 are given 
in Figures 1 and 2. The n(r) curves are expressed in the dimensionless form 
n|ny. ‘The h(r) values are negative and can be converted to magnetic fields (in 
gauss) by multiplying —A(r) by 26:3. 

A study of these curves for different values of ym) shows further that both 
sets cannot simultaneously be reduced to one another by some scaling factor. 
The general features are as follows: (a) For a fixed value of yn) the maximum 
value of h is independent of the value of €. (6) An increase in the value of € (at 


_ fixed yn) makes the ‘pinch’ more marked, by which we mean the value of r at 
_ which n(r) falls to its half value; this point is emphasized in Table 2. (c) A - 
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characteristic feature of all the n(r) curves, which is apparent in the detailed 
calculations and which can also be seen in the curves, is the increase in the 
absolute value of (1/n)dn/dr followed by a slow decrease. Since (1 /n)dn/dr is 


negative, the curves have the form of a decreasing exponential (as a function — 


of r), after a certain value of r. (d) The h(r) curves rise linearly at first, then 
more slowly to a maximum value, after which the decrease is gradual. This 
maximum value occurs at practically the same value of 7 as the maximum value 


of | (1/n)dn/dr|. 


0 0:2 0-4 06 0-8 1-0 0 0:2 6-4 06 08 1-0 
. r (cm) r (cm) 
Figure 1. The radial variation of the Figure 2. The variation of the magnetic 
electronic density. field as a function of radius. 


a: €=0°5; 6: €=1:1255 “es €=2-03ad 3 €=4-5. Invall cases yx5— — 39-0 and Tj Tp = 160. 


The values of |h,,,,| for yap = —390-0, —39-0 and —3-9 are 14-2, 4:5 and 
1-4 respectively. The position of |,,,,| as a function of 7, for the various € 
values, can be obtained with fair accuracy by multiplying the r values in Table 2 
by 1-60. 


Table 2. ‘ Half value’ Points in n(r) Curves. 


yNo —390-0 SSO SR, 
é 1 ORO 250 aan0-20) 38) 2-08 125° 10-50 48 8 4 
xy (cm.) 0:090 0-175 0:455 0:065 0-145 0-330 0-560 0-019. 0-11650-222 


We consider next the validity of two of the assumptions made in the 
calculation: (i) It is assumed that one can neglect the difference in density of 
the positively charged particles and of the electrons. This point can be checked 
by calculating dy/dr (i.e. the electrostatic field) from (10), obtaining its derivative 
numerically and then inserting in Poisson’s equation. In most cases, the value 
of dn/dr was found to be too small to be calculated with any accuracy; but the 
upper limit of the difference in charge density is small enough to justify the 
assumption amply. (ii) ‘To simplify the numerical work, it was assumed that 
ya and y;, in (4) could be replaced by unity. This assumption has been examined 
in a typical case €=2, T,/T,=160. The calculations for m(r) and h(r) were 
carried through again, using the curves for y,(h), y,(h) given by Tonks (1939). 
Other curves for dy,/dh and dy,/dh, which were needed, were constructed from 


these y,, 7, functions. A comparison of these calculations, with the previous ones | 


Ving ae 
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is made in Figure 3. _ It will be seen that the assumption made does not materially 
change the conclusions drawn from the results. The change in the curves is 


small, for example, compared with the effect of changing € from 2-0 to 4:5 or 
from 2-0 to 1-125. 


(Q) < S 
T 0 2 
oo | | ! | { eel Laie 
0 001 0-02 003 0:04 0:05 006 007 008 009 
MH (B) - 
ies 7 1-0 (€) 
= 
— 
S oe peer (ee BG 2 SiGe 2) 
= 0 005 010 O15 0:20 0-25 030 
2 = Hl 
3 S rr (d) 
S 3 7 40 
x 
~ Ss 
FS aS 09 es SS BOT 2 
2 0 Oo 02 03 04 05 O06 
a) 
a Hl 
x 
g 
~ T 00 ; o 2) 
1 0 09 fs | 1 L | Neal 
0 0:08 016 024 0 O01 02 03 04 O05 O06 O7 08 09 
7 (cm) r (cm) 
Figure 3. The radial variation of the Figure 4. The radial variation of 
electronic density (I) and magnetic U=(~—u,)/D,.€. a: €=4:5; b: €=2:0; 
field (II): comparison of the approx- c: €=2-°0; d: €=1:125; e: €=0°5. 


imate with the more accurate calcu- 
lation. In curve a, yg=1=yp, in 


curve b, y, and y, are treated as 


functions of h. The parameters are 
€=2:°0, ynu= — 39-0. 


Inc, y, and y, are treated as functions 
of h, in all other cases y,=1=yp. 
yNy = — 39-0 in all cases. 


We can further obtain from the n(r), h(r) curves the value of the electron 
velocity (u,), in the longitudinal direction. This is best done in the form 
—u,/D.€ which has the value unity at r=0. This quantity was calculated for 
all the cases quoted in Table 1. It was found to be nearly independent of 7, the 
maximum variation being of the order of 10°. A typical set of curves is shown 
in Figure 4. That the conclusion arrived at is independent of the assumption 
y¥,=1=y, is shown in Figure 4(c). 


§4. THE CASE OF A+40 

In the case AO the solutions of equations (8) and (12) do not show the same 
character as for A=0. This is illustrated in Figure 5 where the n(r) curve is 
plotted (curve I). ‘The parameters here are = 1-125, B=43-2 (cf. equation (14)), 
T.[T,=160. -The value of B is deduced from an ionization cross section of 
0-58 x 10-16 cm? (which is about one-eighth of the collision cross section). 

It will be seen that the n(r) curve will cross the r axis (i.e. produce negative 
values of the density). This point has been checked in two independent calcu- 
lations (for one of which we are indebted to Mr. K. Skwirzinski). The 
calculation has also been repeated with the more exact form of the differential 
equations obtained by dropping the approximation y,=1=y,. This alters the 
result only to a small extent; the curve shown is that obtained in this last 
calculation. That the type of curve obtained is independent of the choice of yi 
and of € has also been checked by investigating two special cases. 
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It is clear on physical grounds that the solution we have obtained must differ 
from the correct solution in the region of small n/n, values. ‘The reason for this 
can be found by examining the value of dy/dr (i.€. the electrostatic field) which 
was done using the most accurate of the numerical calculations. This-shows. 
(as do the less accurate calculations) that the electrostatic field rises sharply at 
small values of m/m) and would tend to infinity as approaches zero, In this 
case the assumption m, =”, breaks down; an examination of the original equations 
shows that the solutions would be markedly changed when x, and n, differ 
appreciably. , 

The numerical calculation for this case proved to be so complicated that it 
has not been thought worth while to show that the solution would yield values. 
of 2,(r), ,(r) which are always positive. There is, however, no reason to suppose 
that one would not obtain a tailing off of the ,, 7, curves in the region r =0-25 cm. 
in Figure 5, in a manner similar to that found for A=0. 
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Figure 5. The radial variation of the electronic Figure 6. Comparison of the radial density 
density (I) and U=(—u,)/D,€ (II) when variation with that for the function 
A40. ymy=—39-0, €=1-125, B=43-2, n=n(1-+a9r?)—*. Curve a refers to this 
A=0:074. function. In curve b€=2:0, yny= — 39:0, 

A=0; y, and 7, are treated as functions 


of h. 

The A(r) curve shows a normal character, i.e. it rises linearly at first and 
then more slowly, though a maximum value is not reached at the largest value 
of 7 used. 

From this curve one can obtain the longitudinal velocity in the form 
—u,/D.€. This differs from the previous set of such curves (see Figure 4) in 
that the values rise sharply with r (Figure 5, curve II). ‘Though the occurrence 
of the extreme values are related to the effect discussed above, a significant part 
of the rise occurs in a region in which equations (8) and (12) should apply. 


§5. COMPARISON WITH OTHER CALCULATIONS 
A solution of the problem with which we are concerned has been given 
(for A=0) by Bennett (1934), on the assumption that the radial distribution 


function can be deduced from the equilibrium state. This treatment, which is 
relativistic, leads to the density function 


n(r) =19(1 + 97?) Rey 5), 
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The same formula has been obtained in a recent paper by Thonemann and Cowhig: 
(1951), using a more physical approach. Their fundamental equation for n(r) 
is obtained by equating the diffusion ‘pressure’ to that produced by the electro- 
static and magnetic fields. The problem is again treated as an equilibrium 
case, with A=0, and the longitudinal equation (3) of Morse and Allis is replaced 
by the assumption that there is a fixed longitudinal velocity independent of the 
radial coordinate r. The constant «) is expressed in terms of the total current, 
the longitudinal velocity, the electronic temperature and the radius of the tube. 
It is further shown that there is good agreement between theory and experiment.. 

Formula (15) has been compared in a particular case with the radial distribution 
calculated from the differential equations given above. The case chosen is that 
of €=2, yn) = —39-0, for which the more accurate numerical integrations were: 
performed (see §3). The value of aw» was fitted at small values of 7, Here (15) 
becomes 2 =m (1 — 2a r?) which is identical with (13) if b=2a 9. The agreement 
obtained can be seen in Figure 6; it is reasonably good, the extreme difference 
being 13%. Further, it has been pointed out (§3), that the longitudinal velocity 
is nearly independent of r, the maximum deviation being of the order of 10°%, 
(see Figure 4(c)). 

A further paper which deals with self-magnetic fields was published recently 
by Schliiter (1951), after the work reported here had been completed. This 
differs from the present work mainly in that the initial equations do not take 
into account the influence of the fields themselves on the velocity distribution. 
These initial equations are hence simpler. With the same condition as above of 
zero radial current, the radial variation of the density is described by a single 
differential equation for A40. A boundary condition is imposed on this equation, 
that the density of electrons should be zero at r=R, where R is taken to be larger 
than the radius of the tube. The n(r) distribution is obtained numerically. 
The above boundary condition is seen to hold automatically in the AA0 case: 
(see §4) where the n(r7) curves always tend to cut the axis. As far as can be judged 


_ from the published curves, the (7) distributions seem to be similar to those given 


above (Figure 5, curve II). 
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ABSTRACT. A train of regularly spaced rectangular pulses of voltage, amplitude- 
-modulated at a frequency f and applied to the control grid of a commercial valve working 
-conventionally, caused space current changes of up to 0:01 amp/cm? from the oxide-coated 
cathode. An unmodulated pulse train, interlacing the first with a delay of tg, sampled the 
-anode current; the samples contained a component at the frequency f,, of low level. The 
component decayed as fg increased, with a time constant of about 1 msec. It was sometimes 
in phase with the component obtained when tg=0, and sometimes in antiphase. The 
-magnitude and phase were often dependent on the history of the valve. Thus an operation 
known to poison the cathode was usually followed, temporarily, by a large antiphase 
component. 


An in-phase component denotes an enhancement of emission following a positive 


-pulse; an antiphase component denotes fatigue. There is evidence that the impedance 
at the interface between cathode core and coating is not responsible for the effects. The 
fatigue may well be largely a surface, rather than a bulk or matrix, effect. It bears some 
-relationship to effects reported elsewhere with pulses involving current densities of 
10 amp/cm?. 


$1) INTRODUCTION 


OOMES (1946) has reported that Ramsay of Bartol has observed two 
types of decay in the pulsed emission of commercial oxide-coated 
cathodes. ‘The first had a time constant of the order of ten microseconds 

but was rarely found in very active cathodes; the second had a time constant 
of the order of milliseconds and was, presumably, found more generally. More 
recently other workers have noted similar effects. Wright (1949) for instance 
has emphasized the disparity which exists between the pulsed emission of oxide- 
‘coated cathodes and the D.c. (continuous) emission; he attributed the disparity 
to a decay effect. He found the more rapid of Ramsay’s decay to be rare in his 
‘commercial cathodes, but the second to be almost always present. He found 
no time dependent effects at ‘low’ current densities (‘low’ would seem to mean 
less than 100 ma/cm?); he was in general, however, observing differences of 
currents of up to ten to one and might well not have concerned himself with 
changes less than, say, 1%. 

Feaster (1949), after having noted that many valves do not show a rapid decay 
of emission during pulses of duration of about 10 ysec., has described how valves 
with anodes coated with sodium chloride show much decay during these pulses 
when the current density at the cathode is about 5 amp/cm?. 

Very recently, Forman and Rouse (1951) have studied the effects caused by 
the presence of small quantities of magnesium in the nickel of cathode sleeves; 
they noted, when pulsing carefully made diodes at current densities of the order 


of lamp/cm”, that the current passed fell appreciably during pulses of duration 
250 psec. 


| 
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The effect of pulsing on the space current of a valve now reported arose when 
Flood and Tillman (1951) attempted to confirm, experimentally, new analyses 
-of the interchannel interference (crosstalk) which arises, in multichannel systems 
-of telephony using time division and amplitude modulation, from the deficiencies 
-of the common transmission path. Early measurements showed an additional 
-source of crosstalk which was traced to a change of, perhaps, 0-1°% in the space 
current, in a valve in the path, following a pulse of additional current. ‘The 
valve had a conventional indirectly heated, oxide cathode. The space current, 
-except during pulsing, was never more than 10 ma. (corresponding to 10—20 ma/cm?) 
_and the pulse never more than 10 ma. The change of current, though small, bore 
some resemblance to those reported by Coomes, Wright, etc. The original 
-apparatus was first modified to minimize the effect, but was later adapted to 
measure it. 

The magnitude of the effect could cause some concern in the design of telephone 
systems using time division and amplitude modulation. Thus it might well be 
-significant in a ten-channel system with a pulse repetition frequency of 8 kc/s. 
(channel spacing 12-5 sec.) and a duration of channel pulse of, say, 8 psec. 
Designers of good quality systems might well have to find ways of avoiding it. 
“The mere collection of data, over a wide range of small valves, is, however, 
probably not worth while. On the other hand, the effect may be of importance 
in the physics of the oxide-coated cathode, particularly as it can be partially 
correlated with the treatment and operation of-the valve and the better known 
properties of the cathode. 

§2. EXPERIMENTAL 
2.1. General 


In the apparatus used to measure the effect, shown in the Figure, the electrodes 
-of the valve under test, e.g. a pentode type CV138 (6F12) or CV1065 (SP61) 
-received conventional voltages. The potential of the first grid V,, was adjusted 


to give a space current J, of 8 ma.; of this current about 2 ma. flowed to the 


screen grid G2 and about 6 ma. (=J,9) to the anode A, whose circuit contained 


a resistive load, R,. G1 received a pulse of voltage every 1,000 usec., whose 


unmodulated amplitude v was about 0-5yv. and whose duration was 100 psec. 


or, for a few experiments, some shorter period). The pulse could be of either 
-polarity. 
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It would have proved very difficult to measure a small difference between. 
the space current immediately before a pulse and that at some time later merely 


by periodically sampling the two currents and comparing them after D.c. amplifi- 
cation. Instead, the train of pulses P,, was modulated by a sinusoidal signal, 


of frequency f,, (usually about 220 c/s.), whose magnitude was adjusted to give 
a depth of modulation m of about 80%. Thus the anode current during a pulse 
was I,)+ gv{1+mcos(2nf,,t+¢)} according to the polarity of the pulse, where: 


g is the mutual conductance of the valve (usually about 5 ma/v.) and ¢ is a constant. 


A selecting gate pulse P, also occurring every 1,000 usec. had a duration adjustable 


from 10-100 psec.; its leading edge could be delayed with respect to that of Py, 


by any time ¢, up to 900 sec. It sampled the voltage across R,, as a measure of 


the anode current J,. 
When ¢,=0, the voltage after sampling was Ry[J,) + gv{1 + mcos (27f,,t + $)}] 


at times common to pulses of both trains Py, and Pp, R,J,9 for any time by 


which a pulse P,, was prolonged after the end of a pulse P,, and zero otherwise. 


The waveform of this voltage is known to contain a component of frequency f,,,. 
which after amplification was measured by the wave analyser as a voltage Vo. 
When ¢,>100pusec., however, the pulse trains Pp and Py, never overlapped, so: 
that the voltage sampled was, nominally, R,J,) and no component V, of frequency ~ 
f,, was expected to appear at the input to the wave analyser. Although the 


inevitable restriction of the bandwidth of the transmission path from A to B 


introduced crosstalk (i.e. gave V, a finite value), component values were so- 
chosen that the crosstalk attenuation V)/V4 was at least 90 db. (30,000 times) 


for all values of ¢, greater than 120 usec.; crosstalk other than that introduced 
by the transmission path made an even smaller contribution to any residual 
value of V,. 

None the less, when #, was, say, 200 ysec., V/V was usually 0-001 to 0-0001, 
corresponding to an attenuation of only 60-80 db.; it far exceeded any expected 


residuals. ‘The attenuation increased with increase of t,, being a few decibels. 


more at ¢, = 500 psec. 


Because the time constant of recovery from the effect noted is of the order 


of 1 msec., it might seem unwise, or at least leading to the loss of some of the 


finer points, to use a repetition frequency f, as high as 1 kc/s. But there are 


reasons why a change to a much lower value of f, is both undesirable and 
unnecessary. First, measurements become more difficult. Second, experiments 
showed that the effect following one pulse is not obliterated by a second pulse; 
VY, at any time arises not only from the modulated pulse immediately preceding 


the selecting gate but also, to.an ever-decreasing extent, from the pulses. 
occurring at earlier times, whose modulations differ progressively. Corrections. 


could be made therefore to the measured value and phase of V4, to deduce with 
sufficient accuracy the value and phase of V, resulting solely from the pulse 
immediately preceding P,. 


When the-correction to phase was made, V4 was either in phase with Vo. 


(and hence V/V could be spoken of as positive) or in antiphase with V, (when 
V,/V, could be spoken of as negative). 


When P,, was a train of positive-going pulses, V/V) was, in general, greater 
if it was negative and smaller if it was positive than when the train was. 


negative-going. 
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‘Triode-connection of the valve gave results usually only slightly different from 
pentode-connection. A few experiments showed V/V, to be substantially 
independent of the depth of modulation and not a rapid function of the 
unmodulated pulse height. It was, however, often dependent on the quiescent 
cathode current J,; V4/Vo, if positive, increased with increase of J,, but if 
negative behaved in no consistent way from valve to valve, sometimes 
becoming more negative with increase of [,, sometimes even becoming less 
negative or positive. 

As described in a later section, V, was often affected by treatment of the valve, 
through temporary change of operating conditions. 


2.2. Results for Valves before any Special Conditioning 


V4/Vo at t;=200 psec. was only occasionally less than 0-0001 for new valves. 
For the type CV 1065 and for valves made in the Thermionics Group of this 
establishment to the specification of the CV 138, VY, was almost always in phase 
with the modulation on the pulse of current passed (i.e. V/V) was positive) 
when the heater voltage V;, was the rated value (6-3 v.); for a minority of them 
V, became in antiphase (i.e. V/V) became negative) when V,, was reduced to 
4-0v. or less. The reversal of phase could lead to values of V/V, as small as 
0-00003 at some intermediate heater voltage. Most of the commercial valves 
type CV 138 tested gave negative values of V,/V5, irrespective of V,,. 

Some valves of the type CV 1065 which had already run for 60,000 hours 
were available for test. ‘They retained an emission comparable with that of 
new samples, but the interface between the coatings and the sleeves of their 
cathodes had developed impedances of some tens of ohms, whereas samples 
from current production had impedances of only an ohm or two. ‘These much 
caged valves gave however amongst the smallest values of V,/V) noted. 


2.3. The Effects of Short Term Conventional and Unconventional 
Usages of the Valves 
‘The more important conditions applied to the valves from time to time were: 
A with V,=V.,.=150v. and J,,=5 ma., often for many hours or days. 
as diodes with J,,=5 ma., often for many hours or days. 
as B but [,,=25 ma. 
for 30 sec. with V,=0, Vy =200V., Ig =20 ma., thereby much exceeding 
the maximum rated dissipation of G2. 
for 30 sec. with V, =V.=200V., J, ~ 16 ma. and J. > 4 ma. 
with V, =V..=0, V;,=4-0v., V.,=7v. for 5 minutes. The initial cathode 

current p fell to a final current q. 
as F, but with V,,=7-0v. and lasting 1 hour. 
with V,=V.=50v. and [,,=25 ma., for 10 minutes. 
as F but V,,=+13 v. 
with V,=V 4 =V.2= —9vV., for many hours. 
as J but for 1 hour. 

Condition A is typical of everyday use and does not normally cause rapid 
changes in the state of the cathode coating. It rarely affected Vy, unless some 
previous treatment, e.g. F, had seriously changed Vy, when A would usually 
restore V, to a more stable value, often that of the valve when new. 
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Conditions B and C, though not typical of everyday use, would hardly be 
expected to cause rapid damage to the cathode coating; indeed because they 
involve no positive ion bombardment of the cathode, they may well undo some 
of the slow damage resulting from prolonged, more conventional, use. Both 
conditions had an influence on V’,, even if they followed a long period of A, C much 
more so than B. Negative values of V,/V, were either much reduced or were 
changed to small positive values. Positive values were less affected; some were, 
however, significantly reduced. 

Operation D often caused temporary softness in the valves, particularly 
those of type CV 1065, and usually some change in V;. The softness gradually 
disappeared within a minute or so, and with it most of any change in Vz. When, 
before D, V/V, was positive, D usually reduced it and occasionally reversed 
its sense; when, before D, V/V, was negative, D often increased it. 

Operation E is less severe on the valve than D; it caused only small changes 
in’ V3. 

Operation F is thought (Metson 1949) to be able to dissociate those compounds. 
present on the control grid whose heats of formation are less than 7ev. The 
dissociation gives rise to gas which finds its way to the cathode; usually the 
emitting surface is poisoned and gq is less than p. The emission, if so reduced, 
can usually be restored substantially to that applying before F, by running the 
valve for a minute or so with conventional voltages applied to its electrodes. 
When q was less than p/4—-and for some valves it was as little as p/20—V,/V, was: 
usually much changed, for a few seconds at least; it was then always negative 
and often large, e.g. 0-003. When q was greater than p/2—and for some valves 
it was as much as 0-9p—positive values of V/V) were usually reduced but 
negative values changed little. If V4/V) had become large and negative 
following F, its recovery to a more normal value could be retarded by reducing 
the voltage of the supply to the heater of the valve when it was in the apparatus 
measuring V,/V5. 

Operation G carries operation F to greater lengths as far as ridding the control 
grid of those compounds mentioned above is concerned. But the higher heater 
voltage applied throughout ensures that recovery, or reactivation, of the cathode 
takes place so rapidly that the reduction of emission resulting from the poisoning, 
which is also going on, is only just measurable. G sometimes affected V/V, 
driving it negative or more negative. It much reduced the effects of any further 
operations F. 

Operation H had little or no effect on V/V 5. 

Operation J resembles F but can dissociate some compounds unaffected by 
F. It affected V4/V, in the same sense as did F, but usually much less so. 

Condition K suppresses all space current. It has been said to be capable 
of accelerating valve ageing, possibly by facilitating the growth of a poorly 
conducting interface between the coating and the sleeve of the cathode. A few 
hours or days of running under condition K had little effect on V4/V>. Time 
did not permit of longer periods, although several months if not years of 
condition K are normally necessary to promote a substantial interface. 

Operation L was to J much as G was to F._ It never resulted in any important 
changes-to V,/Vp. 

The effects of H, J, K and L were not investigated for many valves. 
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§3. DISCUSSION 

Because the effect noted could be either a value of V, in phase with V, or- 
a value in antiphase, it ought perhaps to be looked upon as two effects, both 
being usually present, but one exceeding the other. ‘The argument for simul- 
taneity is strengthened by the results for some valves for which V,/Vy was positive 
though small for values of t, less than about 200 usec., and negative though small 
for greater values of t,, or vice versa. 

Discussion is simplified if it is confined to the effects following positive: 
pulses of cathode current; it enables a direct comparison to be made with the- 
results of the workers using pulses of high current densities. Everything said 
will, however, have its counterpart for negative pulses. 

A positive value of V,/V) can reasonably be called an enhancement of 
emission and a negative value a fatigue. 

There is clearly an important correlation between the effects and the state- 
of the cathode; the enhancement was most significant when the cathode had. 
been so treated that poisoning had been Sanu and fatigue most significant 
when poisoning had been intense. 

Although the effects have time constants of recovery to normal of the same- 
order as those of decay observed by Ramsay, etc., there are at least two reasons 
why they must not be assumed necessarily to have the same origins as those 
resulting from pulsing with current densities of 1 amp/cm? or more. First, 
the experiments with high densities of current deal almost entirely with fatigue. 
Second, the effects now reported seem to be as much a function of the state of the 
cathode before the pulse, brought about by previous treatment, as of any change 
of state effected by the pulse. 

The causes of the effects cannot here be more than the subjects of brief 
speculations. ‘They can hardly arise from the mere presence of an interface 
between the sleeve and the coating of the cathode. The time constant is several 
powers of ten removed from that of interface impedances; moreover the aged 
valves, with much increased impedances, showed no more effect than did new 
valves with very small impedances. 

The causes may lie with the body of the oxide layer, e.g. in changes of 
conductivity. ‘There are, however, such conflicting views as to the mechanism 
of conduction (Loosjes and Vink 1949, Hannay, McNain and White 1949) and 
doubts about the importance of electrolysis and electrolytic conduction (Sproull 
1945, Wright 1950) in controlling emissivity that speculation is dangerous. 
Because increased activation of the cathode can sometimes be achieved by drawing 
a large current through it, the explanation of the enhancement found may merely 
be that of this means of activation. ‘The enhancement cannot, however, be due 
to temperature changes during the pulse, which are of the wrong sign when 
the heat lost (or gained) due to the change of emission during a positive (or 
negative) pulse is taken into account. 

Thus in the experiments described the change of dissipation produces about 
10-8 cal. per positive pulse ita the oxide layer if its resistance is as much as 5 ohms, 
and about 10-7 cal. are required to account for the increased space current 
(assuming a work function of 2 volts). Because the thermal capacity of the 
oxide layer, whose thickness is about 10~? cm., can hardly be less than 10~? cal/°c., 
the mean temperature change is hardly likely to exceed 10-*°c. On the other~ 
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hand the relationship between space current and cathode temperature was, 
under the conditions of test, about 0-02 ma/°c., so that the effects noted—changes 
.of current of 0-0004 ma. or more—would have demanded a temperature change of 
at least 10-? °c. é‘ 

It is tempting to explain the fatigue primarily as a surface effect; if the surface 
has been poisoned, each pulse, if positive, tends to exhaust the cathode, so that 
when it has passed, the emission preceding it is not immediately forthcoming ; 
a negative pulse on the other hand rests the cathode momentarily, enabling it 
to pass more current. The experiments with operation F afford the most striking 
illustration for this explanation. Wright favours surface effects as the cause of 
fatigue during pulses of high current density. 

The electrical effects noted are perhaps merely further manifestations of the 
continuous metamorphosis of the oxide-coated layer. ‘They may, however, be 
-of use in understanding the operation of the layer if only in conjunction with 
other experiments. ‘They may also be of use in determining the state of the 
cathode; if they can be shown to be of use in predicting the useful life of a valve, 
they may have an advantage over the measurement of total emission, which must 
necessarily involve either a cathode temperature or a space current very differen 
from that of normal usage. 


§4. CONCLUSIONS 


The space currents before and after the passage of a positive or negative 
pulse of current of a few milliamperes differ. For some valves the change has 
the same polarity as the pulse, when it can be called an enhancement; otherwise 
it has the opposite polarity and is called a fatigue. The history of the valve’s 
usage often determines the polarity and magnitude of the difference. Two 
‘causes seem necessary to explain the results. It seems hardly wise at this stage 
to link the causes closely with those postulated by workers observing decay of 
emission during pulses of very high current density, although the time constant 
of recovery from the effects noted is much the same as those of the decays 
réported elsewhere. 
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ABSTRACT. Measurements of ionospheric interaction depend on special radio 


» transmissions which are only available on a limited number of occasions, and only on 
» some of these occasions are ionospheric conditions steady enough for direct comparison 
» between experiment and theory. This paper describes the use of averaging methods to 


extract some information from observations made when conditions are not steady. The 


+ results confirm the existence of a ‘ low-frequency anomaly’ observed by Huxley, but 


show no exceptional interaction for a variation of the frequency of the disturbing wave 
around the nominal gyromagnetic resonance frequency. 


; §1. INTRODUCTION 
HE theoretical analysis of Bailey (1937 a, 1938) has been confirmed by 
observations at Birmingham (Huxley et al. 1948, Huxley 1950) and at 
Cambridge (Ratcliffe 1948) in the sense that during favourable conditions 
of observation the relation between transferred modulation and frequency of 


» modulation is of the correct form in magnitude and phase; but there are a number 
» of secondary points which need investigating. Huxley and Ratcliffe reported 
' that no substantial correlation between the interaction and the direction of 
_ propagation relative to the earth’s magnetic field had been found, and the 


experimental evidence on the existence of a gyromagnetic resonance effect is 


+ so varied as to appear inconclusive. It ranges from the qualitative positive result 
teported by Bailey (1937 b), based partly on comparison between the effects of 
* transmitters so far separated in frequency (1,195 kc/s. and 877 kc/s.) that 
+ they might have been effective at different heights in the ionosphere, through 
the results of Cutolo et al. (1946), who found extraordinarily strong interaction 
' for quite low transmitter power, to the completely negative observation at 


Birmingham in 1949 (see below) in which, however, the frequency variation may 


* have been too small to sweep through the resonance curve.* ‘The second point 
| is that in some conditions the intensity of transferred modulation is less at a 


modulation frequency of 60 c/s. than at rather higher frequencies, which is 


' contrary to the theoretical analysis as usually presented. 


The objects of the present paper are therefore as follows: (a) to suggest 


» methods of discriminating between random fluctuations in the observations 
- and significant trends, on occasions when the results would otherwise have 


been discarded on account of excessive variability; (b) to assemble sufficient 


_ evidence of the existence of the ‘low-frequency anomaly’ (the second point 
- above); (c) to record briefly the negative result of the Birmingham observations 


on a gyro-resonance test on 8th July 1949; (d) to suggest the type of modification 
to the theory which might lead to an erplanation of the low-frequency anomaly. 


* Huxley and Ratcliffe (1950) disagree with Bailey’s estimate of the magnitude of effect to be 


! expected near the gyrofrequency. i 


| 
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§2. INTERPRETATION OF SCATTERED OBSERVATIONS 


In all but the most favourable conditions, the depth of transferred 
modulation at a given modulation frequency will vary from minute to minute, 
in a way similar to the general fading of intensity of a signal transmitted via the ; 
ionosphere. (‘Depth’ of transferred modulation here refers to percentage | 
modulation, so that the variation in question is additional to and separate from | 
fading of the carrier.) 

The observations reported in this note were made on transmissions in} 
which the modulation frequency was changed at half-minute intervals, so that 
several minutes were required to cover all frequencies; and the random} 
variations could therefore be comparable with the systematic effect of change | 


of modulation frequency. In spite of this, and other difficulties*, the author; 
considered that the rarity of special test transmissions for ionospheric investigation } 
justified an attempt to extract some information from a necessarily imperfect | 
experiment. On the nights in question the random variations were considerable, , 
so that it is necessary to examine the dispersion of the readings and attempt to: 
find an average value. 

The simplest method is to take a set of observations which appear: 
homogeneous, arrange them in order of magnitude, and plot them at uniform : 
percentage intervals on probability graph paper: if the distribution is gaussian, . 
this plot will be a straight line. Figure 1 shows a plot of this kind of a group : 
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Figure 1. Distribution of observations for a single modulation frequency (120 c/s.), showing the! 
proportion of observations which would be expected to fall above a given value of) 
transferred modulation T if the distribution were gaussian. Curve A: distribution of) 


percentage transferred modulation. Curve B: distribution of logarithm of percentage; 
transferred modulation. 


of readings of transfer of 120 c/s. modulation, the group appearing reasonably | 
homogeneous internally and being obtained during a period when the mean}, 
carrier level was fairly constant. From curve A, of the points plotted as|| 
percentages, it is clear that the lowest four points are not homogeneous with the} 
main group, and that the distribution of the main group is not gaussian. 


Most of the original team of experienced observers had left Birmingham. 


ee en ee en 
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Owing to successful use of decibel units elsewhere, it was decided to re-plot 
the main group of these observations on a logarithmic scale of depth of 
modulation. This has some theoretical support, since the transfer of modulation 
is related to an absorption effect which varies exponentially with length of path 
and absorption coefficient, and therefore a random variation of length and 
absorptivity of the path through the interaction region should give a random 


variation in the logarithmic measure of depth of transferred modulation. The 


central part of this curve (B in Figure 1) is perhaps more linear than the 
corresponding part of A, but both very low and very high values are less frequent 
than they would be in a gaussian distribution. The shortage of very low values 
may be due to the background of noise and interference, but the scarcity of very 
high values appears to be inherent in the distribution. But curve B does suggest 
a stable distribution, from which one could reasonably take an average value for 
comparison with the averages of other groups of observations. The three different 
means which one could take have values as follows: (a) arithmetic mean of all 
observations plotted on curve B, 0-875%; (b) mean of logarithms plotted, 
0-867 °°; (c) 50° mark on curve B, which falls at a value corresponding to 0-871 %. 
The differences between the three means are small, and if the distribution were 
gaussian throughout, the standard deviation would be about +0-05°%, which 
is greater than the difference between (a), (b) and (c); the arithmetic mean is 
therefore sufficiently true. Curve B suggests that if one can collect about 
twenty observations which relate to similar conditions, the mean should be 
reliable to better than +1 part in 10, when the range of individual observations 
is of the order of +40% of the mean. Figure 2 shows a set of (T,,, w) curves 
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Figure 2. Set of average curves for each }-hour period of percentage transferred modulation 
plotted against modulation frequency at values 60, 120 and 300 c/s. 


(T,,=transferred modulation at modulation frequency w/27) for different 


quarter-hour periods during the night. ‘Three modulation frequencies were 
available, 60, 120 and 300 c/s., and each plotted point in Figure 2 for a particular 


) modulation frequency and quarter-hour period represents the average of from 
-2 to 45 readings (one reading only for the very high point 3). The 


significance of the varying shapes of these short arcs of the full (7, w) curves 


_ ean be appreciated by comparing them with the normalized form (Figure 3) 


of the theoretical curve 7',,=7, [1+(w/Gv)?]¥* given by Huxley and Ratcliffe 
(1949). By comparison with Figure 3 it will be seen that a(7\,,@) curve which 
is markedly concave upwards, as the curves of groups 1 and 4 of Figure 2, 
corresponds to w/Gy ~ 2 at 120 c/s. and Ty large, while a curve such as 
group 7 corresponds to w/Gy = 0-5 at 120 c/s. and Ty not so much greater than the 


transferred modulation at 60 c/s. 
4 A-2 
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If 1/72 is plotted against «w? the theoretical formula is a straight line which 
can be extrapolated back to cut the w? axis at a negative value equal to the : 
magnitude of (Gv)? (Huxley et al. 1947); and the data of Figure 2 have been 
plotted in this manner in Figure 4. Only in the case of group 7 do the three : 
points fall on a straight line; in all the other groups the 60 c/s. point lies below } 


to 
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Figure 3. Normalized curve for variation of transferred modulation (relative to value TJ, for zero 
frequency) with ratio of angular frequency « of modulation to collision parameter G7. 
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Figure 4. Extrapolation of relation between 1/Tm* and w? to determine Gp. 


the line through the 120 and 300 c/s. points. In these cases there is no unique! 
line on the basis of internal evidence, but external evidence provides two reasons}; 
for drawing the line through the 120 and 300 c/s. points and ignoring thei! 
60 c/s. point: (a) this yields values of Gy comparable with that given by the ond! 
consistent group (7) and with the values reported elsewhere for observatioril| 
under more favourable conditions, and (b) although the present case, with three) 
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frequencies only, does not indicate which is wrong, observations on other occasions 
have revealed a tendency for the points at higher modulation frequencies to 
lie on a straight line with that at 60 c/s. off the line (see § 3 below). A value of 
Gi for each group has therefore been calculated from the 120 and 300 c/s. readings, 
and these values are plotted in Figure 5 against time for comparison with the 
simultaneous figures of average carrier level. 
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Figure 5. Comparative plot of carrier strength and estimated value of Gy against time. 


The mechanism suggested by the data in Figures 2 and 5 is as follows: 
‘The reflecting layer which was operative at the start is rapidly losing its ionization 
(Gv =2,100 in the first period but only 1,030 in the second) and at about 
0108 hours is almost penetrated by the wanted carrier. This gives high 
attenuation and large transferred modulation with unsteady behaviour. After 
0115 penetration is complete, and a new layer is operative with Gi fluctuating 
around 1,500, until 0230. ‘There is then along period of very disturbed conditions 
until we return to the initial carrier level and to Gy=2,120 at 0350, probably 
due to new ionization following sunrise in the ionosphere. 


§3) THE ‘LOW-FREQUENCY ANOMALY’ 

The scatter in the readings with step-by-step change of modulation 
frequency is almost always so great that one can select a number, perhaps a 
majority, of observations which fit the theoretical curve, while there are others 
which depart widely from it. In such cases the evidence for a departure from the 
theoretical curve must be mainly statistical, in contrast to one experiment 
(Huxley 1950) which established absolutely the existence of the low-frequency 
anomaly in qualitative form, when using simultaneous modulation of the 
interfering carrier by two frequencies. * 

On the lines described in § 2, the whole of the observations for the period 
0120 to 0415 B.s.r. on the morning of 27th July 1949 were plotted as 
distributions on probability graph paper. (The ‘wanted carrier’ was radiated 
by Lisnagarvey on 1,050 kc/s. and the modulated disturbing. signal from 


* Huxley’s report of this experiment was not available to the author when these experiments 
were carried out and analysed. 
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Droitwich on 200 ke/s. with reception at Birmingham. The modulation 
frequencies used were 60, 120, 300, 700, 1,000, 1,250, 1,500 and 2,000 c/s. during 
successive half-minute intervals, and the modulation-frequency output from the 
receiver was measured with a wave analyser having a band width of approximately 
4 c/s.). The observations relating to modulation frequencies of 60, 120, 300, 
700 and 1,250 c/s. are reproduced in Figure 6 (a), where it is obvious that the 
transferred modulation at 60 c/s. is consistently below that at 120 c/s., and not 
greater than that at 300 c/s. It is also noticeable that the readings for 60, 120 and 
300 c/s. divide fairly sharply into two different ‘populations’ while those for 
higher frequencies show more continuous distributions. Finally, the frequencies 
of 1,250, 1,500 and 2,000 c/s. in Figure 6 (b) seem to show a lower limit of 
0-05°%, which probably represents the background of noise and interference. 
(Since the interference level from atmospherics etc. is a variable quantity, and the 
signal—noise ratio depends also on the carrier level which has been eliminated 
from the readings before plotting the percentage modulation data in 
Figure 6, one does not expect an absolutely constant level of background: hence 
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Figure 6. Distribution of percentage transferred modulation. 


a few points fall below the average cut-off level of 0-05°%.) For these modulation 
frequencies the best mean value of 7., is that taken from Figure 6 as the value 
where the line through the points above 0-05°%, modulation (extrapolated if 
necessary) cuts the line marked 50° on the probability scale. 

In a linear plot of 1/7, against w? undue weight is given to readings with 
large values of w and small (and therefore less certain) values of TJ\,, and 
a logarithmic scale has therefore been used, with no attempt to adjust reference 
points, since a preliminary linear plot showed that when the low-frequency 
anomaly was present, a large part of the data (i.e. readings for the higher 
values of w) approximated to a straight line which differs from the theoretical 

Te Saving an intercept on the w axis small and apparently positive. The averaged 
60 C/s sure 7) are in the following form: the averages for ‘all’ the experimental 


Canoe tes s. plus the average for the ‘high’distribution and the average for 
under mo. 
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the ‘low’ distribution for those cases where the total distribution, as plotted 
in Figure 6, appears to include two separate groups of points. Figure 7 shows 
that, on any interpretation of the data, the interaction effect is less at a frequency 
of 60 c/s. than at 120 c/s. 

It was at first suspected that the low-frequency anomaly might be due to some 
instrumental error (Huxley’s double-tone results were not known at the time), 
and confirmation was sought from other receiving stations. Figure 8 is plotted 
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Figure 7, Figure 8. The low-frequency anomaly 
deduced from data supplied by the 
B.B.C. receiving station. 


from data supplied by the B.B.C. Research Department* of reception from 
Tatsfield on 18th May 1949, 0035-0350 B.s.t. ‘The wanted carrier was 
Lisnagarvey, 1,050 ke/s., and the disturbing signal Droitwich, 200 kce/s. 
Individual readings were recorded in decibels relative to 5°% modulation, 
and for any one modulation frequency show a spread of about 6 db.; averaging 
and conversion to percentage modulation gave the smooth curve Figure 8. This 
confirms both the reduced interaction at low frequencies and the incorrect slope 
at higher frequencies which was deduced from Figure 7. 

As far as can be gathered from the small number of tests at present analysed, 
the conditions which lead to a low-frequency anomaly in the (7, w) curve include 
(a) summer night conditions in the ionosphere and (4) low-frequency disturbing 
signal acting on a medium-frequency wanted carrier. A theoretical consideration 
of the anomaly is given in § 5. 


§4. GYRO-RESONANCE TESTS 6TH and 8tH JULY 1949 

In the gyro-resonance tests of 6th and 8th July 1949 the wanted carrier 
was radiated from Westerglen on 767 kc/s., and the disturbing carrier from 
Ottringham, its frequency being varied in 5 kc/s. steps at quarter-hour intervals 
from 1,325 kc/s. to 1,355 kc/s. and back again. ‘The intensity of interaction 
certainly varied considerably during the night, but behaved similarly on 
27th July when the stations used were Lisnagarvey wanted (1,050 kc/s.) and 
Droitwich disturbing (200 kc/s.); these results indicate that the variations are 


* The original readings were lent to the author by the B.B.C., and it must be emphasized that 
the averaging and interpretation is the responsibility of the author only. 
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due to changes in the ionosphere rather than to changes of carrier frequ 
Intensity of 60 c's. transferred modulation, relative carrier level of wanted 
and frequency of disturbing signal are plotted against time in Fi igure 9. 
There is some doubt whether a frequency swing of + 15 ke/s. is sufficient 
include an appreciable part of the gyro-resonance curve if there is in fact 
significant effect of such resonance on the relationship between interai 
intensity and carrier frequency. One might then consider whether the average 
intensity of interaction observed with the 1,340 ke/s. disturbing signal was much 
greater than with other frequencies. But this would be a very doubtful 
comparison since the use of other frequencies might mean interaction in a 
different part of the ionosphere, and the ionospheric conditions are also liable 
to vary from night to night. However, the general level of interaction shown 
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Figure 9. Variation of received carrier strength and transferred modulation during frequency 
variation of the interfering carrier about the gyro-resonant frequency. 


in Figure 9, 1-3°,, is comparable with that observed on other occasions | 
(Huxley et al. 1948). It is therefore concluded that these observations show no | 
positive evidence of enhanced interaction near the gyro-resonant frequency, | 
although they do not exclude the possibility of such resonance in some conditions. . 


ws 


5. POSSIBLE MECHANISM OF LOW-FREQUENCY ANOMALY 

A cause of irregularity at the high-frequency end of the (T,,, w) characteristic : 
is the existence of interaction at two different regions in the ionosphere (Huxley 
1930), with a phase difference of 27l/A, where / is the path difference and A the? 
wavelength of the modulation frequency. This mechanism seems impossible } 
for an anomaly at a frequency as low as 60 c/s., since the path difference would be} 
excessively large. 

The most likely source of a time constant of the order of 5-10 milliseconds } 
(since the low-frequency anomaly shows maximum interaction occurring at) 
100-200 c s, instead of at zero frequency) is in the rate at which the free electrons | 
in the ionosphere absorb and lose the additional energy which is imparted to) 
them by the disturbing wave, and it seemed to the author that there might be an} 


} 


Ionospheric Interaction in Disturbed Conditions 1061 


analogy with the question whether a sound wave would be propagated 
adiabatically or isothermally in a given medium.* 

The usual method of calculating the interaction (see, for example, Huxley 
and Ratcliffe 1949) is to use the approximation that the collision frequency can 
be taken as a constant p during the modulation cycle for the purpose of calculating 
the rate of loss of energy by collision (although it is in fact the variation of v during 
the modulation cycle which causes interaction); write 

(d/dt)O(t) + GAQ(t) — Qo] = P(t) 1 
rate of change of K.E. + dissipation = power absorption ies 
and assume that P(f) is obtained simply by substituting &(t) = &(1+ M cos wf) 
into a steady-state formula P, =f(v, pp) x &p?, where pp is the (angular) carrier 
frequency of the disturbing signal and &) its field strength. 

The general result of replacing v by »(t) in equation (1) will be to reduce the 
intensity of interaction, for, when the electron energy rises during the modulation 
peaks, » will also increase (being proportional to the square root of the total 
agitational energy), and this will increase the dissipation term Gy[Q(t) — Qo], so 
limiting the rise of O to a lower value than it would have had with v remaining 
constant at the mean value. But since v* is proportional to the total agitational 
energy of the electrons, this effect can only be important when the increase due to 
absorption of wave energy is an appreciable fraction of the thermal energy; 
and the evidence provided by the linearity of the relation between transferred 
modulation and power of disturbing wave (Huxley and Ratcliffe 1950) gives 
no positive support to the idea that the wave energy is comparable with the 
thermal energy. Moreover, although the mean wave energy may rise slightly 
at the lower frequencies as it follows the modulation, the author has not been 
able to find any quantitative modification of the quasi-equilibrium treatment 
of the process which will account for an actual reduction in transferred 
modulation at 60 c’s. relative to 120 c/s. 

It therefore seems that the ‘power absorption’ in equation (1) must. be 
regarded as a function of the instantaneous disequilibrium between the field 
strength and the electron energy, and not merely of the equilibrium condition 
which would eventually be attained at a given power level, and that the — 
modulation component of fundamental frequency, as well as the distortion, is 
modified by this. The general experimental evidence indicates that w <Gv when 
w/2z=60 c’s., so that at this frequency the rate of loss of energy from the electron 
should be able to keep pace with the falling part of the modulation cycle. The other 
factor is the rate of absorption of energy during the rising part of the cycle. 
In this connection it should be noted that the approximate formula for power 
absorption (Huxley and Ratcliffe 1950) assumes that the squared collision 
frequency is small compared with p,? already defined. But since Gy is of 
the order of 1 to 2 x 10, while estimates of G range from 1-3 to 2-6 x J0-, it follows 


* The point was expressed in this form in a report of the observations which was circulated to 
the Radio Research Board in August 1949. A similar idea was included in a note (unpublished) 
communicated by L. G. H. Huxley to a discussion in Birmingham in November 1949, where 
he comments on the limitations of the formula Tw=T,[1+(w/Gy)?]-1: “ The formula in 
question is obtained when the collisional frequency v in the factor Gv is replaced by the mean 
value > and treated as a constant in this factor.” This would account for observations showing 
the transferred modulation reduced at frequencies below 100 c/s., and “‘ The formal rigorous 
solution [of the equation of energy balance of the electrons in the field] is 


0-0,-| exp( = si Gr at) |x sh w(t) [ex( |) Grat) | dt-+const. xexp( -|, Gy at). 
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that » is not much less than 10% and therefore comparable with the angular 
frequency of a long-wave transmitter (e.g. 200 kc/s. or 167 ke/s. carrier 
frequency).* As noted at the end of § 3, the low-frequency anomaly seems to be 
dependent on the use of a low-frequency disturbing signal, so that v° may be 


comparable with Pp)’. 
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Heavy Beam Loading in Linear Electron Accelerators. 
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ABSTRACT. 'The influence of the beam current on the length, the obtainable particle 
energy and the efficiency of a linear accelerator is considered. ‘The beam loading is especi- 
ally important if the accelerator waveguide has low losses. The efficiency, defined as the ratio 
of H.F. power transformed into kinetic energy to the total H.F. power fed to the accelerator, 
can have a high value, 70-80%, but the particle energy obtainable drops considerably with 
heavy beam loading. 


§1. INTRODUCTION 
N theoretical investigations of linear accelerators it is common to neglect the 
| effect of the beam loading on the accelerator performance. This is justified 
when only the field configuration is being considered. However, the 
influence of the beam loading on the attenuation, and consequently on the 


* Huxley and Ratcliffe state that “ Under some important conditions the assumptions which 
allow us to replace equation (2) by equation (3) (i.e. v/(v?-+pp2) by v/pp”) are not valid ”’, 
+ Now at Chr. Michelsens Institute, Bergen. 
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obtainable particle energy and on the accelerator length may well be appreciable. 
‘This is especially important in the case of low-loss accelerator tubes and it may 
be necessary, for instance, in dielectric loaded tubes (Shersby—Harvie 1948 b) to 
take the beam loading into account even for fairly small beam currents. Further, 
high beam currents are essential for some applications, e.g. in accelerators for 
medical purposes from which a high x-ray yield is required, and in accelerators 
for electron sterilizing experiments on food, etc. 

In this paper the influence of beam loading will be investigated. In the 
calculations it is assumed that the field configuration is unchanged by the loading, 
an assumption which is valid, as already mentioned, in all practical cases, and 
which has been confirmed by a separate investigation. 

As the electron velocity is close to the velocity of light along most of the 
accelerator, we shall confine ourselves to the case in which waveguide parameters 
are unchanged along the accelerator. The calculations are, therefore, not 
strictly correct for the bunching section, but this will hardly affect the overall 
performance of the accelerator. 


§2. ESSENTIAL ACCELERATOR PARAMETERS 

In this section we shall briefly summarize the main parameters which 
determine the operation of a linear accelerator. | 

The waveguide used as accelerator tube is best characterized by its equivalent 
shunt impedance per unit length and its equivalent series impedance. The shunt 
impedance per unit length is defined as 

Mae eee en ee (1) 
where E£ is the amplitude of the axial electric field on the axis of the waveguide 
and w is the H.F. power dissipation per unit length, when there is no beam 
loading. (Sometimes different definitions are used, but there is always a simple 
relation between them.) Approximate values for 7 have been worked out or 
found experimentally for the most commonly used waveguides (Shersby—Harvie 
1948 a, b, Mullett and Loach 1948, Slater 1948), and therefore 7 will be assumed 
to be known. 

The series impedance is defined as 

cae Tel dst Lag nen We allie ail Soa shee (2) 
where E has the same meaning as earlier and P is the H.F. power flow. € is also 
assumed to be known, as it can be worked out fairly accurately for the waveguides 
in question. Walkinshaw (1948) has calculated it, for instance, in the case of the 
metal corrugated waveguide. 

If an external feedback is used, it is characterized by the feedback ratio, which 
is defined as (Shersby-Harvie and Mullett 1949) 

att PMY, cau Wh lee Aw mld A eee ere (3) 
where P, is the power flow at the input of the accelerator and W, is the power 
delivered by the H.F. source. 

In the calculations the ratio £,/E), where the indices 0 and 1 indicate the value 
at the input and output respectively, will be used as a parameter. If the feedback 
ratio and the losses in the feedback circuit are known this ratio can be worked out. 
If no feedback is used, it is determined by the amount of H.F. power one can 
afford to lose in a matching load, or by practical considerations which limit the 


accelerator length. 
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The degree of beam loading is most conveniently expressed by the 
dimensionless quantity ne cot 


where 7 and E, have the same meaning as earlier, J is the beam current and 
dy is the angle by which the stable point precedes the peak of the wave. The 
physical significance of the quantity m becomes more clear if (4) is rewritten 
as m=E,I cos $9/(E,?/n) from which it is seen that m is the ratio of power 
absorbed by the beam per unit length to waveguide losses per unit length, all 
at the input end of the accelerator. 


§3. THE ACCELERATOR LENGTH 
From the definition of the shunt impedance (1) the loss per unit length is 
w=k?/n. In addition, the beam absorbs an amount of H.F. power 
wy EL c08:P5) Sve) «0 Seal Ue (5) 
per unit length. The total conversion of H.F. power per unit length is then 
w,=w+w,=(1+y7l cos ¢/E)E?/n or, taking (4) into account, 
©. =(l mE Ei. ee ee (6) 


The series impedance is defined by (2) and this, inserted in (6), leads at once 


to the following differential equation for the power flow: 

dP= —(E/n)|[P+m( PoP)? ]dz |. 22 ae ae (7) 
the solution of which is 

2=(2n/€) in[(mE1)/e E/E ee eee (8) 


where z is measured from the input end, where the power flow is Py. The 
length of the accelerator is given by 


EL/2n = In[(m+1)/(m+E,/E)). wee (9) 


It is difficult to build a feedback bridge for feedback ratios much different 
from unity and the two most important cases are, therefore, this case and that in 


which no feedback is used at all. Two examples are shown in Figure 1. In 


curve A, no feedback is used, and it is assumed that 10% of the u.F. power is 
dissipated in the matching load, which gives (E,/E))?=0-1. In curve B, the 
feedback ratio is assumed equal to unity and the losses in the feedback circuit are 


neglected, “This gives (#,/£,)°-=0-5. 
Firstly we notice the fact pointed out by Shersby-Harvie and Mullett (1949), 


that a feedback reduces the accelerator length considerably, and secondly, the 
point which is of main interest in this connection, that the accelerator length 


decreases rapidly with increasing beam loading. 


Often, especially if there is no feedback and the beam loading is small, the 


frequency sensitivity and other practical considerations may limit the 
accelerator length (Shersby-Harvie 1948 a). It may then be convenient to 


rewrite (9) as E,/E,=(m+1)exp(-€L/2n)-m a... (9 a) 


and to substitute this value for £,/E) in the expressions given in the next two 


sections. 
§4. THE PARTICLE ENERGY 
The particle energy at a point z is given by 


eure | E cosy dz, = an beet ee (10) 
0 


Heavy Beam Loading in Linear Electron Accelerators 1065 


where Up is the injection energy of the electrons. If (7) and (2) are taken into 
account, this gives 


U— Uy) =(2n/£)Ey cos $o{1 — (E/E) —m In [(m+1)/(m+E/E,)]}. .....- (11) 
The particle energy at the output is given by the same expression when FE =F. 
_A more convenient expression, however, is obtained by substituting 


este ( Baer) Ve AG rete Timur Sd kame (12) 
which is derived from (2) and (3). For the gained particle energy this gives 


Uy 0 f 
fae =(1+7)*{1 — (E,/Ey) — m In [(m+ 1)/(m+ Ey/Ey)]}. ....(13) 


This expression has been plotted in Figure 2 for the same two cases as earlier. 
It is seen that the obtainable particle energy depends strongly on the beam 
loading when the length is assumed to vary according to Figure 1. Already at 

m=1 the energy is only about half the theoretical value for zero beam loading. 
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Figure 1. L=accelerator length, 7=shunt Figure 2. U,—U,=particle energy gain, 
impedance per unit length (ohm/m.), W,=available H.¥. power, m=degree of 
&=series impedance (ohm/m?), m=degree beam loading. 


of beam loading =7/ cos ¢y/Eo, where I is 
the beam current, E, the maximum 
accelerating field on the axis at the input, 
and ¢, is the angle by which the stable 
point precedes the peak of the wave. 


At first sight it may appear surprising that the curve in Figure 2 for the 
accelerator without feedback always lies above that for the accelerator with 
feedback ratio equal to one. The reason for this is mainly that the length, 
according to Figure 1, is much shorter for the feedback accelerator. However, 
from the definition of m it is clear that if we are interested in a given beam current, 
its value for the accelerator with feedback is (1 +7)~!? times that when no feedback 
is used. This means that, in this case, the accelerator without feedback has higher 
loading than that with feedback, and we shall not necessarily get higher particle 
energy by choosing an accelerator without feedback. ; 

In addition to this comes the fact already pointed out at the end of the previous 
section, that an accelerator with no feedback and small beam loading often cannot 
be made as long as assumed in Figure 2. If, for this reason, L is fixed and not 
E,/E,, the dependence of the particle energy is, of course, smaller than found 
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in the examples shown in Figure 2. In that case it is convenient to insert from 
(9 a) into (13). This gives, in the case of no feedback 


U,—U, ies = \s ‘9 
(Wa[E)?In cos Fy =(1+m)[1—exp (—€L/2n)]—mEL/2n, ...... (13a) 
which shows that U,— U, in that case is a linear function of m, whose slope 
depends on the parameter €L/2y only. If this parameter is much smaller than 
unity the slope is small. From Figure 1 it is noticed that if the length has to be 
fixed smaller than the value given by curve A, €L/2r is usually considerably smaller 


than unity. 
§5. THE EFPEICIENCY 


The efficiency of a linear accelerator is defined differently according to the 
purpose for which the accelerator is designed. Ifas high particle energy as possible 
is required from the accelerator, a definition is chosen which is closely related 
to the shunt impedance per unit length (Shersby-Harvie 1948 a). Sometimes, 
however, it is of interest to see how much of the H.F. power is transformed into 
kinetic energy in the beam, and it may be of importance to make this figure as 
high as possible. In this case it is convenient to adopt the definition of efficiency 
current in engineering practice, viz. the ratio of H.F. power transformed into 
kinetic energy, to the total H.F. power fed to the accelerator. This can be written 
as 


2=1K0, -U) ke See (14) 
Substituting from (13), (4) and (12) we obtain the final expression 
n=2(1+7r)m{1 —(£,/Eo)—m In [(m+1)/(m+ E,/E,)]}. ...... (15) 


In Figure 3, 2 has been plotted for the same two cases as before. For high beam 
loading the efficiency is, as expected, highest when a feedback is used, but the 
difference is small. From (15) it can be shown that tends to unity, as the loading 
increases, if the losses in the feedback are neglected. If no feedback is used 7 tends 
to 1—(£,/E,)*. Again the reason why the difference between the two curves 
in Figure 3 is so small is mainly that the feedback accelerator is so much shorter 
than the accelerator without feedback. If the latter, for practical reasons, has 
to be made shorter than found from Figure 1, its efficiency drops accordingly, 
as seen from (9a) and (15). 


1-0 


Curve A r=0 £,/Ey=10"? 
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Figure 3. n=efficiency = 1(U,—U))/W», m=degree of beam loading. 


The interesting thing to notice is that the efficiency defined in this way 
increases so rapidly with m that there seems to be little difficulty in transforming 
most of the H.F. power into kinetic energy in the beam. If the accelerator is 
designed for the purpose, 70-80%, should not be too optimistic a figure. However, 
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as shown in the previous section, a high efficiency can be obtained only at the 
expense of the particle energy, if the power from the H.F. source cannot be 
correspondingly increased. 

§6. CONCLUSIONS 

The calculations in the previous sections have shown that in all respects the 
beam loading affects the performance of a linear accelerator, usually unfavourably, 
but beneficially for some applications, as shown in §5. It is only if the degree 
of beam loading, defined by equation (4), is considerably smaller than unity that 
theories worked out for zero beam loading apply. To see more closely what this 
implies, a numerical example will be given. 

Let the available H.F. power be 2 Mw., the wavelength 10 cm. and ¢)=35°. If 
a hole diameter of 4 cm. is used a reasonable value for 7 is about 30 MQ/m. With 
this hole diameter it can also be shown that €=3 MQ/m?. This gives 
m= 10 I(1+1r)-1?, which shows that if there is no feedback, m reaches unity for 
a beam current of 100 ma. (in each pulse). If a unit feedback ratio is used, 
m=1 when [=140 ma. 

If we now assume that 7 = 1 and neglect the losses in the feedback, the following 
results can be obtained for an accelerator using the waveguide and H.F. source 
specified above: (i) If the beam current is neglected it is found from Figure 2 
that U, — U) ~ 16 mev. and from Figure 1 that L ~7 m. (ii) If the beam current is 
increased to 50 ma., the same figures have already decreased to U, — Uy = 12 Mev. 
and L>5m. About 30% of the H.F. power is in this case transformed into 
kinetic energy. (iii) If we increase the beam current to about 300 ma., we get 
only U,—U,=5 mev., L=2m., but now a has increased to about 70%. 
The mean current in a linear accelerator can be up to about one thousandth 
of the peak current. From this example, it can therefore be concluded that 
beam currents usually required will not affect the performance of a linear 
accelerator very much, and that the theory worked out for zero beam loading can 
be used. That is true for all waveguides used in accelerators built so far. 
However, if waveguides can be made which have a much higher 7, e.g. 5 to 10 times 
as high as for the all-metal waveguide, the beam loading must be taken into 
consideration even for much smaller currents than those assumed in (il). 

For some purposes, such as in the use of a linear accelerator as a neutron 
source or as an accelerator for sterilizing experiments on food etc., it may be 
essential to use the highest obtainable beam current. In these cases, the loading 
must be taken into account in the design of the accelerator, and the calculations 
carried out here indicate that satisfactory performances can be expected. 
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Currents in a High-Voltage Ion Accelerator Tube 
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Communicated by D. R. Chick; MS, received 13th April 1951, 
and in final form 27th June 1951 


ABSTRACT. It is well known that a large electron current passes through a vacuum tube 
when used for the electrostatic acceleration of positive ions. This paper deals with experi- 
ments which show that many of the electrons are produced by neutral particles bombarding 
the lower end of the tube. The particles are formed in the region of the ion source canal 
and have alow energy due to the extraction potential. It has been found possible to suppress 
all the electrons produced below the bottom of the accelerator tube. 


§1. INTRODUCTION 

HEN describing the performance of vacuum tubes for the electrostatic 

acceleration of ions, many authors have stated that the total current 

passed by the tube is from two to ten times greater than the ion current. 
Craggs (1942) showed that many fast neutral atoms enter the accelerator tube 
through the canal of a glow discharge ion source. However, the positive ion 
focusing efficiency of his accelerator tube (Craggs 1942) was such that a large part 
of the ion current was lost from the beam, therefore, although he concluded that 
most of the unidentified current was due to electrons produced by bombardment 
of the metal vacuum fittings at the base of the tube, he was unable to distinguish 
between secondary emission due to charged and uncharged particles. 

The current in excess of the target current is a useless load on the generator 
supplying the accelerating voltage and the electrons give rise to a high intensity of 
x-rays. Also they are probably instrumental in causing breakdown in the tube. 

Experiments described in this paper have been carried out on a tube with a high 
focusing efficiency and it was therefore possible to show that the excess current 
was mainly due to the bombardment of the base of the tube by neutral particles. 
This was confirmed by measuring the current passing through the tube when the 
ions were deflected back to the top of the tube immediately after entering it from 
the ion source canal. It was found also that nearly all the electrons originating 


below the last electrode could be prevented from entering the bottom of the 
accelerator tube. 


§2. THE ACCELERATOR TUBE AND ION SOURCE 

The accelerator tube consisted of three porcelain cylinders making a total 
length of 54inches. In these there were twenty aluminium accelerating electrodes 
having an inside diameter of 7} inches (Figure 1). Because of the short length of 
the tube, in comparison with the electrode apertures, practically all the particles 
from the ion source canal passed unobstructed through its full length. 

The accelerating voltage applied to the tube was obtained from an electro- 
static generator and could be set at any value up to 500kv. A nearly uniform 
potential gradient along the tube axis was obtained by connecting the accelerating 
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electrodes to suitable tapping points on a bleed resistor. A series of helical 
Perspex sleeves fitted outside the accelerator tube served to support the elements of 
this resistor, and also eliminated corona losses. 

The ion source was similar to that described by Thonemann et al. (1948). The 
extractor voltage used was about 3-5 ky. except when it was varied to check the effect 
of this parameter on the electron current. 'The gas pressure in the discharge could 
be varied from 5 x 10-*mm. Hg to 3 x 10-?mm. Hg without producing much 
change in the pressure in the accelerator tube, while the pressure in the tube could 
be varied from 5 x10-'mm. Hg to 10-*mm. Hg causing little change in the 
pressure within the ion source envelope. The proton content of the ion beam 
varied with time, but it was always between 80% and 90%. 
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Figure 1. Accelerator focus electrodes and secondary electron suppressing electrodes. 


The inside of the top plate of the tube was coated with willemite, so that 
electrons striking the plate produced a visible image which could be seen in a 
mirror placed in the vacuum chamber at the base of the tube. 


§ 3. THE VACUUM CHAMBER AND SECONDARY ELECTRONS 
SUPPRESSING ELECTRODES 


The porcelain cylinders were supported ona steel chamber (Figure 1) which was 
at ground potential. Mounted on insulators across the bottom of the chamber 
there was an aluminium plate A, witha 1-inch diameter hole in it to allow the beam 
to pass to the target below. A brass sleeve B was supported from the plate by 
means of an insulating flange. ‘The inside diameter of B was 14 inches so that 
particles from the anode of the tube could not strikeit. When A and B were earthed 
high velocity particles striking A or earthed surfaces below B produced electrons, 
many of which were deflected into the accelerator tube by the weak field in the 
neighbourhood of A due to the potential on the lower accelerating electrode. If A 
were charged sufficiently positive, the field at its surface was such that electrons 
produced there were returned to it. However, electrons produced below A were 
accelerated through the central aperture when B was earthed. When B was 
charged to — 300 volts with respect to earth no electrons passed up through it. 
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§4. MEASUREMENT OF TARGET CURRENT 


The target was mounted 20 inches below plate A at the end of a 2-inch diameter _ 


tube made up of flexible bellows, brass tube and a length of glass tube. 

The ion beam reaching the target was focused into a spot 3mm. in diameter. 
In order to measure the ion current incident on the target it was necessary to ensure 
that the secondary electrons emitted returned to the target. This was done by 
supporting on the glass tube, } inch above the target, a large diameter metal ring 
which was charged to — 200 volts with respect to the target. As a check on the 
effectiveness of this method of electron suppression the target current was observed 
with negative potentials of up to 1,000 volts applied to the ring. It was found that 


— 100 volts returned approximately 85° and — 200 volts returned approximately _ 


95° of the electrons. 


§5. MEASUREMENT OF THE EXCESS CURRENT IN THE TUBE 


| 


In order to measure the excess current in the tube, it was necessary to measure _ 
not only the total current supplied by the electrostatic generator, but also the | 


currents down the bleed resistor and to the target. The difference between the 
generator current, measured by a meter in the higher voltage enclosure, and the 
resistor plus the target currents was referred to as the excess current, in other 
words the excess current was the electron current plus any unfocused ions. 
Before starting the ion source, the generator was run up to the voltage at which 
observations were to be made. The current supplied by the generator was then 
found to be equal to the resistor current, showing that no corona occurred from the 
accelerator tube and that there was no current through the tube when the ion source 
was off. 
With the electrodes A and B earthed, a series of measurements was made 
using a constant beam energy of 350 kev., while the pressure in the tube was varied 
by throttling the pump. ‘Target currents of 10a. and 20a. were used, and the 


ratio of excess current to target current was determined. The variation of this . 
ratio with pressure is plotted in curve 1 of Figure 2. It was reduced to the values — 


plotted in curve 2 by charging plate A to + 4 kv. and sleeve B to — 300 volts. 
Thus, when plate A was earthed, the electrons leaving it produced in the tube a 


current represented by the difference between the corresponding ordinates of © 
curves 1 and 2, while curve 2 represents the upper limit of the ion current lost from | 


the beam, and striking plate A. Hence, it follows that the electron current was 
greater than twenty times the ion current lost. 
With a fixed value of target current the excess current varied rapidly with the 


pressure in the envelope of the ion source. ‘This is shown by the points a and } » 


which lie off curve 1 of Figure 2, points a being for a lower pressure and b for a 
higher pressure than the value maintained for the other points on the curve. This 
will be discussed later. ; 
When the electrodes A and B were earthed, the electrons produced a fairly 
uniform fluorescence over an area about 6 inches in diameter on the plate at the top | 
ofthetube. Withthe suppressing voltages applied there was no fluorescence when | 


the pressure was less than 10-> mm. Hg but there was faint glow at a pressure of || 


5x10-?mm. Hg. Thus it was of value to know how much of the excess current | 


represented by curve 2 was unfocused ions and how much was unsuppressed | 
electrons. An x-radiation meter was therefore mounted beside the lower end of | 
the tube in a position where x-rays from the top were detected after passing through | 


} 


| 


| 
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¢ inch of steel. It was then assumed that all the excess current observed, when 
the tube pressure was 5 x 10-6 mm. Hg with the electrodes A and B earthed, was 

_ due to return electrons, and this was used to give a calibration of the radiation 
meter in terms of electron current. The radiation observed when the suppressing 
voltages were applied to A and B indicated that the electron current was negligible 
when the gas pressure was less than 10->mm. Hg, and one fifth of the target 
current when the pressure was 5 x 10-> mm. Hg; that is, about one half of the excess 
current shown by curve 2 was due to electrons. The other half is assumed 
to be due to ions scattered from the beam during transit through the tube. 
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Figure 2. 
Note.—Errors shown are due to limit of reading meters. 


§6. POSSIBLE SOURCES OF THE ELECTRONS 


As was stated above, the electron current was greater than twenty times the 
ion current lost from the focused beam. ‘This large ratio indicates that all the 
electrons could not have been produced directly by ion bombardment (see Hill ef 
al. 1939). However, they may have arisen in two ways. First, the ions incident 
on plate A released some electrons which may, as suggested by Craggs, have been 
multiplied, in their passage to the anode, by successive collisions with the electrodes. 
In fact this did not occur because the resistor current was unaltered by suppressing 
the electrons, even when the excess current and the resistor current were each 
twice the beam current. Also, the high x-ray flux produced after filtering due 
to 1 inch of steel showed that the electrons must have been accelerated through 
more than half the total potential. Secondly, neutral particles entering the tube 
from the ion source may have liberated electrons from plate A. If neutral particles 
which entered the tube with an energy due to the extractor potential of the ion 
source were responsible for the production of the electrons, the current would have 
been independent of the beam energy. If, on the other hand, they were produced 
by the lost ions, followed by some unobserved multiplication, the magnitude of 
the effect would probably not have been independent of beam energy, although 
two compensating variations may have tended to make it so. ‘The excess current 
was therefore measured with beam energies of 100kev., 350 kev. and 450 kev. 
In each case the same current ratios were observed, a fact which is consistent with 


the theory that neutral particles gave rise to the electrons. 
4 B-2 


nn 
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§7. ELECTRONS DUE TO NEUTRAL PARTICLES 

To find how many electrons were produced by neutral particles, the ions were 
deflected back to the top of the tube immediately after entering it from the ion 
source. This was achieved by charging the top two electrodes to +4kv. with 
respect to the extractor anode. Then when plate A was earthed through a meter, 
the current in the tube was equal to the electron current measured previously 
with the beam passing down the tube. These measurements were made over the — 
limited pressure range of 5 x 10-§ to 2 x 10-° mm. Hg. As would be expected, with 
plate A charged to suppress the electrons, there was no current through the tube. 
Thus only neutral particles could be responsible for the electrons. The current to | 
plate A was 4% greater than the tube current, since this fraction of the secondary 
electrons did not re-enter the accelerator tube. Again, no multiplication of 
electrons at the electrode surfaces was detected. 


§8. DISCUSSION, AND THE EFFECT IN OTHER TUBES i 

Since the cross section for ionization and charge exchange between the ions and © 
the gas varies with the energy of the ions and since the cross section is known only 
over a limited energy range, it is impossible to estimate accurately the fraction of © 
the beam which would be lost due to these processes. Taking the value of approxi- — 
mately 2 metres suggested by Craggs (1942) for the mean free path for charge - 
exchange at a pressure of 5 x 10->mm. Hg the expected loss would be of the order - 
of 30° which agrees well with the value of 20°, obtained from curve 2, Figure 2, 
after allowance is made for the electrons. Hence, it is reasonable to assume that ° 
curve 2 represents the current due to the ions formed in the tube, by charge 
exchange and ionization, plus some secondary electrons which were not suppressed. _ 

The shape of curve 1 can be explained approximately as follows. At a / 
pressure of 5 x 10-¢mm. Hg neutral particles produce an electron current equal [ 
to twice the target current. Ata pressure of 5 x 10-°mm. Hg, 0-2 of the beam is + 
neutralized by charge exchange so that in order to maintain a constant target : 
current the beam from the ion source was increased to 1-2 times its value at q 
5x10-*mm. Hg. Also the neutralized ions strike the bottom of the tube, giving + 
a further increase, which cannot be estimated, in the electron current due to neutral | 
particles. Thus, at a pressure of 5 x 10->mm. Hg, the electron current due to) 
neutral particles would be slightly more than 2-4 times the target current. The :) 
ions formed in the tube would be mainly diatomic and consequently would produce # 
about six electrons per ion (Craggs 1942) giving a further electron current of 1:2 y 
times the current to the target. Therefore the total electron current expected at a 
pressure of 5 x 10-°>mm. Hg is slightly more than 3-6 times the target current. 
The value observed from curve 1 is 4-4 after allowing for the scattered ion current. 

It is concluded that when the tube pressure is less than 10-5 mm. Hg, the excess 
current in a tube having good positive ion focusing is due to electrons produced at ' 
the base of the tube by neutral particles formed in the region of the ion source 
canal. 

The neutral particles may be formed by protons being neutralized at the walls} 
of the canal as Oliphant (1929) demonstrated in the case of helium. They are,} 
however, more likely due to charge exchange with the gas in the canal, where the 
pressure is about 10™* mm. Hg, and the mean free path is therefore about the same 
as the length of the canal. This would explain why the number increases with 
the pressure in the ion source, so displacing points a and 4 from curve 2, and why the| 
beam current shows a corresponding drop. 
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If an arc type ion source is used, the neutral particles will have a very low energy 
and may not give rise to electrons. 

A very great improvement in the performance of the accelerator can be realized 
by the suppression of the electrons. For example, the electrostatic generator used 
can only supply 160 wa. of current at 500 kv., of which 50a. pass down the resistor. 
Originally, only 25 ua. of beam could be accelerated, but when the electrons were 
suppressed it was possible to increase this to 100 za. At the same time the x-ray 
intensity in the target room was reduced from 1:0 r/8 hr. to less than 0-07 1/8 hr., 
and erratic movement of the beam was eliminated. 
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The Regenerative Deflector for Synchro-cyclotrons 
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ABSTRACT. Tuck and Teng have proposed to extract the beam from a synchro- 
cyclotron by modifying the normal magnetic field in such a way as to build up the 
amplitude of radial oscillations until protons can escape from the magnet. In the present 
paper the method is studied analytically and formulae are derived which express the 
performance in terms of the magnetic field perturbations : these formulae have been 
evaluated numerically in sufficient detail to guide practical design. 


: Se LN R OM CI Ow 
Te new method recently proposed by Tuck and ‘Teng (1950, 1951) for 


extracting the proton beam from a synchro-cyclotron promises an unusually 

high efficiency of operation. Their proposal is to modify the normal 
magnetic field at the edge of the cyclotron in such a way as to build up the amplitude 
of the usual radial oscillation of the protons about their mean orbit to a value 
sufficient to allow protons to escape from the cyclotron on the outward swings. 

The arrangement is shown in Figure 1(a). Normal synchro-cyclotron 
acceleration continues up to the radius 7, at which beam extraction is to commence. 
Outside r, the vertical magnetic field H decreases with radius in a region, of angular 
width 6,, called the ‘peeler’ and increases with radius in a region of angular width 
6, called the ‘regenerator’ : these deviations of H from the value H, of the normal 
field at r, are shown in Figure 1(6). "Tuck and Teng have shown that the required 
fields can be produced by blocks of iron fixed between the poles of the electro- 


magnet. 
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Tuck and Teng have proved the feasibility of their deflector by tracing the 
paths of particles through the system, in steps of a few degrees only, for several 
combinations of peeler and regenerator field strengths. ‘The method is in principle 


precise, though in practice the lengthy graphical computation may introduce — 


some numericalerror. The disadvantage is that the procedure is too laborious to 
permit the examination of sufficient arrangements to determine the optimum 
values of the six disposable quantities, namely r,, 6,, 9,, the angle d, and the peeler 
and regenerator field strengths, which enter into the design. 


Regenerator 
Ts Tic 
Normal Field 
Peeler 


(b) 


(a) 


Figure 1. (a) Peeler—regenerator coordinates. (b) Radial variation of vertical field H from 
radius rz, outwards. 


In connection with the design of a deflector for the Liverpool 156-in. synchro- 
cyclotron, the problem has been studied by purely analytical methods which 
readily yield a complete solution if the particle equations of motion are linearized. 
This introduces only negligible errors and the analytical solutions in good agree- 
ment with the numerical values computed by Tuck and Teng in particular cases. 
The present method allows a general survey of the problem and is particularly well 
adapted to determine the best combination of angles and field strengths. 

The main results and final formulae are in §§4, 5, 6 and 8; additional 
mathematical details are given in three Appendices. 


§2. NEED FOR LINEAR RADIAL VARIATION OF DEFLECTOR FIELDS 


The azimuth 6 and vertical displacement z of a particle are measured from the 
centre of the peeler and from the median plane of the cyclotron respectively. If 


the magnetic field is represented by a vector potential A in the tangential direction, 
the Hamiltonian for the particle is 


KH =(p, + p,”)/2m + (py/r—eA/c)?/2m 


where the angular momentum p,=mr?d0/dt+eA/c is a constant of the motion. 
Thus (p,/7 — eA/c)? is an effective potential for the motion inv and z and must bea 
minimum on the synchronous circular orbit of radius r = cp/eH. 

Let p and z denote small radial and vertical displacements of the particle from its 
synchronous orbit. It is desired to build up a radial oscillation p of large amplitude 
without at the same time increasing the vertical oscillation, which would tend to 
reduce the intensity of the extracted beam. This can be achieved if the develop- 
ment of the effective potential in powers of p and z contains no terms in pz, pz or 
ps” which couple the radial and vertical oscillations. Mathematically this requires 
that the derivatives 0°4/dr?dz and 0°A/drdz? vanish on the synchronous orbit; 


these conditions are better stated as 0*H,/dr?=0 and 0?H,/érdz=0, where. 


H,,=0A/dz is the radial component of the magnetic field. The first condition is 
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automatically satisfied since by symmetry H,, vanishes everywhere on the median 
_ plane, the second is more conveniently written in terms of the vertical field 
component Has 0?H/dr?=0. Thus (‘Tuck and Teng 1950) to minimize the coupling 
between the induced radial and the vertical oscillations the strong peeler and regenerator 
fields must vary linearly with radius as indicated in Figure 1 (b). 
When this condition is satisfied it is sufficient to use the first order Kerst—-Serber 
(1941) equations of motion for p and z; these are 
d*p r oH d*z r 0H 
Tet (1+ wae )P=9 ae - (Har )e=o eheleher sie (1) 
and may be derived from the above Hamiltonian or from the Lorentz equations of 
motion. 


§3. SOLUTION OF THE EQUATIONS OF MOTION 
‘The equations (1) must now be applied to the motion of a particle through the 
deflector, starting from the revolution at which it first enters the deflector flelds 
with synchronous radius cp/eH equal tor,. ‘Then as indicated in Figure 1 (a), p is 
to be understood as the radial displacement of a particle measured from r,, the 
increase (of the order of 0-0025 inch per turn) of the synchronous radius during the 
few (of the order of 10) cycles needed for extraction of the particle being neglected 
in comparison with the displacements p, of several inches, produced by the deflector. 
If the unperturbed magnetic field in the neighbourhood of r, varies as 7" and 
its second radial derivative is neglected, equations (1) take the more familiar form 
2 
Te A i eee (2a) 
where tr an an Bae (2b) 
are the characteristic frequencies, measured in ‘cycles per orbital revolution’, of 
the radial and vertical oscillationsingapsdandf. Similarly, for the motion through 
the peeler 


= + Q%p =0, 


cede ap —1)p,, d*2idu- +p 2—. eae (3 a) 
with Pare CEO Ori catia |e hin" bao Feke = (3 b) 


representing the logarithmic derivative of the total field in the peeler region and for 
the motion through the regenerator 


dp[d?+(q?+l)p=0y dzx/dP=qeze . ...:.. (4a) 
with oe RIT ICEL OP rere Me eta cee (4) 


representing the logarithmic derivative of the total field in the regenerator region. 
Because of the assumed form (Figure 1(b)) of the deflector fields, these 
equations (35), (4) apply only if the displacements in peeler and regenerator are 
outward. In practice the particles approach the deflector with the radial oscil- 
lation, of amplitude perhaps half an inch, induced during the previous acceleration 
and the outward maximum of this oscillation precesses in azimuth at the frequency 
1—Q=0-03 cycles/rev. Deflection of a particle begins when this precession 
brings the outward swing of its existing radial oscillation into the region occupied 
by the peeler and regenerator fields. ‘The expansion of synchronous radius 
during a complete cycle of precession is only about 30 x 0-0025 =0-075 in., so that 
every particle may be assumed to pass through the conditions stated at the beginning 


of this section. 
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The solutions of these equations are most conveniently written in terms of the 
2 x2 transfer matrices P, D, Q, F which give the radial displacement and velocity 
at exit from peeler, gap d, regenerator and gap f respectively in terms of their 


values at entry to these structures and of transfer matrices P, D, Q, F which — 


describe the vertical motion in the same way. ‘Typical examples are, for gap d, 


(\eaS0), | Soe 6) 
P/ exit P/ entry ae sin Od* cos Od* p entry 


where d* =d—14(6,+0,); and for the regenerator 


ie 
(5) =) C = emer eae t (;) Tek (6) 
<i & entry 
. hg eu qsinhg@, coshgé, < 


Then the resultant effect of one revolution through the deflector is expressed by 
the relationships (p=dp/d0, z =dz/d0) 


shen a By > AEQDPG i Sica eters (7a) 
ra es es ah pxAFODPasl hem ecm (75) 


between the coordinates of the particle at the (m+ 1)th and mth entry to the peeler. 
The effect of m revolutions, 


(Gs Quam Ge nt 


can easily be worked out by reducing A and A to their principal axes. 

Each of the above matrices D, Q etc. has determinant unity and therefore the 
same is true of the resultant matrices A and A, defined by (7). If follows that the 
product of the eigenvalues of A is unity and the product of the eigenvalues of A is 
unity. 

Let the eigenvalues and corresponding eigenvectors of A be e4, (u, v) and 


e*, (w, & so that 
A (*) =e (*) and A Ga =e oe (*:) a 3 (9) 
) v ) v 
Then since any initial radial displacement py and velocity p» can be expressed as 
@ ) =a (‘) +y (°,) ahaha (10) 
P/o @ e 


/ ! / . 
Pow — Pott q Pow — Poll 
7 7 an A ; ane 
UO =—OU Cu —UVUU 


with = 


the effect of m revolutions is derived from (8) as 


(‘) = Amy 6. +Amy (a = yp @) + yerms Ga Ieee (11) 


The deflector must be built so that one of the eigenvalues, say e4, is greater than 
one and then the other e~ is necessarily less than one. After a few revolutions 
e* is negligible and 

Pm =e up =e a eee (12) 


} 
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Thus, independently of the initial conditions of entry into the deflector, if the eigenvalues 
of A are real the particle is captured into a motion which grows in radial amplitude by 
the factor es per turn, and has the phase of the characteristic oscillation (u, v). 

Itis, of course, necessary to check that, in the characteristic oscillation (wu, v), the 
radial displacement in the peeler and regenerator is outward. 

A similar analysis may be applied to the vertical motion, with z, 3, A replacing 
p,p,A. Then if either of the eigenvalues of A is numerically greater than one the 
amplitude of the vertical motion will increase steadily; to avoid this both of the 
eigenvalues must have modulus unity since their product was shown to be one. 
Thus the vertical motion remains finite after any number of turns if and only if the 
eigenvalues of A are of modulus unity, in which case they take the form e*"’. 


§4. CHARACTER OF THE INDUCED OSCILLATIONS 

After capture in the deflector, the radial motion at the mth passage through the 

gap d has the form (Figure 1 (a)) 
p= Re sine 60) wee en ae (13) 

where y is a phase angle determined in Appendix B from the eigenvector (uw, v) of 
A and & depends on the conditions of entry. Since p must be positive at both 
peeler (# = 0) and regenerator (6 =d), the peeler regenerator angle must be less than 
7; in practice values less than 7/2 are considered. 

If the angular widths 6,,, 0, of the deflector fields are small, (say <20°), the radial 
gain e‘ and phase ¢ are given (see Appendix B) by 


ce ie F 
cosh \=cos Q(f+d)~ 7 sin Q(f+d)— xopsin Or siiOdy ee (14) 
_ eA—cos Q(f+d) —(T/Q) sin Qf cos Qd 
os cot = Sn OF +4)—(T]Q)sin OfsinQd ey 


where (2? = 1 — and, for narrow fields, S and 7 are the integrated gradients of the 
field modifications* introduced at peeler and regenerator, 

S=—O,{(7/H)OH/0r} cers T=O,{(7/H)OH [Or }regen+ +--+ ++ (16) 
The slight alterations required to deal with wide deflector fields are detailed in 
§ 8 below. 


Similarly stability of the vertical motion requires (see Appendix A, equation 
(13)) that 


cosA=cos w(d+f)+ 


2w 


should be numerically less than unity. 

Although (17) ensures that the vertical motion does not increase indefinitely, 
some vertical expansion of the beam is inevitable and this must be limited. Let 
2h, be the vertical thickness of the beam at entry into the deflector, then the 
maximum vertical displacement which occurs in the deflector for any particle is 


(see (A 10), (A11)) 


th; : 
sin w(d +f) — sasin Psi de hook (17) 


tea igelA + COS e)/(1 — cose) ye SORES (18) 
and the root mean square thickness 2A is given by (A 13) as 
JES — cry TUNE On eile eee (19) 


* In a synchro-cyclotron, the normal field gradient is very small so that the derivatives (30), 
(4) of the total field are hardly distinguishable from the derivatives (16) of the field modifications. 
The distinction would be important in application to a machine, such as a synchrotron, with a 
large normal field gradient. 
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The vertical defocusing parameter « is given by 
er sin A sin w(f+d) (20) 
SINS = T= Gos Ncos w(f Pa) Ul Sila") eta 
which is derived in Appendix A. 
Values of the parameters A, #, cos «, which fully characterize the motion, have 
been calculated as functions of the peeler regenerator strengths S and 7 and of the 
angle of separation d and may be read from the graphs (Figure 4). 


§5. REQUIRED RADIAL GAIN AND PHASE ANGLE 
To escape completely from the magnetic field, the beam must pass through a 
magnetic channel in which H is reduced. The channel may be arranged as 
indicated in Figure 2, and Figure 3 shows a cross section of the mouthof the channel. | 


Magnetic 
Cha pea wee 
ae, Regenerator 


Figure 2. Possible magnetic channel Figure 3. Cross section of channel entry, showing 
arrangement. necessary radial gain eA=OA’/OA. 


In the last turn before entry into the channel, the radial displacement p must 
increase by more than the wall thickness, otherwise all particles hit the channel 
walls. For example, in Figure 3, a particle which at the mth revolution just 
grazes the outer wall at A must enter the channel at some point A’ on the next 
revolution. ‘Then the radial gain, 

mn OA 1 AA AD 
OA OAR OAS 
must exceed the minimum value fixed by the wall thickness. 

After sufficient revolutions, every particle in the beam must pass through the 
area AWXA'YZ in Figure 3, of height equal to the mean beam thickness 2h 
given by (19), and of width WX = AA’ =ZY determined by the radial gain. The 
extraction efficiency is obviously given by the fraction of this area occupied by the 
channel entry BA’ and to achieve a good efficiency A’ must lie towards the extreme 
edge C of the channel entry. 

‘These considerations indicate the desired radial gain e4, values of about 1-4 
should be sufficient in practice and the graphs show that these can be produced by 
realizable deflector fields. 

The area WXYZ is equal to 2h.AA’ so the capture efficiency is inversely 
proportional to f, that is, it is proportional to 1/sin e. 

An upper limit to the allowable value of the vertical defocusing parameter 
cos € is given by the requirement that 2h,,,,, given by (18), must be less than the 
dee separation. The initial thickness 2h, is uncertain, but cose=4, giving 
hyyax = 1°7 hy, should be tolerable. The graphs show that adequate radial deflections 
can be achieved without exceeding this defocusing. 


a a 


ae (21) 
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The distance OA to the edge of the channel must be quite large, perhaps 3 in., 
otherwise, according to (21), an excessive radial gain would be required and there 
would also be a danger of the stray fields of the channel upsetting the normal 
orbits which should extend up to the radius 7,. If @ is the azimuth of the channel 
entry, (13) shows that particles will achieve a radial displacement 


OARin tO) cgtaetee Eien (22) 


in the last turn before entering the channel. The permissible value of this dis- 
placement is limited by the room in the tank and more seriously by the fact that it 
tends to occur in the regenerator, so that the field gradient has to extend up to this 
distance. "Thus sin (s+ Q@) must be close to unity. 

The azimuth @ of the channel entry fixes the direction of the extracted beam 
relative to the position of the peeler, and in typical buildings 6 must be small if the 
beam is to emerge in an acceptable direction without the awkward necessity of 
placing the peeler structure inside the live dee. 

These considerations indicate the desirable values of the phase %. With 
@=20°, giving reasonable clearance between peeler and channel as indicated in 
Figure 2, 4=50° giving sin (s+ Q0) =sin 70° =0-93 would be very satisfactory. 
The phase angle ¢% is sensitive to the angle d between peeler and regenerator, the 
graphs show that suitable values can be realized. 


§6. CONCLUSIONS 

The graphs of Figure 4 cover all deflector strengths and angles which are 
likely to be of interest: additional points can be calculated from the formulae of 
§4. 

Examination of the graphs drawn with 2 =0-0625 and d=30°, 60°, 90° shows 
that: (1) Given deflector strengths S, T produce greater deflections and more 
vertical defocusing at the larger peeler regenerator angles d. (11) For any 
given angle d, the greatest ratio of radial gain to vertical defocusing occurs at the 
peaks of the cose contours and this is the desirable combination of peeler and 
regenerator strengths. ‘These peaks occur for peeler strengths weaker than the 
corresponding regenerator strengths because the peeler action is assisted by the 
normal field gradient. (iii) For field strengths chosen according to (11) suitable 
phases % occur for d=60°, though a somewhat larger value would be satisfactory 
and may have to be used if it proves difficult to produce strong fields. For 
d=90°, the phase is rather too small. (iv) The ratio of radial gain to vertical 
defocusing decreases with d. 

Comparison of the graphs drawn with d=60° and n=0-0625 and n=0-09 
shows that in the latter case slightly stronger fields are needed to produce a given 
gain e4, Presumably this is because for n=(-0625 the natural frequency 
(1—n)!? of radial oscillation is closer to the frequency unity of the induced 
oscillation. Comparison of these two graphs also shows that the mathematical 
approximation made in §3, which amounts to neglecting the variation of » from 
point to point of the orbits, cannot introduce serious error. 

The deflector strengths indicated by the graphs correspond to practicable 
fields. For example, for d=60°, T=0-5, S=0-3 would be satisfactory; with | 
r,=75in., H,=18,000gauss these give 46,(0H/dr)=—72gaussrad/in. and 
6 (@H/Ar) = 120 gauss rad/in. which might be achieved by field gradients of 
—720 gauss/in. and 1,200 gauss/in., each acting for 0-1 radian. 
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§7. COMPARISON WITH RESULTS OF TUCK AND TENG 
Tuck and Teng (1950) give results of computations for two deflectors with 
fields of angular width 0-1 radian: 


(0H/Or) pecter (OH /0r) regen S L d_ Radial gain 
II — 1000 gauss/in. 750 gauss/in. 0-426 0:32 90° 1-3 
Th — 1500 gauss/in. 900 gauss/in. 0-64 0:38 60° 1-5 


the radial gain being estimated from their figure 8.15. For these values, and. 
n=0-0625, the graphs give e4 = 1-3 and 1-4 respectively, in reasonable agreement. 


Ee Voie) ‘ e45156 37° " A 
/ 
th + 
746°/7)5157/ 0-137 31° $ 
ie Zot] eM stay 3 2G eA cpaas wks 


0-4) 
0:2 / 
i FS 
aa Ss Se Y 1 1 
“06 -04 -02 0 0:2 0-4 06 0:8 1-0 12 


ee ae a eS 


Figure 4. The radial gain e4, phase % (degrees) and vertical defocusing cos € in terms of the peeler 
and regenerator strengths S and T. 
(a), (b), (c) : n=0:0625, d=30°, 60°, 90° respectively. (d): n=0-09, d=60°. 


Tuck and Teng’s points do not satisfy our condition (ii) for the minimum 
vertical defocusing. This is attributed to the fact that the defocusing of any parti- 
cular particle depends on the phase of its oscillation at entry to the peeler and that 
the present calculation considers an average over all phases whereas that of Tuck 
and ‘Teng does not. 

The present criterion (ii) indicates peeler strengths weaker than those required 
by Tuck and Teng. As the peeler field is difficult to make, this point may be 
important for construction of a working deflector. | 


38. ALLOWANCE FOR WIDE PEELER-REGENERATOR FIELDS 
Since only limited field gradients can be produced experimentally, it may turn 
out that the required deflector strengths S, 7'can most conveniently be achieved by | 
using wide peeler and regenerator fields. In this case the simple definitions (16) | 
must be replaced by (see Appendix C) 
sinh 6,,(p? — 1)? 
Ca 


sinid (gq? + bye? 


S=(p?—n) (erie? cere ee 


NES past) 
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for the radial motion and by 
sin pé,, : sinh g@ 
Pe aing 2 TD eee oie: (24) 
for the vertical motion. 

Now —p? and q? were defined by (3d) and (4d) as logarithmic derivatives of the 
total field, — is the logarithmic derivative of the normal field. Therefore 

(p?—n) = —{(7/H)OH/Or}yecrer and (P+n)={(r/H)OH/Or}regen ++ +++ (25) 
are precisely the logarithmic derivatives of the field modifications which previously 
appeared in (16). Equations (23) and (24) differ from (16) only by the factors 
p' sin p@, etc., which modify S, T by +4S?0, and +47?6, approximately ; these 
differences are negligible for S, T>4 and 0,, 6,<}. 

With the definitions (23), equations (14) and (15) still describe the radial 
motion if f and d are replaced by f’ =f+(T0,?— S0,,)/24, d’ =d+(T0,? — S0,,7)/24 
to allow for the effective thickness of peeler and regenerator. Similarly, equations 
(17) and (20)still describe the vertical motion if fand dare reduced by(T6,?— S6,,?)/24. 
These corrections to f and d are always negligible, so the graphs can be used for 
wide fields provided that S, T are interpreted according to (23) and (24). 


SoSAPPLICATION LO SYNCHROTRONS 

The application of the regenerative deflector to synchrotrons has not been 
examined in detail. The analysis developed in this paper is applicable and suggests 
that in its present form the deflector will not work with the large value of found in 
synchrotrons, but that a successful modification might be possible. 

Consider n = 0-6 giving Q =(1—n)1? =0-63 and w =n? =0-78 for the character- 
istic frequencies of radial and vertical oscillation : thus 27Q, = 227° and 27w =280°. 
Then (14) gives 

cosh A = — 0-68 + 0:57(T— S) —(TS/2Q?) sin QfsinQd. ...... (26) 
Since cosh A greater than unity is required for radial deflection it is best to choose 
S =0, that is to rely on the large normal field gradient for the peeler effect. Then 
(26) requires T greater than 2:9 for radial deflection to occur. However (17) 
gives cosA=0-18—0-637T—(TS/2w?)sinwfsinwd —...... (27) 
so that the vertical motion would be unstable with such a strong regenerator. 
Similar conclusions hold form =0-5 andn=0-7. ‘The reason is that the free vertical 
oscillations of frequency w more nearly resonate with the frequency induced by the 
deflector than do the free radial oscillations : this situation is just the reverse of that 
in the synchro-cyclotron where is small. 

It may be possible with large values of nm to design a regenerative deflector to 
build up the vertical oscillations of the particles and to deliver them to a magnetic 
channel placed above the normal position of the circulating beam. At the same 
time the radial oscillations must not increase. ‘This change in the role of the 
radial and vertical motions requires only obvious modifications of the analysis: 
the most important are that in (14) cosh A must be replaced by cos A which is 
required to be numerically less than unity and that in (17) cos A must be replaced 
by cosh A, where e’ is the ‘ vertical gain per turn’, 
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APP EN Dieses 
THE VERTICAL MOTION 


For narrow deflector fields the peeler and regenerator matrices, introduced in 
§ 3, take the simple form 


‘al few 
P= i me o-(5 oF cae (Al) 
with S, T given by (16). Then 


A=(* - bier (A2) 


(i Sit 
with « =cos w(d+f) — 5 sineo(d +f)+ * sin wf cos wd— —z Sinwf sin wd 
i: Lay 
f= ~, Sinw(f +d) + —3 Sin wf sin wd 
y=—wsinw(f+d)—Scosa(f+d)+ T cos wf cos wd — ~* cos wfsinad 


5=cosw(f+d)+ Eerie 
(63) 


The eigenvalues e+" of A are determined by 
ot S644 ott 2 COs ai ad, ea) aie (A3) 


which gives equation (17). ‘The two eigenvectors are given by 


A (, ) =exi? (oa) bait leek ocareea ee (A4) 


Then, if the initial vertical displacement and velocity of a particular particle are 
expressed as 


(),-e( Si) ter() tera( a). 
where the phase angle ¢ is fixed by 
(3/%)o =(cosb—sin btan ¢)/apocy Wee ees (A6) 
the motion after m turns is 
(: ) = & ie pA Mo ye ae oS (A7) 
The amplitude of vertical oscillation after 7 turns through the deflector is defined by 
h, 2 22-7 (2, [ey 2 On ee Te eens (A8) 


and may be calculated from (A7) as 


C2 
h,2= Se (1 + aw?) {1 + cose cos 2(A+u+mnr)} 


2awsinb sin 2b 
where sin «= ——,—_ and_ tan 20 = ———_—___ 
14+ aw? Ye wae + cos 2b 


Thus h,, is related to the amplitude hy at entry by 


h,* _1+cosecos2(P+ndA+p) 
gt Theos C068 DEE (A10) 
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Ifa family of particles all enter the deflector with given amplitude fp, the greatest 
_ vertical displacement of any member is h,,,,, given by 


hnax |Mo = (1 + cos )/(1—cos e) Serer ane (A11) 


which is equation (18) above. 
For any given particle, the mean square amplitude of oscillation averaged over 
many turns is 


ho2/{1tcosecos2(Pt+p)} eee (A12) 


which for different particles varies between h,2/(1+cose) and h)?/(1—cose), 
depending upon the phase of vertical oscillation at entry to the deflector. For 
some particles the mean amplitude of oscillation is increased by the deflector, for 
others it is reduced; this compensation is a consequence of Liouville’s theorem. 
The geometric mean of these two limits is 


= he il coset = hse, eet ot oe (A13) 


which was used in (19) above to represent the mean square amplitude of the whole 
beam. 
It remains to determine the characteristic angle « in terms of the deflector 
strengths. Equation (A4) gives 
asinX —sina 
B ya 


and then (A9) yields sin e =2 sin A/(Bw —y/w) which leads to (20) if the values (A2) 
of 8 and y are substituted. 


sinb= 


Ge PereNeDITEX 3B 
THE RADIAL MOTION 


For narrow deflector fields, the peeler and regenerator matrices are 


P=(% igh Q=(_5 iar Need (B1) 


with S, T given by (16). Then 
5S & _) PY (B2) 


_ is given by equations (A2) with w replaced by and S, T replaced by — S, —T. 
Equation (14) for the radial gain follows immediately and in (9) one may set u=B, 
v=er—a, 

From (12) it follows that, at the mth exit from the peeler, the radial displacement 


and velocity are 


e m 10 os md as — hems sin ip 
(5) = Ke As ne) =e ‘ee =ke ee) ees (B3) 


say, with cot =(Su+v)/Qu=(SB + e4 — a)/QB which is given explicitly by (15). 
Equation (B3) leads immediately to (13). 
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APP Ec NtDEb x6 


REDUCTION OF THE DEFLECTOR FIELDS 
The matrices P, Q, P, Q introduced in § 3 to describe the deflection produced by 
fields of any width may be simplified by an analysis similar to the reduction of a 
thick lens to its principal planes. 


Let cos Od o>? pint m) 


cals (2 QsinQd cos QP 
represent the transfer matrix for the radial motion through azimuth¢. Thus, in 

Then any of the deflector matrices such as (6) may be reduced to the simple 
form of (A1) or (B1), by means of the general formula 


(7 2) =na9 (| 9) MGB. eee (C2) 
with the effective strength x given by ¥=c+bQ*" ~~ ~ | | “alae (C3) 
and a thickness ¢ given by 

O-* sin 06 +x0-* sin? 4O¢=p-—— ee (C4) 
Thus P-IGs,) PGs). C= Ns. ON (C5) 
with P= (‘5 ar Q=(_F ‘a my (C6) 


where S and T are given by (23). Because of the group property 
II(d*)U(3¢) = (d* + 34), 
there results A=NPOUNGIR” ~ i eee (C7) 
where A’ = II($¢,,) AII(—3¢,) relates the radial coordinates of the particle at its 
(m+ 1)th and mth passage through the centre of the peeler, and 
ee =f+ 2($py ae am Oe oe 84); d'=d+ (> + $4 - 0, < 0). peri (C8) 

f and d’ are given explicitly by equation (25). | 

Actual calculations are much simplified by working with the matrices P’, 
Q’, A’ and all the formulae have been derived in this way. 


A similar reduction applies to the vertical motion, Q being replaced by ow. 
In place of (C6) there results 


ea pase a pele ets 
P-(_5 ae o'= (7 ) ee (C9) 
with S and T given by (24). 
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LETTERS TO THE EDITOR 


The Luders Deformation of Mild Steel 


The writer has presented evidence (Hall 1951) which seemed to show that the Liiders 
bands appearing during the deformation of mild steel can each be characterized by a uniform 
shear front which propagates uniformly along the specimen. ‘The shear displacement at 
the edge of the band is straightened by the grip constraints, and in addition to these two 
modes of deformation, shear and bending, a certain amount of creep extension also appears 
to be present. 

In a tensile test at room temperature, the shear displacement will be straightened. 
plastically, for the yield stress of the material which has just sheared over at the edge of the 
Liiders band is certainly no greater than the applied (lower yield) stress. Thus, a specimen 
removed from the testing machine before the end of the Liiders extension appears straight, 
apart from a small angular kink of the order of 2° at the band edge. 

The Figure (a), shows a plastically deformed thin strip specimen of mild steel laid 
beside a straight edge. The dimensions of the strip are 135 x 2x $ mm., and the wide front 
face is shown in the photograph. No evidence of bending can be seen. 


il mn i | 


0) | eae 
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Strip specimens of mild steel deformed (a) at room temperature; (5) at 200° c. 
Strain rate: 0-045/min. 


Recently, however, these experiments on the deformation of mild steel have been 
extended to temperatures in the ‘ blue-brittle’ range, around 200° c. Here the deformed 
area of the specimen just behind the Liiders front ages very rapidly during the test, so that 
the yield stress in this region is much above the applied stress. ‘Thus, the bend induced in 
the specimen is only straightened elastically, and on removing the specimen from the machine, 
a marked curvature remains in the test piece. 

(b) shows a specimen with the same mechanical properties as (a), deformed at the same 
strain rate, but at 200°c. instead of at room temperature. A distinct curvature is seen in the 
plane of the wide front face of the specimen. ‘Thus, mild steel may be bent by straining it 
in tension. 

The radius of curvature of the bent specimen is very nearly 1 metre, so that there is a 
macroscopic strain inhomogeneity of about 0-2°/, between the two sides of the specimen. 
The inability of large specimens to accommodate such bending moments must undoubtedly 
play a part in determining the form of the Liiders band in the test pieces, and is perhaps. 
the reason why in specimens of large cross section under a high load, the Liiders bands 
present a very complex appearance, while in these thin strips or wire specimens, usually only 
one or at most two Liiders fronts are observed (Hall 1951). 
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This work forms part of a project on the deformation of mild steel in the blue-brittle 
range, which was undertaken for the British Iron and Steel Research Association. The 
writer is indebted to the Association for permission to publish this letter. 
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of these experiments will be published elsewhere. 


Cavendish Laboratory, 
Cambridge. 
10th October 1951. 


Hatt, E. O., 1951, Proc. Phys. Soc. B, 64, 742. 


Ultrasonic Absorption in Supercooled Liquids 


In December 1949 my attention was drawn to work on the existence of a discontinuity 
in the temperature coefficient of viscosity of di-phenyl ether in the critical region near the | 
melting point (Dodd and Hu Pak Mi 1949). 
changes in the liquid, a marked discontinuity may be expected in a suitable function of the 


If this phenomenon be due to structural 


ultrasonic absorption in a liquid on its being supercooled. 


Figure 1 shows preliminary results, obtained in this laboratory, for the ultrasonic 
absorption and wavelength in menthol at a frequency of 68 Mc/s. 
absorption only in salol and di-phenyl] ether. 


A full account 


E. O. Hare 


Figure 2 shows 


I have termed ‘classical’ that part of the 


absorption which is due to the Stokes’ shear viscosity term. This has been calculated 
from the measured wavelength and viscosity, the variation of density being ignored. 
While the measurements reveal no obvious discontinuity when absorption is plotted 
against temperature, it will be seen that all three liquids show a marked increase in 


absorption with approach to the freezing-point. 


It appears likely that the measured 


absorption in menthol passes through a minimum at a temperature rather higher than 


could be attained with the existing apparatus. 


The small temperature coefficient of 


absorption in di-phenyl ether makes the results for this liquid the least reliable: repeated 
experiment with different samples appears to confirm the existence of a minimum, but the 
temperature at which it occurs is ill-defined. The curves of wavelength against temperature 
give more evidence of discontinuity than in the case of absorption and lend themselves 
more readily to accurate measurement. 
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‘It was thought that the observed increase in absorption on approach to the freezing 
point could be due to scattering at crystal nuclei. 
menthol at the higher frequency of 87:5 Mc/s. to test this hypothesis. At all temperatures 
the ratio of the measured absorption at the two frequencies used was close to 1-65, in 


accordance with a square law. 


An experiment was performed in 


If appreciable scattering were present an approach to 


2°72 might be expected, in conformity with a fourth-power law. 
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The three liquids were chosen for convenience of melting point, ease of supercooling 
and availability of viscosity measurements made in them by Dr. Dodd and his co-workers. 
‘The frequency range used gives an ultrasonic attenuation of the order of 20 db/mm., so 
requiring only a few cm® of sample liquid. The pulse technique is used (Pellam and 
Galt 1946, Pinkerton 1949), with separate transmitter and receiver transducers and a 
six-inch fused quartz acoustic delay line (Rapuano 1947). Velocities are measured by 
injection of a phase-locked signal from a generator and observing the regular changes in 
output as the receiver is moved. 

I am much indebted to Dr. Pinkerton for valuable advice with the design of the 
apparatus and to Dr. Dodd for viscosity data supplied prior to publication. It is hoped 
shortly to publish detailed results taken with a new apparatus nearing completion. 


Physics Department, A. N. HUNTER. 
University College of Leicester. 
26th September 1951. 


Dopp, C., and Hu Pak M1, 1949, Proc. Phys. Soc. B, 62, 454. 
Lirovitz, T. A., 1951, ¥. Acoust. Soc. Amer., 23, 75. 

PELLAM, J. R., and Gatt, J. K., 1946, #. Chem. Phys., 14, 608. 
PINKERTON, J. M., 1949, Proc. Phys. Soc. B, 62, 286. 
Rapuano, R. A., 1947, Phys. Rev., 72, 78. 


The Optical Constants of Lead Sulphide and Lead Telluride 
in the Region 0°5—3 microns 


In a recent communication (Paul, Jones and Jones 1951) measurements of the absorption 
characteristics of PbS in the wavelength range 2-12 microns have been described, and in 
addition the reflection coefficients of the material in this region have been measured. In 
this laboratory some experiments are in progress to measure the refractive and absorption 
indices of lead sulphide and related materials in the visible and infra-red spectral regions. 

The method used depends on the measurement of certain reflection characteristics of 
‘specimens of the materials. The ratio of the reflection coefficients for incident light 
polarized in and perpendicular to the plane of incidence is measured for various angles 
of incidence, and, with the aid of graphs computed beforehand, the optical constants can 
‘be derived from these measurements. The method is an improvement on an earlier one 
-due to Collins and Bock (1943), and it is hoped that complete details will be published 
shortly. 

Measurements have so far been restricted to wavelengths below 2:6 microns, and have 
‘been made on polished surfaces of polycrystalline specimens of PbS, and on one polished 
cleavage face of an artificially prepared single crystal of PbTe. 

The refractive indices and absorption constants for two specimens of PbS are shown 
in Figure 1. One of these specimens was pure Sardinian galena, and the other a piece of 
galena of unknown origin. ‘The measurements in the visible region of the spectrum were 
made using a classical polarimetric method. ‘The error in the refractive index values is 
some 4%, or less, except at the extreme long wavelength end of the range, and is 10% in 
the absorption constant values. 

Comparison of these absorption constant values with those of Gibson (1950) for 
chemically deposited layers of PbS, shows that for the present specimens the absorption 
edge occurs at a longer wavelength than for the layers, and that the long wavelength ‘ tail’ 
to the edge, observed by Gibson, is very pronounced in this case. ‘The magnitude of the 
absorption constant at wavelengths below 2 microns is considerably greater than for the 
chemically prepared layers. 

Apart from a slight increase at about 0:8 micron, the refractive index shows no marked 
variation, rising gradually to values of about 4 at 2:5 microns, in agreement with the values 
quoted by Paul, Jones and Jones. 

The surfaces of the specimens were first ground flat using successively finer grades of 
abrasive, and finally polished using a cloth lap and metal polish. ‘The various specimens 
‘were ground and polished in as nearly as possible an identical manner, and the values 
obtained from four surfaces prepared separately agree to within 10-15%. 

AC=2 
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Figure 1. PbS, polished crystals. Figure 2. PbTe, polished single crystal. 


The values of the refractive index and absorption constant for the polished single 
crystal of PbTe are shown in Figure 2. The refractive index rises from a value of 3-8 at _ 
i micron to 5-2 at 1:8 microns, remaining then roughly constant at about this value. The 
values of the absorption constant obtained are slightly higher than those obtained by 
Gibson for layers of the material prepared under various conditions. 


Physics Department, 
Royal Holloway College, D. G. Avery*. 
Egham, Surrey. 
27th August 1951. 


* Now at Telecommunications Research Establishment, Gt. Malvern, Worcs. 


Co uns, J. R., and Bock, R. O., 1943, Rev. Sci. Instrum., 14, 135. 
Gipson, A. F., 1950, Proc. Phys. Soc. B, 63, 756. 
Pau, W., Jones, D. A., and Jongs, R. V., 1951, Proc. Phys. Soc. B, 64, 528. 


The Refractive Index in Electron Optics 


In 1933 Glaser derived a correct and now well-known value for the refractive index 
in electron optics. Although quite generally applicable, his is a particular value in so far - 
as it omits an arbitrary function, and by a rational and indeed inevitable generalization _ 
suggested by Glaser’s derivation, Opatowski (1943) obtained other values. Some of these, 
however, are wrong; and in 1949 we showed which of his generalizations are valid and how __ 
Glaser’s approach must be modified in order to be theoretically sound and to lead always 
to correct results. As Glaser (1951) suggests that our discussion was based on 
misconceptions we should like to make a brief statement. 

(a) As introduced into mechanics the Lagrangian L and the energy E are undefined to. 
within, in general, different arbitrary constants. Opatowski saw, therefore, that the value 
p.—(L--£)/v as introduced by Glaser in 1933 has no better footing than the more general 
expression 4=(L--E-+K)/v where K is an arbitrary constant. The relevant and in fact 
very simple point which Glaser does not seem to appreciate is that either expression is. 
correct only if E and L are so normalized that the numerator equals 2q;dL/@q;. This. 
normalization is not automatically guaranteed and it is not explicitiy stated by Glaser. 

(b) ‘That the time-dependent principle of Hamilton could not directly give a refractive 
index was stated clearly by de Broglie in his foundation papers on wave mechanics; he saw 
that Maupertuis’ principle must be used. Glaser (1951) appears now to recognize this. 
point but he makes no reference to the fact that in his original paper (Glaser 1933), which 
he still upholds, he refers to Hamilton’s principle only. There he eliminates from L the 
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time-dependent quantities by means of the energy integral so that Ldt=.ds—Edt; pisa 
function of EF, of the coordinates, via the electric and the magnetic potential, and of s, the 
unit vector in the direction of the trajectory. Then he continues: ‘‘ Hence (5) (i.e. 8 fLdt=0) 
becomes equivalent to d5fuds=0 as (5) is not varied at the limits and thus 5{Edt=0.”’ 
According to this statement [uds must also be varied so that at the limits the time as well 
as the coordinates are kept constant. But the whole point of Fermat’s and Maupertuis’ 
principles is that the variation is carried out regardless of the time at the limits. Hence, 
the » derived by Glaser is not a refractive index at all, and his claim of having derived it 
correctly is unfounded. 

(c) Contrary to the statement of Glaser we made no suggestion either in connection 
with Opatowski’s paper or otherwise that the general refractive index could not be 
specialized to a case of axial symmetry. Since this special (‘isotropic’) refractive index 
can be derived from the general refractive index, Glaser is not justified in calling the 
former a ‘ generalization’ of the latter and his deduction from first principles can add as 
little to the elucidation of the situation as can the particular derivation of Maupertuis’ 
principle he proposes. 


Birkbeck College, W. EHRENBERG. 
University of London. R. E. Sipay. 
14th September 1951. 


EHRENBERG, W., and Sipay, R. E., 1949, Proc. Phys. Soc. B, 62, 8. 
GLasER, W., 1933, Z. Phys., 80, 450; 1951, Proc. Phys. Soc. B, 64, 114. 
Opatowskl, I., 1943, Phys. Rev., 65, 54. 


The fact that the Lagrangian is definite with the exception of a total differential and 
the transition from Hamilton’s principle to the principle of Maupertuis belong to the 
well-known fundamental theorems of analytical mechanics. We refer in this respect to 
any textbook. Therefore on this point a difference of opinion would be hardly possible. 
That our derivation (Glaser 1933) of refractive index in electron optics is based on these 
well-known principles, was fully shown in the paper quoted (Glaser 1951). However, 
we are not in the position to assume responsibility for the correctness of papers, published 
later by others, which claim to be based upon our work. 


Technische Hochschule, W. GLASER. 
Vienna. 
25th September 1951. 


G.asER, W., 1933, Z. Phys., 80, 450; 1951, Proc. Phys. Soc. B, 64, 114. 
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Obituary Notices 


LESLIE JOHN COMRIE 


THE death of Dr. L. J. Comrie, F.R.S., in London on 11th December 1950 at the age of 
fifty-seven has removed an original and independent figure from the ranks of British 
scientists. 

Leslie John Comrie was born on 15th August 1893 at Pukekohe, Auckland, New 
Zealand. He was educated at Auckland Grammar School and Auckland University College, 
where, after a promising career, he graduated M.A. in 1916, with First Class Honours in 
Chemistry. He served in the First World War with the New Zealand Expeditionary Force, 
incurring severe wounds, with consequent amputation of a leg. It was doubtless this 
incapacitation, together with the handicap of increasing deafness, that caused him to recon- 
sider his career, for after the war he remained in Britain as a New Zealand Expeditionary 
Force Scholar,.studying astronomy and computation at St. John’s College, Cambridge, where 
he was an Isaac Newton Student from 1920 to 1923. In 1923 a thesis, devoted to the history 
and detailed examination of methods of predicting the occultation of stars by planets, 
gained him the degree of Ph.D. 

Two years in America followed, the first as Assistant Professor of Mathematics and 
Astronomy at Swarthmore College, Philadelphia, the second as Assistant Professor of 
Astronomy at Northwestern University, Evanston, Illinois. In 1925 he was appointed 
Deputy Superintendent of H.M. Nautical Almanac Office, Greenwich, and in 1926 entered 
upon his duties there. He was no stranger to Greenwich, having already, in 1920, spent 
two months at the Royal Observatory as part of the conditions of his Isaac Newton student- 
ship. In 1930 he became Superintendent, a post which he held until 1936. Without 
delay he used the opportunity to embody in the next issue of the Nautical Almanac, for 1931, 
certain changes which he had without doubt considered for some years, notably the use of a 
standard equinox. The Almanac of 1931, breaking away from the tradition of almost a 
hundred years, incorporated, among other radical changes, the new ephemerides. 

Comrie’s imagination, essentially direct and pragmatic, had been captured, from his 
first contact with calculating machines, by their unexplored possibilities. He saw that with 
the rapid improvement in calculating machines many accepted procedures, legacies from the 
pre-machine era, were wasteful in time and effort and should be superseded or re-drafted in 
the newer idiom. He reconsidered such processes as interpolation, numerical integration and 
subtabulation from first principles, and devised simple and powerful repetitive techniques, 
requiring the minimum of copying and capable of being mastered and understood by 
operators of only moderate mathematical training. His ingenious mind ever busied itself 
with inventing methods for turning the latest devices, such as transfer levers, split keyboards 
or twin machines, to practical use. It was characteristic of him, for example, to harness 
together National Accounting Machines in order to build up and print tabular values from 
differences of high order, or to use punched cards and the Hollerith machine, as he did in 
1928 with E. W. Brown’s Tables of the Moon, to sum large numbers of harmonic terms. His 
methods, severely simple and simply described, are explicit in the prefaces or introductions 
to the various tables published under his editorship, notably in Interpolation and Allied 
Tables, reprinted from the Nautical Almanac for 1937. 

In 1936 he resigned from the Superintendentship and began a new and independent 
venture, the founding on business lines of a computational firm, the now celebrated Scientific 
Computing Service Limited. 'This, at the time, was a bold step, though soon justified by 
increasing success, by requests for computational help on a large scale by the universities, 
by industry and by Government departments.. With the outbreak of the Second World 
War the demand, official and unofficial, for such services was extraordinarily accelerated. 
The office in Bedford Square served as training ground for many computers, who entered 
with moderate qualifications and passed out as experts. By 1944 the staff numbered thirty; 
and after the war, with the business now firmly established and enlarged, Comrie found the 
need of an assistant manager, of similar interest and authority in all matters concerning 
computation. ‘This was Dr. J. C. P. Miller, who now succeeds him. 

Comrie had now for many years been recognized as an authority of international 
reputation on mathematical tables. He had been the first Director of the Computing 
Section of the British Astronomical Association, Secretary to the Mathematical T'ables 


\ 


: 


Obituary Notices I0gI 


Committee of the British Association, 1930-36, member of the Committee on Mathematical 
‘Tables of the National Research Council, Washington; and from 1948 to 1950 he was a 
member of the Mathematical Tables Committee of the Royal Society. He was elected a 
Fellow of the Royal Astronomical Society as early as 1919, and served on its Council 
from 1929 to 1932. He was also a Fellow of the Royal Statistical Society, and of the Physical 
Society. He was an active member of the International Astronomical Union, attending all its. 
conferences from 1928 onwards, and was elected President, in 1932, of its Commission on 
Ephemerides. His native country recognized him by honorary fellowships, of the Astro- 
nomical Society of New Zealand in 1933 and of the Royal Society of New Zealand in 1945. 
In March 1950 he was elected a Fellow of the Royal Society, but lived only ten months to 
enjoy this culminating honour. 

His papers, though they number about a hundred, represent only a small part of his’ 
work. A proportion of these, mostly appearing in the Monthly Notices of the Royal 
Astronomical Society, are on astronomical matters, but the larger proportion (often appearing” 
in such journals as the Office Machinery Users’ Association Transactions) deal with table- 
making and the application of various types of machine to problems of computation. Of- 
great importance are his reviews and notes on errata in existing mathematical tables, a 
characteristic example being the article in Mathematical Tables and Aids to Computation, 1 
(1945), pp. 391-399, on Jahnke and Emde’s Tables of Functions. His erudition in this field 
was immense. 

The tables for which he was personally responsible bear on every page the marks of his 
exacting standards, imposed not least upon himself. First, by rigorous systematic 
checking and counterchecking, the values must be impeccably accurate; they must also be 
interpolable with the maximum of convenience, readable with the maximum of ease, and 
typographically elegant. No one can doubt, perusing, for example, the third and fourth 
enlarged editions of Barlow’s Tables, issued in 1930 and 1941, or, perhaps most con- 
spicuously, the new Six-Figure Mathematical Tables in two volumes, logarithmic and natural, 
published by Chambers in 1948, that such aims have been achieved to the full. These last 
tables call for special mention; they resulted from the decision of the publishing firm, on the 
centenary of the old tables in 1944, to take the advice of a modern expert in regard to their 
re-issue. They are rightly described in the Editor’s preface as “‘ one of the greatest of all 
British table-making enterprises ’’. 

In character Comrie was intolerant of imprecision, direct and outspoken to a degree: 
sometimes startling, though this was really a concomitant of his essential simplicity and 
singlemindedness. ‘The writer of the present notice was never a fellow committee-man with 
him; but it is no secret that his abruptness of manner in delivering an opinion, his refusal to: 
use any of the ordinary arts of compromise, did not always meet with a favourable response. 
This was, however, only one side of the man, reflecting his inveterate honesty. 

Perhaps one who was bound to Comrie by ties of friendship, of birth in New Zealand at" 
about the same period, and of similar service and similar though less severe experience in 
the First World War, may conclude on a specially personal note. Comrie as a young man,. 
before departing for France and suffering so severely, was physically vigorous and energetic, 
a good rifle shot and tennis player, pursuits in which his qualities of precision and decisive- 
ness found a ready outlet. It is more than probable that the disabilities that permanently 
precluded such active outlets produced a certain frustration, which in its turn was overcome 
by the unremitting energy which Comrie poured into everything that concerned astronomy, 
machines, tables and computation; for, as has already been stated, the amount of work he 
did in his lifetime cannot by any means be summed up in terms of publications. Not least 
was the help he unobtrusively and willingly gave to others in the field, of which a fitting 
memorial exists in the Index of Mathematical Tables of Fletcher, Miller and Rosenhead, 1945, 
already classical. New Zealanders visiting London, and astronomers and computers from: 
America and Europe can testify ‘o the warm hospitality he gave them, ably seconded by his 
wife ; but very few knew of his many acts of secret and private generosity, some of which have 
only come to light since his death. The picture that lingers most in the mind is of Comrie as 
host, cut off by deafness, it is true, from the cross-currents of conversation, but himself 
sitting quietly happy at the evident enjoyment of his guests. It is a question whether, in 
giving full rein to his proclivities, he did not drive himself too hard along a single track; but 
his work is of great, permanent and visible value, and he has set standards, in all that concerns: 


table-making and computation, that can hardly be surpassed. 
A. C. AITKEN. 
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LEONARD GEORGE GRIMMETT 


I worKED as-assistant to L. G. Grimmett for 18 months, 1934-35, in the Physics 
Department of the Radium Beam Therapy Research at The Radium Institute, London. 
I have had as transient a contact with many others, yet they have passed out of my 
recollection and interest. In contrast, Grimmett made a lasting impression. I have always 
felt a surge of pleasure on meeting him again, or on hearing of his exploits, and was very 
grieved to learn of his sudden death from heart failure on 27th May 1951 in Houston, 
‘Texas. He had a unique and admirable personality. 

The two pillars of Grimmett’s career were his love of science and of music; he was a 
pianist of the first rank. On leaving the Holloway County Secondary School he paid the 
expenses of a physics degree course at King’s College from his earnings in restaurant and 
theatre bands. He made a voyage as ship’s pianist to Buenos Aires. Undoubtedly this 
double life was a heavy strain. To me he deplored that, at a time when he was engaged 
in this struggle, his Professor was completely discouraging, believing that a degree could 
not be taken in this eccentric fashion. Nevertheless it was taken (in 1926), and so well that 
he was encouraged to do three years’ research. In 1930 he became Assistant Physicist 
at the Westminster Hospital Radium Annexe, working under Professor H. T. Flint, and 
together they developed one of the earliest hospital physics departments. Grimmett was 
a pioneer in the development of radium beam units for the treatment of malignant disease; a 
two-gramme unit designed by him in the early 1930’s is still in use. Methods of auxiliary 
physical measurements, such as those of dose distribution round these units, which are 
now routine, also claimed his attention. 

He collaborated with Dr. F. G. Spear (1932) in an early quantitative study of the 
influence of gamma-ray intensity on the inhibition of cell division in tissue cultures (Spear 
and Grimmett 1933). This work caught the eye of Professor McLennan, Chairman of the 
Executive Committee of the Radium Beam Therapy Research, who persuaded him to 
join that body. He served it, and the Radiotherapeutic Research Unit of the Medical 
Research Council which assimilated it, until 1944. 

Grimmett excelled in the design of apparatus; he was an artist in temperament. The 
common view that if a scientific apparatus works it is satisfactory was foreign to him. 
It must not only work but be well designed, elegant, and properly finished. Furthermore, 
each stage in its development must reach the same standard. It was not sufficient that 
the drawing should be adequate for the mechanic’s use; the penmanship must be excellent 
also. He saw everything in a big way. The idea of improving the collimation of a radium 
beam by fitting a solid gold nozzle to the teleradium unit was typical, and he made it a 
reality (Flint, Grimmett, Rock Carling and Cade 1934). He spent time and money 
without inhibition or fear to get the results he wanted. He supplied most of the ingenious 
ideas for the 5-gramme radium unit (later 10 grammes), with pneumatic transference of 
the radium (Grimmett 1937), replicas of which have been installed in hospitals all over the 
world. He was one of the first to realize the importance of the total absorption of radiation 
by a patient (now known as ‘ volume dose’), as distinct from the absorption in the 
malignant growth itself. He tackled its measurement practically by the construction of a 
full-size human shape of parallel celluloid sheets with narrow air gaps which formed 
ionization chambers (Grimmett 1942). 

He believed in the need for higher voltages in x-ray therapy, and in 1939 began the 
design of what ultimately became the first high pressure Van de Graaff generator to operate 
in Britain. Harassed by the war, the work was moved from the bombed Radium Institute 
to the Imperial College of Science and then to the Hammersmith Hospital, advancing 
slowly, since each stage must reach perfection without short cuts. Concurrently he 
supervised the physical side of a long clinical experiment which compared the relative 
merits of radium beam and 200 kv. x-ray treatments, with the fields of both types arranged 
to have closely matched dose distributions. At Hammersmith he built up a fine team 
and a magnificently equipped physical laboratory. 

Then in 1944 he disappeared from his usual scenes, reappearing briefly from time to 
time from Paris, or wherever else his service for the British Council, UNESCO, and UNO 
had taken him. Finally it was reported that he had settled in Houston, Texas, where he 


was developing a radiobiological research laboratory in the really big way which suited 
his temperament. 
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Grimmett’s career might suggest the dominant, aggressive personality, yet he was 
kind, gentle, always soft-spoken, and quite imperturbable except in sympathy for others. 
As we walked down Regent Street in 1934 we passed seedy, out-at-elbow violinists playing 
in the gutter. He was filled with distress—‘“ I know this kind of man; I have played with 
them. They are good musicians displaced by the ‘Talkies’. He declined to assist in 
the atom bomb development. ‘‘ The M.R.C. were approached”’, he said, ‘‘ but I asked them 
to squash the idea’’; then, with a wry grin in oblique explanation, “‘I don’t mind killing 
‘Germans in odd ones’”’. 

Grimmett showed that, with few initial advantages other than talent, courage, and 
tenacity, he could be a ‘ 20th century’ adventurer in the best sense, and at the same time 
make substantial contributions to human welfare. All who knew him will miss him 
greatly, and feel deep sympathy with his widow. 

JOHN READ. 


Fiint, H. T., Grimmett, L. G., Rock Car.inc, F., and Cane, S., 1934, Brit. Med. ¥., 14 April. 
Grimmett, L. G., 1937, Brit. ¥. Radiol., 10, 105; 1942, Ibid., 15, 144. 
Spear, F. G., and Grimmett, L. G., 1933, Brit. #. Radiol., 6, 387. 


JOHN R. LOOFBOUROW 


Dr. LoorsouRow, who had been a Fellow of the Physical Society since 1941, died in 
Boston, Massachusetts, on 22nd January 1951. At the time of his death he was a professor 
-of biophysics at the Massachusetts Institute of Technology, where he was in addition 
Executive Officer of the Biology Department and Chairman of the Faculty. 

Loofbourow, the son of John Wilson and Henrietta Loofbourow, was born in Cincinnati, 
‘Ohio, on 1st November 1902. He was educated at the University of Cincinnati, where he 
obtained the degree of Bachelor of Arts in 1923. He served as an instructor in Physics there 
from 1925 to 1929, and as a research associate for the following six years. In 1936 the 
University of Dayton in Ohio awarded him the honorary degree of Doctor of Science, and 
he became professor of Biology at that institution. From 1935 to 1940 he was a research 
professor at the Institution Divi Thomae in Cincinnati, and in 1940 he was called to the 
Massachusetts Institute of Technology, where he was made an associate professor of bio- 
physics. 

During the war years Dr. Loofbourow was Executive Secretary of the Radar Division of 
the National Research Committee. He served as a special adviser to the Atomic Energy 
Commission in 1947 and 1948, in which. year he was awarded the Presidential Certificate of 
Merit. 

Dr. Loofbourow was a Fellow of the American Association for the Advancement of 
Science, the American Physical Society, the American Academy of Sciences, the New York 
Academy of Sciences, the Optical Society of America, the Beta Theta Pi. He was also a 
member of the Chemical Society, the Faraday Society and the Biochemical Society of 
England, as well as the Physical Society. 

Dr. Loofbourow’s untimely death from a heart complication was tragically followed by 
that of his wife, a practising physician, from the same cause within ten days. ‘They are 
survived by a son, John W. Loofbourow. 

Professor Loofbourow was a frequent contributor to the literature of biophysics, special- 
izing in the field of absorption spectroscopy, particularly in the ultra-violet region. At the 
time of his death he was engaged in important work on the purines and pyrimidines, and was 
engaged in developing new techniques for absorption studies at extremely low temperatures. 
With Professor Richard C. Lord and the undersigned, be was the author of the text Practical 
Spectroscopy published in 1948. 

With his unassuming and friendly manner, breadth of interest, and complete grasp of his 


-subject, Professor Loofbourow had won a host of friends among both colleagues and students. 


As a member of a committee on improving educational methods and as Chairman of the 
Faculty he had contributed much to the academic community. His tolerance and quiet 


competence left a deep impression on all who came in contact with him, and contributed 
-greatly to the present feeling of loss among his colleagues, students, and many friends. 


GEORGE R. HARRISON. 
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EDWARD MALLETT 


Epwarb MALLETT, born in 1888, was educated at Southend High School and University 
College, London. He served an apprenticeship with Messrs. Crompton & Co., Chelmsford, 
and entered the Post Office Engineering Department through a Civil Service competition 
in 1909. Following a short period of training, he wes posted to the West Section of the 
Scotland District and remained there until the outbreak of war in 1914. He was interested. 
in the Territorial Army and as a result, as soon as war was declared, was posted to the 52nd 
Lowland Divisional Signals and: saw service in the Middle East area of hostilities, finally 
leaving the Service with the rank of Captain. 

After the war, Mallett was posted to the Transmission Section of the Engineer-in-Chief’s. 


Office and here his grasp of electrical fundamentals combined with his abounding energy _ 


found outlet in dealing with the problems of telephone transmission, which was then entering” 
into a new phase due to the introduction of the thermionic valve. There isno doubt that 
Mallett, had he remained in the Post Office service, would have had a brilliant future, but a 
friendship made during the war years turned his energies in another direction and he left 
the Post Office service in 1920. It was during this period of service in the Post Office that 
he carried out his researches into the vibrational modes of a circular diaphragm which led 
to the work for which. he was eventually granted his doctorate. 

In September 1921 he joined the academic staff of the City and Guilds College, South. 
Kensington, 2s Assistant Professor of Electrical Engineering, which post he held until 
September 1932. 

During these eleven years he was developing and expanding the ‘light current” 
enginecring section of the department. A considerable expansion was taking place under 
the direction of Professor Fortescue, and Mallett, with his active and fertile mind, backed by 
his considerable experience in the Post Office Engineering Department, found scope for 
developing the instruction and research in electrical communication. 

He had an intuitive engineering sense and considerable academic skill, and these com-- 
bined produced a lecturer of inimitable style and a research worker capable of producing” 
results from the minimum of apparatus. He was extremely popular with both staff and 
students and joined in many of their social activities. 

The title of Reader was conferred on him in 1927 by the University of London, and from. 
that date until his death he remained a member of the Board of Studies in Electrical. 
Engineering. He was for several years Secretary to the Board of Studies, acquiring an 
intimate knowledge of the activities of the Board and of the University, which knowledge 
served them well over the period of his long membership. 


In addition to these activities he still found time to write several papers on a variety of 


subjects of interest to the communications engineer which he presented to the Institution 
of Electrical Engineers and to the Physical Society. 


In July 1926 he was awarded the Degree of D.Sc.(Eng.) as an Internal Student of 


University College, London, the thesis title being Forced Oscillations, Electrical and Mechan- 
ical. He wes frequently consulted by outside bodies and the University on. electrical 
engineering matters, and in 1929 he produced his textbook on Telegraphy and Telephony 
which became, and still is, a valuable foundation for students of electrical communication. 

It was at the height of his ecademic activities that he left the City and Guilds College to 
become Principal of Woolwich Polytechnic, butalthough he became more and more immersed 
in administrative studies, he never lost contact with the College or the University of London. 

As Principal of Woolwich Polytechnic he was intensely interested in the University 
aspect of the Polytechnic but at the same time took a very active interest in the large non- 
university side of Polytechnic activities. He was a keen believer in the National Certificate 
Course for part-time and evening students, and had a place in his heart for the craftsman,. 
ensuring that he was provided with proper shops and tools. It was during his Principalship 
that the Sandwich, Course in Engineering was started, and this eventually developed into the 


Engineering Degree Course of today. This was soon followed by Degree Courses for B.Sc.. 


(Chemistry, Physics, Mathematics), and the full-time School of Art Course. 
The Polytechnic Union came into being during his time. He saw that the union of the 
diversified clubs and societies of the Polytechnic into one corporate whole could be of great 


value to the students and members, and the passing of time has proved him correct. The: 


_<$_ __ 
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success of the Union in all its aspects is a tribute to him: he gave it a free, unfettered 
hand, and held it in high esteem. 

Extensions to the Polytechnic were made during his time, and future developments 
were being planned by him when his fatal illness intervened. 

He was a member of the Institution of Civil Engineers, the Institution of Electrical 
Engineers, the Institute of Physics and the Physical Society (since 1921), and a Fellow of 
University College. 

Edward Mallett, D.Sc.Eng.(Lond.), M.Inst.C.E., M.I.E.E., F.Inst.P., Fellow of 
University College, London, died on the 10th December 1950. C. R. URWIN- 


HANTARO NAGAOKA 


PROFESSOR HANTARO NaGaoka, who had for nearly forty years occupied a worthy 
place among the few Honorary Fellows of the Physical Society, died on 10th December 1950 
at the age of 85. In the English-speaking world his name is probably most familiar in 
connection with ‘ Nagaoka’s formula’ for the calculation of inductance, and indeed his. 
work on the subject has been of great value to workers in the radio field and to physicists 
concerned with precision measurements in inductance. For example, in addition to his 
formulae he published, with S. Sakurai, tables which give, with very little computation, 
values for the self-inductance of single layer coils, and for the mutual inductance of 
coaxial circles to six-figure accuracy. This work, however, is hardly representative of his 
main interest; he was essentially an all-round physicist in the Rayleigh, Kelvin, 
J. J. Thomson tradition, and one who could refer to physics as ‘the Study of Nature’. 

He was born in 1865 and was of honourable lineage. He graduated in physics at the 
Tokyo Imperial University in 1887, a time when, as he recalled later, the library in the . 
physics department occupied but a single shelf six feet wide, and the only periodical 
available was the Philosophical Magazine. Nevertheless he must have immediately 
embarked upon experimental research in magnetism, for papers published by him in 1888, 
1889 and 1890 are discussed in some detail by Ewing in his Magnetic Induction in Iron 
and Other Metals. Among other things Nagaoka showed that in certain circumstances. 
the magnetization of nickel is reversed by torsion, an effect so surprising that both Kelvin 
(then W. Thomson) and Weiss had it confirmed in their own laboratories. Sixty years 
later Nagaoka remarked that the reason why still remained quite unknown. He continued 
his work in Tokyo for six years, took his doctorate in 1893, and then proceeded to 
Germany, where he studied for three years, attending lectures by Helmholtz, Boltzmann, 
Planck and others, and working on magnetostriction. During this time he visited the 
Physikalisch-Technische Reichsanstalt and was received by Lindeck who, having recently 
translated Ewing’s book into German, knew of his work. Lindeck asked him when he had 
left Ewing’s laboratory at Cambridge, and when he learnt that the work had been done 
several years earlier in Japan was sufficiently impressed to show him the laboratory in 
detail. Nagaoka with a pleasant touch of humour said later that he realized then that 
a bit of good work is a much better introduction than a distinguished visiting card. 

Returning to Japan, he was appointed Professor of Physics at his old university, and for 
the next fifty years and more he flourished as an active and versatile leader of physics in 
Japan, occupying positions of increasing responsibility and influence, which in view of 
his own sentiment mentioned above need not be set out in detail here. 

Much of his early work was done on geophysics with Professor A. 'Tanakadate who 
was at that time Director of the Physical Institute of the university; a little later he returned 
to magnetism and with Professor K. Honda investigated the magnetostrictive properties 
of nickel-steels of various nickel contents; they showed that changes of volume, as well 
as of length, occur. The geophysical work included the determination of g at Tokyo, 
comparisons of g at Tokyo and Potsdam, and a systematic study of the elasticity of rocks 
in relation to the propagation of seismic waves through the earth’s crust. His interest in 
geophysics remained with him until the end of his life, but round about 1900 he seems. 
to have been strongly attracted by two widely different scientific developments : atomic 
physics arising from the discoveries of the electron, x-rays and radium on the one hand, 
and wireless communication on the other. His interest in the atom ultimately gave rise 
to a very considerable body of detailed work on spectroscopy, hyperfine structure and 
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the Zeeman effect, carried out with Dr. T. Takamine. His interest in ‘ wireless ’ seems 
to have led to his work on formulae for inductance calculation, already mentioned, and 
to a consideration of the effect of the ionosphere, familiar to him from its association with 
terrestrial magnetism which he had investigated earlier, on radio communication, 
He foresaw that the solar eclipse of 1914 would provide an opportunity of testing the 
Kennelly—Heaviside theory, if the ionized layer arises from the photoelectric effect of 
solar radiation. He made a calculation of the effect that should be observable on suitable 
wireless signals. Unfortunately the first world war prevented the making of the 
observations, but the detailed observations made in the U.S.A. at the eclipse of 1925 
showed that the idea was sound although, in the light of later information, Nagaoka’s 
original calculation was of course only the first rough approximation. 

It is estimated that Nagaoka published in all some 300 papers, and it will be realized 
from what has been said that they are so diverse in character that they do not lend 
themselves to brief effective description. What does stand out very clearly is his life-long 
devotion to physical science : we see him as mathematical physicist, an experimenter—he 
devised a magnetograph, a silica pendulum for gravity determinations and spectroscopic 
devices—a patient observer—he himself is said to have made fifty determinations in the 
gravity survey that he organized—and an original thinker. On this last point it must 
be noted that he had the honour of being the first to suggest a planetary structure for the 
atom. Rutherford in Radioactive Substances and their Radiations wrote “ This type of 
‘ Saturnian ’ atom was first suggested in 1904 by Nagaoka, who investigated mathematically 
the stability of an atom composed of a central attracting mass surrounded by rings of 
rotating electrons. He showed that such a system was stable if the attractive force was ~ 
large and the electrons in rapid motion’’. When it is remembered that it was in this 
year also that J. J. Thomson was still investigating the model consisting of electrons 
inside a positively charged sphere, the originality of the conception can be better 
appreciated. It was spectroscopy rather than radioactivity that was in Nagaoka’s mind; 
his aim was to account for the discrete spectral lines characteristic of the atom, and he 
thought that they might be associated with vibrations of the electrons in and perpendicular 
to the plane of the ring. We now know that the vital key to the problem was still missing, 
and that Bohr was to find it some years later, but Nagaoka’s contribution commands our 
respect. 

Nagaoka came to occupy an eminent position in the scientific life of Japan. He was 
elected a member of the Imperial Academy in 1905, became Director of the Institute of 
Physical and Chemical Research in Tokyo (now the Scientific and Research Institute) 
in 1926 and was among the first recipients of ‘ Bunkwa Kunsho’, a medal of great 
distinction, in 1937. He was President of the Imperial Academy in 1939 and 1948. 
He visited England on two or three occasions, was awarded an honorary degree at Cambridge, 
and read papers in London before the Physical Society. He was elected an Honorary 
Fellow of the Society in 1912. That this tribute to an outstanding fellow worker in the 
new school of physics then rapidly rising in Japan has survived the shock of two world wars 
may serve as a welcome reminder that physics still stands as a co-operative human enterprise. 

L. H. 


J 


FRANK BENNETT YOUNG 


Ir is recorded with regret that Frank Bennett Young died on 29th January 1951. He was a 
Fellow of the Physical Society from 1916 to the time of his death, and a founder Fellow (1920) 
and life member of the Institute of Physics. 

He was born in 1877 in Bristol where he received his early education. In 1896 he went 
to the Borough Road College, Isleworth, where he was awarded a scholarship enabling him 
to spend a year at Jena University. In 1899 he obtained the B.A. degree of London Univer- 
sity. At this time he became interested in science from a philosophical point of view and 
began to study physics, returning to Bristol to teach physics in a Higher Grade School. 
During this period he attended evening classes at University College, Bristol, before the 
foundation of the University. It was there he met Professor Chattock who had an important 
influence on his career. His enthusiasm for research in physics increased and he spent all 
his spare time (evenings, week-ends and holidays) on research under Chattock’s guidance. 
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His published papers at this time included ‘“‘ The field at the point of an electrified needle ’” 
(Phil. Mag., 1907, 13, 542) and “‘ Critical phenomena of ether ”’ (a study of the opalescent 
phenomena near the critical point) (Phil. Mag., 1910, 20, 793). On the foundation of the 
University of Bristol he took advantage of one of its early ordinances which permitted him 
in 1910 to present these publications for the degree of B.Sc. by research. In 1915, soon after 
the outbreak of the 1914-18 war (Prof.) A. M. Tyndall obtained for him the job of localizing 
shrapnel bullets in the body by x-rays at the Bristol Royal Infirmary, a problem quite beyond 
the medical man of that day. This was the first time a physicist had entered its doors to 
give service. 

In 1916 he joined the Admiralty, for anti-submarine research under the Board of Invention 
and Research (W. H. Bragg and Rutherford were on the controlling committee). He worked 
in Admiralty research establishments at Aberdour (Fife), Parkeston Quay (Harwich), 
Shandon (Gareloch) and, after the war, at Teddington (Middlesex). In the earlier days he 
was mainly concerned in research work on hydrophones and in underwater acoustics. 
Some of this work was published after the war, viz. ‘“‘ On light body hydrophones and the 
directional properties of microphones ”’ (Proc. Roy. Soc. A, 1921, 100, 252) and ‘“‘ On the 
acoustic disturbance produced by small bodies on plane waves transmitted through water, 
with special reference to the single-plate direction finder’ (Proc. Roy. Soc. A, 1921, 100, 
261). ‘This had an important bearing on the design of directional hydrophones. At a later 
stage of the 1914-18 war he carried out important researches on. the generation of electric 
currents in the sea (a) due to the presence of steel ships, and (6) due to the tidal flow of a 
conducting medium (the sea) in the earth’s field. Some of the experiments and conclusions. 
reached in the latter were published in Phil. Mag., 1920, 49, 149, in an article entitled 
““ Electrical disturbances due to tides and waves’’ by F. B. Young, H. Gerrard and W. 
Jevons. At this time he was also interested in the development of a new method of locating 
faults in submarine cables—using alternating currents and silver-chloride electrode search 
gear. This method has since been extensively used by cable companies and by the G.P.O- 
Soon after the end of the 1914-18 war he obtained the degree of D.Sc. (Bristol 1919). 
On the removal of the Admiralty research establishment from Shandon (Gareloch) to 
Teddington he worked on course plotters for ships and developed what is now known as the 
A.R.L. (Admiralty Research Laboratory) Course Plotter. Another line of research which 
he directed related to the preservation of the rubber insulation of cables in ships. _ It will 
be understood that much of his research work for the Admiralty remained unpublished on 
grounds of security. In July 1934 he became Superintendent of the Admiralty Research 
Laboratory, Teddington, a post which he filled until his retirement on reaching the age limit 
in October 1937. Whilst in this post he was awarded the O.B.E. 

In 1939 he returned to Admiralty work in Portsmouth soon after the outbreak of the 
1939-45 war, taking charge of an important section of research on the degaussing of ships. 
as a protection against enemy magnetic mines. He finally retired from Admiralty service 
in 1941, taking on various unpaid wartime duties in Bristol. 

He eventually settled in Portsmouth after the war where he revived his old-time interests 
in music and in sketching and water colour painting, He was of a retiring disposition, very 
reliable, thorough and conscientious in all his work, and generous in giving more than their 
share of credit to others who worked with him. He will be missed by those who, like 
myself, knew him well. A. B. WOOD. 


ARNOLD SOMMERFELD 


On 26th April 1951, Arnold Sommerfeld, Honorary Fellow of the Physical Society, died 
from the consequences of a street accident at the age of eighty-two. 

He was born in 1868 at Konigsberg, East Prussia, and attended school and university in 
his native town. In 1893 he went to Gottingen and came under the influence of Felix 
Klein. Here he had athorough mathematical training and acquired the skill in the theory of 
differential equations, analytical functions and complex integration which he later applied 
with great success to physical problems. The direct fruit of bis co-operation with Klein was. 
their famous book on the spinning top (1895-1910). 

In 1895 Sommerfeld became lecturer in G6ttingen, in 1897 professor at the Mining 
College in Clausthal and in 1900 he was appointed to a chair of mechanics at the Technical 
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College in Aachen. From 1906 until his death he was professor of theoretical physics in 
Munich. He declined several attractive positions, as for instance the chair of Planck in 
Berlin. f 

Sommerfeld’s scientific work covers a vast field of which only some major points can 
here be mentioned. His contributions to classical physics were in hydrodynamics, by work 
on the stability of laminar motion, and in optics, which he enriched by ingenious rigorous 
solutions of some diffraction problems. But soon he directed his attention to the new ideas 
of relativity and quantum mechanics. Combining these two subjects he developed Niels 
Bohr’s theory of atomic spectra by taking into account the relativistic variability of mass, 
and obtained in 1915 his celebrated fine structure formula for hydrogen-like spectra. 
Sommerfeld’s fine structure constant dominates not only these phenomena but the whole of 
quantum electrodynamics. In 1917 be hegan to write his great book on the structure of 
atoms and spectral lines which has so immensely contributed to the understanding and 
general acceptance of quantum theory. The advent of quantum mechanics led him first to 
write a short complementary volume which was later extended to a comprehensive repre- 
sentation containing a wealth of information. The study of wave mechanics in combination 
with Pauli’s principle led Sommerfeld to another great success, its application to the electron 
theory of metals. 

Sommerfeld was a brilliant teacher and attracted pupils from all over the world. When 
in 1938 he retired from his Chair, he began to prepare the publication of his lectures, of which 
all but one volume have appeared, and it is hoped that this last one will also be published. 
His style and power of representation is such that his books read like novels. But still — 
greater was the charm of his personality which won him the devotion of his pupils, amongst 
whom were Pauli, Heisenberg, Brillouin, F. London and many other outstanding men. 

Sommerfeld was a member of many academies and learned societies, amongst them the 
Royal Society. He received the Max Planck Medal of the German Physical Society and 
many other honours. He was at one time Karl Schurz Professor at the University of 
Wisconsin and visiting professor at the Californian Institute of Technology and at several 
other American universities. MAX BORN. 


STUART ARTHUR FRANK WHITE 


STuAaRT ARTHUR FRANK WHITE, who died on 12th January 1951, aged 80 years, was 
elected a Fellow of the Physical Society in 1895. He was for many years Professor of 
Mathematics at King’s College, London, Dean of the Faculty of Science and latterly 
Vice-Principal and Fellow of the College. He retired in 1933 after occupying his Chair for 
thirty years. 

White’s education began at Portsmouth Grammar School from which he proceeded 
(with a scholarship or exhibition, I believe) to Wadham College, Oxford, at Michaelmas 
1887. In his student days he seems to have been as much interested in physics and 
astronomy as in mathematics, and before going to King’s College he held an appointment 
at Oxford under Professor Clifton. In 1893 he was appointed Demonstrator of Natural 
Philosophy—strangely enough, on the recommendation of the Engineering Committee— 
under Professor Grylls Adams, and in 1903 he succeeded W. H. Hudson in the Chair of 
mathematics. ‘The appointment to this Chair of the demonstrator of natural philosophy, 
rather than of a mathematician vom Fach, appears to have been due to the desire of the 
Council to infuse a more vivid and practical spirit into the teaching of mathematics to the 
students of such subjects as engineering, physics, etc., whose interest lay less in mathematics 
for its own sake than in its practical applications. For such students White was ideal. 
He was a most able teacher, lucid and always interesting.. He was moreover whole-heartedly 
devoted to his teaching duties and to the progress and welfare of his students, whose affection 
and esteem he undoubtedly won. 

White was a keen golfer and skier and for a long period spent part of every Christmas 
vacation in Switzerland. He was intensely devoted to the college he served and his wise 
counsel and guidance were greatly valued, not only by the Professorial Board of the College, 
but also by his mathematical colleagues in the university. WM. WILSON. 
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Miss JESSIE S. CHREE 
Fe.L.Lows of the Physical Society will be sorry to learn of the death of Miss Jessie S. 
Chree, sister of the late Dr. Charles Chree, F.R.S. 


Dr. Chree was at one time President of the Physical Society, and the Society’s 
Charles Chree Medal and Prize was instituted in memory of him through the generosity 
of Miss Jessie Chree. 


ALLAN FERGUSON 


WE record with regret the death on 9th November 1951 of Dr. Allan Ferguson who was 


the Papers Secretary of the Physical Society from 1928 to 1938 and President from 1938 
to 1941. 


An appreciation will be published with next year’s Obituary Notices. 
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ABSTRACTS FOR SECTION A 
The Temperature Effect of Extensive Air Showers, by A. L. Hopson. 


ABSTRACT. Extensive air showers have been recorded at sea level over a period of 
nine months. The data show a fairly good correlation between the shower rate and the 
air temperature at about one ‘cascade unit above the apparatus, e.g. at the 950 mb. level. A. 
temperature coefficient of —(0:38+0-11)% per °c. at this level is indicated. This may be 
explained by the dependence of the lateral spread of an extensive shower on the air density. 


The Disintegration of Light Nuclec by Meson Capture, by A. C. Ciark and 
S. N. RUDDLESDEN. 


ABSTRACT. The disintegration by 7~-meson capture of the elements helium and 
carbon is discussed. ‘Two types of meson theory, scalar and pseudo-scalar, are considered. 
The relative probabilities of the various modes of disintegration of helium are calculated 
and the energy spectra of the most likely one, viz. proton emission, are shown. The 
disintegration of carbon into two a-particles and one singly charged particle is considered, 
using a simple o-particle model. The energy spectra and angular correlations of the 
a-particles obtained by a direct transition and by a transition via an excited state of *Be 


are compared. 


Range-Energy Relation for Protons in Ilford C2 Dry Emulsion, by F. A. EL-BEDEWI1.. 


ABSTRACT. ‘The range-energy relationship for protons of energies up to 12 Mev. in 
Ilford C2 dry emulsion has been determined with an accuracy better than 0-5% in energy. 
The particles were produced by bombarding a thin beryllium foil with deuterons having 
an energy of 8mMev. The energy of the incident deuterons was obtained by a method 
which involves the measurement of the two angles of emission at which the two proton 
groups of longest range from the reaction °Be(d, p)!°Be have the same energy. 


The Differential Cross Section for High Energy Nucleon—Nucleon Collisions and 
the Mean Square Angle of Scatter, by H. S. GREEN AND H. MEssEL. 


ABSTRACT. In an attempt to account theoretically for the lateral and angular develop-. 
ment of the penetrating particles in high energy cosmic-ray showers, we have developed a 
method for deriving the differential from the total cross section, in the laboratory frame of 
reference. An integral equation is set up for the differential cross section, and solved 
assuming that the differential cross section in the centre-of-mass frame of reference is of the 
form R+,S’ cos? 6, where R and S ‘are functions of the incident and scattered energies only. 
Using a form of the total cross section derived from the analysis of high-energy cosmic-ray 
experiments, the mean square angle of scatter in a nucleon-nucleon collision has been 
calculated and found to be inversely proportional to the energy of the incident particle. 
The results are compared with experiments. 


Nuclear Spin and Magnetic Moment of Rhodium ‘3Rh, by H. Kuxyn and 
G. K. WoopGateE. 


ABSTRACT. The radiation from a hollow-cathode tube containing rhodium was passed 
four times through an atomic beam of rhodium. Absorption was found in eight lines, 
four of which showed hyperfine structure doublets of separation 0-023 cm. The splitting 
could be ascribed to the ground state 4d® 5s 4Fy/5. It was concluded with certainty that 
the spin of the rhodium nucleus 1{3Rh has the value J=4. The relative intensities of the 
components show that the magnetic moment is negative. 'The observed splitting leads. 
to an estimated value of “= —0-10 nuclear magneton. The values of I and m indicate 
a configuration (5£9/.)°3p,/. for the last seven protons of the rhodium nucleus in the 
Goeppert-Mayer model. 
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Slow Neutron Scattering by Ferromagnetic Crystals, by R. G. Mooruouse. 


ABSTRACT. Consideration is made of the predictions of the Bloch spin wave model 
within its region of approximate validity, for the magnetic scattering of very slow neutrons 
by a ferromagnetic. It is found that observable peaks of intensity in the angular distribution. 
of scattered neutrons occur due to spin wave absorption. Neutrons are also scattered 
with phonon absorption and by an analogous calculation it can be shown that similar 
peaks occur due to this vibrational scattering. A possible method is suggested of differ- 
entiating between the two kinds of scattering in order to provide a check of the Bloch model. 


Antiferromagnetism by the Bethe Method, by J. M. Zman. 


ABSTRACT. The Bethe method for order—disorder is applied to a lattice of magnetic 
dipoles with antiferromagnetic nearest-neighbour interactions. Magnetization and entropy 
are calculated for all values of temperature and magnetic field, and an anti-parallel ordered 
phase is found below a definite transition curve. The computed curves agree qualitatively 
with those observed experimentally with certain paramagnetic salts at low temperatures 


Asymmetrical Electron Diffraction Effects from Single Crystals of Silver, by 
D. W. PAsHLEY. 


ABSTRACT. Silver crystal surfaces which had been prepared by either electrolytic 
etching or electrolytic polishing were studied by means of the electron diffraction reflection 
technique. It was found that fine structure effects, in the form of split spots or streaks. 
through spots, gave information concerning the topographies of the various surfaces of the 
crystals. ‘The fine structure effects showed certain asymmetries which do not immediately 
follow from the Laue diffraction theory. The possible causes of these asymmetries are 
discussed, and it is concluded that they were probably due to distortion of the silver crystals, 
but that some might have also been due to refraction effects. 


Field Dependence of the Dielectric Constant, by J. J. O’ Dwyer. 


ABSTRACT. A very general calculation has been made by Frohlich giving the dielectric 
constant of a material in terms of the spontaneous fluctuations of the electric moment of the 
material which would occur ina field-free region. This calculation is now extended one step 
further so as to include the first term in the field dependence of the dielectric constant. By 
way of illustration the general theory is applied to Kirkwood’s model for a dipolar liquid. 
(In order to preserve continuity of argument it has been found necessary to reproduce much 


of Fréhlich’s work.) 
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